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Diffuse Scattering of X-Rays from Sylvine. II 


G. G. HARvEy,* Ryerson Physical Laboratory, University of Chicago 
(Received February 16, 1933) 


The absolute intensity of the diffuse scattering of x-rays 
of wave-length 0.71A from sylvine has been measured for a 
range of angles from 5° to 110°. Homogeneous radiation 
obtained by the use of balanced filters was used. Woo’s 
formula taking account of the incoherence of part of the 


scattering was used to calculate the average atomic 
structure factors of the ions Kt and Cl-. The results are in 
good agreement with Wollan’s values for argon as well as 
theoretical values calculated from a Schridinger charge 
distribution by Hartree’s method. 





I. INTRODUCTION 


T has been shown! ? that the diffuse scattering 
from a simple cubic crystal consisting of 
atoms of one kind should be given classically by 


Sctass. =1+(Z—1)f?/Z?— F?/Z, (1) 


where Sciass. is the scattering per electron in 
terms of the Thomson® value in the same direc- 
tion, Z the atomic number, F the atomic struc- 
ture factor including the effect of thermal agita- 
tion and f’ a quantity related to! the true atomic 
structure factor f. Woo‘ has independently de- 
duced a similar formula, taking account of the 
fact that part of the scattering is incoherent. 
Woo’s formula may be written 


S=(f"—F*)/Z+(1—-f"/Z?)(1+<a vers ¢)’, (2) 


where ¢ is the angle of scattering and a=h/men. 
The first term in Eq. (2) represents the coherent 


* National Research Fellow. 

1G, E. M. Jauncey, Phys. Rev. 37, 1193 (1931). 

*G. E. M. Jauncey and G. G. Harvey, Phys. Rev. 37, 
1203 (1931). 

*J. J. Thomson, Conduction of Electricity through Gases, 
2nd edition p. 325. 

*Y. H. Woo, Phys. Rev. 38, 6 (1931). 


and the second the incoherent part of the scatter- 
ing. 

Several investigators’: °:7»*.*® have reported 
measurements on the diffuse scattering of x-rays 
from crystals but in no case has homogeneous 
radiation been used. Although previous measure- 
ments on sylvine’ gave values of the structure 
factors f’ in good agreement with Wollan’s 
values” for argon it was thought that this agree- 
ment was somewhat fortuitous in that no cor- 
rection was made for incoherent radiation (Woo’s 
paper‘ had not then appeared) and the wave- 
length range, though apparently narrow as 
measured by absorption, must have been rather 
broad. Further, it does not seem to the writer 
that an average wave-length as measured by ab- 
sorption is the correct average to use for the 


5G. E. M. Jauncey, Phys. Rev. 20, 405 (1922). 

®°G. E. M. Jauncey and H. L. May, Phys. Rev. 23, 128 
(1924). 

7G. G. Harvey, Phys. Rev. 38, 593 (1931). The values of 
S¢/p on p. 601 of this paper must all be multiplied by 0.024. 

8G. E. M. Jauncey and P. S. Williams, Phys. Rev. 41, 
127 (1932). 

*G, E. M. Jauncey and Paul Ehrenfest II, Phys. Rev. 
42, 907 (1932). 

1° E, O. Wollan, Phys. Rev. 37, 862 (1931). 
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process of scattering. For if \. is the average 
wave-length as measured by absorption and \, 
the correct value to use for scattering and if J(A) 
is the spectral distribution then 


em Hida) f I(\)dd\ = f T(d)e-#™ "dd 
do 


do 


for some particular range of values of ¢ and 


I(d)S(d)d, 


do 


SQ.) f “I(\)d\= 
do 


where (A) is the absorption coefficient of the sub- 
stance used in determining a, 40 the minimum 
wave-length and S is given by Eq. (2). It is ob- 
vious that in general \,+),, though of course 
they may not differ greatly if the spectrum is not 
too extended. It is to be noted, however, that \, 
will also depend on ¢. Recently in the case of 
sodium fluoride Jauncey and Williams* have 
found it necessary to take account of the fact that 
part of the scattered radiation is incoherent and 
although this effect is smaller for sylvine it is 
not negligible. It was hence thought desirable to 
repeat the measurements on sylvine using homo- 
geneous radiation obtained by the use of a 
balanced filter." 


II. APPARATUS AND PROCEDURE 


A diagram of the apparatus used is shown in 
Fig. 1. T is a molybdenum target x-ray tube im- 
mersed in transformer oil in a lead box A. The 
tube was operated at 30 m.a. at a potential of 
about 42.5 kv peak with full wave kenotron 


Cf) ‘es 
we T 
Iw L! 


Fic. 1. Diagram of apparatus. 















1 P, A. Ross, Phys. Rev. 28, 425 (1926). 
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rectification. The x-rays left the box through a 
small celluloid window after passing through 
about 2 mm of oil. The primary beam was defined 
by the “‘slits”’ s; and s2 consisting of small circular 
holes. s; served to eliminate scattering from the 
slit 52 and was slightly larger than the x-ray beam, 
The angular width of the beam was about one 
degree. The balanced ZrO,-SrO filter was placed 
in front of s; and the radiation scattered by the 
crystal C was measured by an ionization chamber 
and a Compton electrometer operating at a 
sensitivity of 5000 mm/volt. The ionization 
chamber was built of copper and fitted with thin 
celluloid windows the absorption in which was 
too small to be measured. All insulation was of 
amber. The chamber was 7.5 cm long and filled 
with krypton at a pressure of 57 cm of mercury 
N.T.P. This results in an absorption of 85 percent 
for the Mo Ka lines. The use of such a small 
chamber not only greatly reduces the residual 
ionization due to cosmic radiation and a-par- 
ticles but results in a rather low capacity elec- 
trometer system with a consequent high charge 
sensitivity. Further details on the construction of 
such ionization chambers have been given by 
Compton elsewhere.” 

The transmission by the filters is shown in Fig. 
2, the unshaded portion representing the part 
used. It is not to be supposed that this represents 
the energy distribution in the spectrum as no 
correction for absorption or second order reflec- 
tion was made and a rather broad primary beam 
was used. This accounts for the peak at 0.8A and 
the very broad Ka lines, though of course does 
not affect the correctness of the balance. 

The procedure was the same as in previous ex- 
periments except that readings were taken with 
first one filter and then the other and the differ- 
ence of the currents taken as a measure of the 
scattering for the Mo Ka lines alone. Due to the 
fact that it is the difference of the readings that is 
required, it is necessary to take a large number of 
readings and in fact several thousand were taken. 
With the chamber set at an angle ¢ the ionization 
current was measured for a series of values of 8, 
the angle between the primary beam and the nor- 
mal to the crystal face. At @= /2 the Laue spot 


was not entirely absent as this would necessitate 


2 A. H, Compton, Rev. Sci. Inst. 2, 365 (1931). 
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Fic. 2. Transmission by filters. 


a more perfect balancing of the filters than is 
possible since the reflected and scattered inten- 
sities are of entirely different orders of magnitude. 
The Laue spot was never very intense, however 
(except of course for \=0.71A). The scattering 
from the crystal at an angle @ was then compared 
with that from a paraffin block at 90°. The ratio 
of the intensity transmitted by the crystal when 
set at 02=¢/2 to that transmitted by the paraffin 
at @=45° was also measured. Because of the 
rather high mass absorption coefficient of sylvine 
(13.7 cm* g~') it was necessary to use a very thin 
crystal, the surface density being 0.0843 g cm~. 
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The surface density of the paraffin was 0.263 g 
cm~*, The range of angles used extended from 5° 
to 110° and at the small angles it was necessary 
to correct for the height of the slit in front of the 
ionization chamber as it subtended an angle of 
about 5°. 

It has been shown" that Crowther’s formula'‘ 
when modified to take account of the fact that 
part of the scattered radiation is incoherent be- 
cause of the Compton effect, becomes 


I,=1,'t(sitseT) sec o/2, (3) 


where J, is the scattered intensity per unit solid 
angle in the direction ¢, J,’ the transmitted in- 
tensity when the scattering block is set at 6= ¢/2, 
t the thickness of the scattering block, s; and Ss 
the spatial scattering coefficients for the unmodi- 
fied and modified rays, respectively, and 


T= (1—e- wt see $12) /(u! — u)t sec @/2, (4) 


uw’ and uw being the linear absorption coefficients 
for the modified and unmodified rays. In Eq. (4) 
we have neglected the difference in the absorp- 
tion of the modified and unmodified rays in the 
air and ionization chamber window. Since the 
absorption in the ionization chamber is not com- 
plete account must be taken of the fact that the 
incoherent portion is more completely absorbed 
than the coherent. As the derivation of the re- 
quired equation has been given by Jauncey and 
Williams® it will not be repeated here. Jauncey 
and Williams assumed all the radiation scattered 
by paraffin at 90° to be incoherent, since they 
used a rather short wave-length (0.39A). For 
0.71A, however, the percentage of coherent 
scattering becomes appreciable and this intro- 
duces a slight modification. Their equation, thus 
modified, becomes 





Ko(f" — F*) 


1-f"/Z2 
—+K,T 














ibe "(ite vers ¢)> 


— (5) (<) (<)(-) (5) 
1+cos? @ ZF eo\W p\ Poof \C 





where Ky, is the percentage absorption in the 
ionization chamber for radiation of wave-length 
0.71+0.024 vers ¢, Cy the ionization current for 
scattering from the crystal at an angle ¢, Ps the 
same for paraffin at 90°, C and P the currents 
produced by the rays transmitted by the crystal 
and paraffin, (pt)c¢ the surface density of the 


crystal, (pt)p the same for paraffin, W and Z the 
atomic weights and atomic numbers, A and B 
the percentages of coherent and incoherent scat- 
tering from paraffin at 90°, and the other quanti- 

13 G. E. M. Jauncey and O. K. DeFoe, Phil. Mag. 1, 711 


(1926). 
4 J. A. Crowther, Proc. Roy. Soc. A86, 478 (1912). 
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ties have been previously defined. K, only varies 
from 0.85 to 0.88 from 0° to 110°. 

To calculate the coherent scattering from paraf- 
fin Wentzel’s formula’ in the form given by 
Compton" was used in a previous” paper giving 


Teon. = 41 e/(1+0.352 x?)4, (6) 


where J, is the Thomson value of the scattering 
from a single electron and x=(sin ¢/2)/X, as the 
coherent intensity scattered by each carbon 
atom in paraffin, it being assumed that only the K 
electrons scatter coherently. Likewise for the in- 
coherent intensity from the K electrons of a 
carbon atom we obtain 


Tineon. = 2 [1 —1/(1+0.352x*)*]. (7) 


Taking the formula for paraffin as C,He2,42 and 
assuming all the scattering from hydrogen to be 
incoherent we obtain for ¢=90° 


Teon./Ttotai = 1.2n/(7.4n +2) =0.16 (8) 


for all values of m above 20. Thus in Eq. (5), 
A=0.16 and B=0.84. Inserting the exp¢rimental 
values of Cy/P and P/C in Eq. (5) and using 
values of F as given by James and Brindley" we 
may calculate values of f’. The structure factors 
so calculated are given in Table I. The values in 


TABLE I. Structure factors. 











ci) x F f Sclass. 
5° 0.061 17.2 (17.2) (0.06) 
10 0.123 14.3 15.15 1.80 
15 0.184 11.8 12.95 2.07 
20 0.245 9.9 11.61 2.62 
25 0.306 7.85 10.00 2.82 
30 0.365 6.45 9.03 2.97 
35 0.424 5.48 8.15 2.82 
40 0.482 4.57 7.61 2.87 
50 0.596 3.40 6.65 2.67 
60 0.704 2.45 5.95 2.53 
70 0.808 1.73 5.6 2.47 
80 0.905 1.12 4.7 2.09 
90 0.995 0.75 4.3 1.98 
100 1.080 0.50 (3.6) (1.67) 
110 1.154 0.35 (3.25) (1.55) 








18 G. Wentzel, Zeits. f. Physik 43, 1, 779 (1927). 

16 A. H. Compton, Phys. Rev. 35, 925 (1930). 

7G. E. M. Jauncey and G. G. Harvey, Phys. Rev. 37, 
698 (1931). In this paper the factor 2 was omitted in the 
exponential occurring in u(r), Eq. (11), and the exponent 4 
in the denominator of Eq. (12). The correct form was used 
in the calculations, however. 

18 R. W. James and G. W. Brindley, Proc. Roy. Soc. 
A121, 155 (1928). 


HARVEY 


parentheses are rather doubtful because of the 
extremely weak intensity. These values of f’ are 
also plotted in Fig. 3. The solid curve represents 
the average of fx+ and f,,- as calculated by James 
and Brindley” from a Schrédinger charge distri- 
bution by the method of Hartree. The crosses 
represent Wollan’s values ' for argon calculated 
from Compton’s formula for the scattering from 
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Fic. 3. Structure factors. Crosses, argon (Wollan); circles, 
sylvine; solid curve, 3}(K*+Cl-)-wave mechanics. 


a monatomic gas'* and the circles represent the 
present values for sylvine. Although it was not 
possible with the wave-length used here to go out 
to as large values of (sin ¢/2)/A as previously’ 
there is an indication that the present values give 
better agreement with theory at large angles. It 
may be mentioned that the f values for argon as 
given by James and Brindley” are slightly larger 
than the average for K+ and Cl- at small values 
of x but slightly smaller for large x. It is seen 
from Fig. 3 that at small angles Wollan’s values 
for argon are above those for sylvine, though the 


19R, W. James and G. W. Brindley, Phil. Mag. 12, 81 
(1931). 
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Fic. 4. Comparison with scattering from argon. Crosses 
= Sgas (Wollan); circles = (S+ F?/Z)crystal- 


difference is perhaps within the combined experi- 
mental errors of the two sets of measurements. 
At large angles the difference is in the same direc- 
tion but the intensity at large angles is extremely 
weak and this difference is less than the probable 
error of the sylvine measurements alone. 

It has been pointed out” that the scattering 
from a crystal, reflection from a crystal, and 


*0G. E. M. Jauncey and G. G. Harvey, Phys. Rev. 38, 
1071 (1931). 
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scattering from a monatomic gas having the 
same atomic number as the crystal are related by 


the equation 


Seas = (S+ F?/Z) cryetal- (9) 


This relation holds whether the S values are 
classical or modified to take account of the Comp- 
ton effect but as Wollan’s values are tabulated as 
Seiass. We Shall use this form. The values of S in 
the last column of Table I were calculated from 
Eq. (1) with the f’ and F values given there. In 
Fig. 4 values of Syass + F?/Z for sylvine are indi- 
cated by circles and Wollan’s values of Sos. for 
argon by crosses. It is seen that the agreement is 
very good on the whole but, again, at small 
angles there is a small difference. 


III. CoNcCLUSION 


Data on the diffuse scattering from sylvine 
when combined with results on reflection yield 
values of the structure factors which within the 
limits of experimental error agree with those cal- 
culated from a Hartree model. The difference 
between the structure factors for argon and syl- 
vine, though hardly significant, is in general in 
the right direction. 

The primary purpose of this investigation was 
to determine the effect of temperature on the in- 
tensity of the diffuse scattering. Experiments at 
the temperature of liquid air with monochromatic 
radiation are now in progress and will be reported 
later. 

In conclusion the author wishes to express his 
appreciation to Professor Arthur H. Compton for 
his interest in the work and for the use of the 
Ryerson Physical Laboratory. 
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Absolute X-Ray Reflectivities of Single Crystals of Calcite, Rocksalt, Rochelle Salt 
and Barite 


PAUL KIRKPATRICK AND P. A. Ross, Stanford University 
(Received February 4, 1933) 


With a two-crystal spectrometer absolute x-ray re- 
flectivity measurements have been performed upon four 
different crystals. In the case of a large calcite crystal of 
uniform appearance, considered a good spectrometer 
crystal, with Ka radiation of silver it is found that the 
theoretical values of percent reflection, integrated re- 
flection, and rocking curve width calculated by Allison on 
the basis of the Darwin-Prins theory of reflection by a 
perfect absorbing crystal are not attained, the experimental 
and theoretical quantities being related by a factor not far 
from two in each case. Rocksalt crystals from a large 
exhibition specimen from Russia were found to be much 
more perfect in structure than any rocksalt specimens 
heretofore reported in the literature of x-ray spectrometry. 


Rocksalt reflectivity data for Ag Ka are given. Artificial 
Rochelle salt crystals are found to possess rocking curve 
widths as narrow as those of high-grade calcite, indicating a 
very uniform structure; at the same time the reflections 
are relatively weak. Rochelle salt reflectivity data for 
Ag Ke are given. Percent reflection, rocking curve width 
and integrated reflection of barite were determined for 
eleven wave-lengths, and scattering factors were calculated. 
The variation of the scattering factor with wave-length in 
the neighborhood of the barium K limit is similar to the 
variation previously found by Wyckoff and by Glocker and 
Schafer for lighter elements except that Wyckoff’s maxi- 
mum does not appear. 





N relatively few cases have the x-ray reflecting 

characteristics of crystals received the direct 
attention of investigators. Calcite has been 
rather adequately studied on account of its great 
interest as a spectrometer crystal. Existing data 
upon rocksalt may be traced partly to a similar 
interest and partly to the importance of this 
crystal as a standard reflector for use in the 
analyses of other crystals. In the performance of 
such analyses spectrometer observations of in- 
tensities of reflection have of course been made on 
many crystals, but while these observations have 
been sufficient for their ends they are often of 
diminished permanent value in themselves either 
because of the slight attention paid to absolute as 
contrasted with relative values, or because the 
reflection data themselves do not attain publica- 
tion. Since the integrated reflection is the only 
reflection characteristic of use to the crystal 
analyzer, his work at best provides no informa- 
tion about the angular distribution of reflecting 
power, nor can this characteristic be deduced in 
the general case from a knowledge of the crystal 
structure, however complete. 

The reflectivity observations presented below 
were obtained with a two-crystal spectrometer 


which has been previously described.' The x-ray 
tube, a dismountable type with changeable 
targets, was operated on constant potential direct 
current and continuously evacuated by a strong 
pumping system. The reflectivity determinations 
discussed below were all conducted according to 
a uniform procedure which in its general aspects 
is familiar through the writings of Davis,’ Alli- 
son’ and others. Having selected the crystal to 
be examined, an entirely similar one in material, 
perfection, and treatment was mounted in the A 
position. The Seemann slit of the A crystal was 
adjusted to an opening of about half a millimeter. 
A narrow self-centering slit was mounted on the 
B crystal table, the crystal and its Seemann slit 
being absent, and the A crystal was adjusted in 
its holding device by screws provided for the pur- 
pose until radiation reflected from it was received 
centrally through the slit on the B table. This 
condition was sometimes observed by using the 
spectrometer ion chamber, but the adjustment 


1P. A. Ross, R.S.I. 3, 253 (1932). 

?B. Davis and W. M. Stempel, Phys. Rev. 17, 608 
(1921); 19, 504 (1922). B. Davis and H. M. Terrill, Phil. 
Mag. 45, 463 (1923). 

3S. K. Allison and J. H. Williams, Phys. Rev. 35, 1476 
(1930). 
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ABSOLUTE X-RAY REFLECTIVITIES 


could be made more rapidly and quite as accu- 
rately by employing a fluoroscope. A vertical slit 
about a millimeter wide and a horizontal slit 
several millimeters high were installed between A 
and B, to limit the beam incident upon B to a 
roughly central area of the crystal. The adjusting 
slit was removed and the B crystal placed in the 
parallel position with its Seemann slit opened so 
wide as to be inoperative. The twice-reflected 
beam was found with the chamber, and adjust- 
ments, often laborious, to insure parallelism of 
the mean reflecting planes of the two crystals 
were carried out. 

Parallel rocking curves were obtained by turn- 
ing crystal B, step by step, through the reflecting 
position and observing the power of the twice- 
reflected beam at each step. Tube potentials were 
maintained at such values as to avoid the produc- 
tion of undesired high-order reflections. Before 
and/or after determining the rocking curve, 
crystal B was removed from its table and the 
power of the once-reflected beam observed. 

From these data were derived, in the case of 
each crystal and wave-length, the integrated re- 
flection R, the percent reflection P, and the half 
width w of the rocking curve at a height half way 
between the maximum and the level of the flat 
background. 


CALCITE 


The two similar pairs of crystals used are con- 
sidered equal to the best spectrometer calcites in 
this laboratory. One pair, 14—1B, has been pol- 
ished on the reflecting surfaces while the other 
pair U,—-T>2, remains as split. First order reflec- 


TABLE I. Reflection of Ag Ka by calcite. 











RX 10° 

w(seconds) Pp (radians) 
Polished surface (1A —1B) 20.5 30.0 9.63 
Split surface (U.—T») 3.3 30.6 2.25 
Perfect crystal (theoretical) 1.63 79 1.50 








tion data for the Ka radiation of silver were ob- 
tained and are listed in the first two rows of 
Table I. 

The theoretical values in the third row are 
those computed by S. K. Allison‘ on the basis of 





*S. K. Allison, Phys. Rev. 41, 1 (1932). 
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Darwin's theory of a perfect crystal as modified 
by Prins to take account of absorption. (The ef- 
fect of the modification is very slight at this 
wave-length.) In common with most of the calcite 
crystals which have been studied heretofore the 
present data deviate from the predictions for 
perfect crystals in a sense which indicates some 
degree of mosaic structure. Evidently these 
crystals are similar in quality to Allison’s pair 
designated II; they typify a high but not perfect 
grade of spectrometer calcite. 

Perfect crystal predictions are of course not 
applicable to crystals which have been polished, 
whatever their former state of perfection. Polish- 
ing in this case had the usual effect of increasing 
the integrated reflection of the crystal at the ex- 
pense of resolving power. 


ROCKSALT 


The rocksalt crystals used were from a fifty 
pound block of high quality, originating in Rus- 
sia. The surfaces used were chosen from a large 
number of freshly cloven pieces by visual inspec- 
tion and rough spectrometer examination. Not all 
of the pieces derived from this block were inter- 
nally clear and superficially plane but the two 
specimens selected for use had the uniform ap- 
pearance of fine calcite. Hence our rocksalt data, 
like those on calcite above, are not to be taken as 
representative of typical crystals of their kind 
but as pertaining to the best obtainable. Results 
for first-order reflection of the wave-length of 
silver Ka were as follows: 


R=18.2X10-5; 


w=87"': P = 17.7. 


The narrowness of the rocking curve shows the 
crystal to be very much more nearly perfect than 
any rocksalt crystal heretofore reported, though 
the width is still more than 50 times as great as 
the value for a perfect rocksalt crystal according 
to Darwin’s® treatment. The present value of 87” 
may be compared with a value of 900” obtained 
by Bragg, James and Bosanquet,* with about the 
same wave-length and with the 300” of Davis 
and Stempel.” 

Though much of the early x-ray spectrographic 
work was done with rocksalt crystals, it is not 


’ A. H. Compton, X-Rays and Electrons, (1926), p. 140. 
® Bragg, James and Bosanquet, Phil. Mag. 42, 1 (1921). 
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highly regarded at the present time for such 
purposes, since advancing refinement in measure- 
ment has increased the demand for higher resolv- 
ing powers. We think that the discovery of 
rocksalt crystals with resolving powers three and 
a half times as high as the previous best should 
do something to restore this crystal to favor in 
investigations where intensity rather than ex- 
treme resolution is the principal desideratum. 
While the resolution is by no means to be com- 
pared with that of good split calcite, the inte- 
grated reflection is eight times as great. 


ROCHELLE SALT 


The Rochelle salt surfaces used were adjacent 
portions of a large (001) face of an artificially- 
grown crystal from a laboratory supply house. 
The material is available in this way because of 
the interest in its piezoelectric properties. The 
face as received was a rectangle 17 cm long and 
9 cm in breadth; it was sawed into three pieces, 
all large for spectrometer use. 

Rochelle salt is a sodium-potassium tartrate 
with the formula KNaC,H,O;—4H.0. The struc- 
ture, apparently, has not been determined. The 
grating space normal to the (001) planes was de- 
termined by observing with a single crystal the 
separation of the silver Ka lines in the first and 
second orders and found to be 5.93A. Reflectivity 
determinations were made for both first and 


TABLE II. Reflection of Ag Ka by Rochelle salt. 











Order of 

reflection w (seconds) a RxX10° 
First 3.5 6.3 0.25 
Second 1.5 1 0.020 








second order reflections of silver Ka and the 
results are summarized in Table II. 

These data are notable for the extreme narrow- 
ness of the rocking curves, resembling those of 
the best calcite, together with the very low per- 
cent reflections and integrated reflections. Ac- 
cording to Darwin’s® theory of reflection by per- 
fect crystals the angular region of perfect 
reflection, a quantity closely related to the width 
of the rocking curve, is given by 


Ad) =n Fe? cot 6/xrmc’, (1) 
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where n is the number of scattering units in unit 
volume, F is the structure amplitude corre- 
sponding to the same scattering unit, and the 
other symbols are of obvious interpretation. The 
only uncertainty in evaluating A@ results from 
the uncertainty in F. Being in ignorance of the 
crystal structure, we can do no more than calcu- 
late an upper limit for F by adding the form 
factors’ of the constituent atoms of the chemical 
formula previously given. Such a calculation for 
the second order reflection leads to the value 101, 
which is the value that F would have for the 
second order reflection if temperature disturb- 
ances of the atomic positions did not occur and if 
the atoms were so located in the unit cell that 
they all scattered in phase. Of course neither of 
these conditions is fulfilled and consequently, on 
both counts, a lower value of F is to be expected, 
but we shall use the value 101 provisionally in 
Eq. (1), and obtain therefrom A@)=11.5X10- 
= 2.37’. Allison® shows that if the diffraction 
pattern of a single crystal is of the Darwin form 
the half width at half maximum of the parallel 
rocking curve should be w=0.60A@. The theoret- 
ical half width at half maximum under the as- 
sumption that F=101 is therefore 0.60 X 2.37 
= 1.4”, in striking coincidental agreement with 
the experimental value, 1.5”. This agreement 
may be taken as an indication that the approach 
of the crystal to perfection is comparable to the 
approach of F to the limiting value assumed. The 
low percent reflections do not, however, favor the 
presumption of a high state of perfection, but 
indicate that the narrow rocking curves result 
from a crystal structure permitting only a low 
degree of cooperation among the many atoms 
contributing to the scattering. It is perhaps un- 
profitable to speculate further upon this situation 
until some information about the crystal struc- 
ture can be obtained. In spite of the low reflecting 
powers observed these gigantic and inexpensive 
crystals with resolving powers comparable to 
those of calcite should not be entirely overlooked 
as possible spectrometer material for certain 
applications. 


7R. W. James and G. W. Brindley, Phil. Mag. 12, 81 
(1931). 
8S. K. Allison, Phys. Rev. 38, 203 (1931). 

















ABSOLUTE X-RAY 


BARITE 


Natural crystals of barite (BaSO,) from 
Felsoebanya, Hungary, with cleavage faces 3 by 
4 cm were mounted on the spectrometer in such 
a way as to utilize face reflections from the (001) 
planes. The length of the unit cell normal to 
these planes is given by Wyckoff® as 7.17A. If 
this value be taken as d in Bragg’s equation, it is 
found that no first order reflections occur, be- 
cause of the complete interference of identical 
planes of half the above separation. The first re- 
flection observed is therefore in the second order, 
and all barite reflections hereinafter discussed 
are properly specified by (002). 

The principal point of interest in connection 
with this crystal concerns the K limit of the 
barium constituent which comes within the 
available x-ray range. As it was desirable to select 
wave-lengths within this range with entire free- 
dom, the continuous spectra of tungsten and 
platinum were used as sources. The rocking 
curves obtained were not all as symmetrical as 
those obtained with good crystals such as calcite 
but bear more resemblance to the curves of rock- 
salt of average quality. On a few occasions 
singular excrescences appeared on the sides or 
feet of the curves, which could be traced to faulty 
areas on the crystal face and removed by install- 
ing horizontal slits or by moving the crystal. 
Table III presents the experimental results for 
the eleven wave-lengths investigated, together 
with data for one additional wave-length (0.71A) 
taken from the work of James and Wood.'® 


TABLE III. Barite reflection data and structure factors. 











(A) w (min.) P . R (rad.) X 108 F’ 
0.260 11.5 0.47 2.46 118 
0.295 7.0 0.50 1.54 99 
0.306 8.0 0.6 2.34 

0.317 10 0.5 2.44 

0.330 5.3 1.6 6.56 

0.350 s 1.8 10.1 109.8 
0.365 6.0 1.9 9.54 108.4 
0.434 6.0 1.5 8.99 114.7 
0.505 6 1.75 8.20 116.7 
0.574 6.5 1.7 7.32 118.2 
0.645 5.0 1.7 5.98 114.4 
0.71 6.63 125 








°R. W. G. Wyckoff, The Structure of Crystals (1931), p. 
289, 

10 R. W. James and W. A. Wood, Proc. Roy. Soc. A109, 
598 (1925). 
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Apparently there are no significant variations 
of the rocking curve width with wave-length, and 
the discussion to follow concerns the values of the 
integrated reflection. These are plotted in Fig. 1. 
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Fic. 1. Integrated (002) reflections of barite (BaSO,). 
a ry point, at the right is taken from James and 
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Fic. 2. Scattering factors of the barite unit cell (of four 
molecules) as derived from the integrated reflections of 
Fig. 1. (The three points of Fig. 1 nearest to the K limit do 
not permit the calculation of scattering factors because of 
uncertainty as to the appropriate absorption coefficients.) 


The abrupt change which occurs at the wave- 
length of the barium K limit is most conspicu- 
ously displayed, though its full magnitude is not 
shown on account of the rounding effect of the 
finite wave band represented by each point. This 
break is chiefly caused by the change in absorp- 
tion coefficient but a simple calculation shows 
that this change is insufficient to produce a dis- 
continuity of the magnitude shown in the figure. 
A discontinuity in the scattering factor is indi- 
cated. 

The last column of Table III and the curve of 
Fig. 2 exhibit the scattering factors of the unit 
cell, a group of four molecules, as computed by 
Darwin’s" equation expressing the integrated 


1 R. W. G. Wyckoff, The Structure of Crystals (1931), p. 
91. 
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reflection by imperfect crystals. This computa- 
tion suffers from the usual uncertainties regard- 
ing the effects of extinction and of thermal agita- 
tion of the atoms, both of which effects have been 
ignored. Consequently on both counts we must 
expect the experimentally derived values of F to 
be too low, and this expectation is borne out by 
an independent calculation of the scattering 
factor of the unit cell, upon the basis of the inde- 
pendently known scattering factors of the consti- 
tuent atoms and the structure of the barite 
crystal as determined by James and Wood." In 
such a calculation no account is taken of the 
special effect of absorption limits; so the calcu- 
lated result is to be compared with Fig. 2 only at 
the longest wave-lengths, where, apparently, the 
effect of the limit is no longer manifested. This 
calculation yields unit cell scattering factors 
from 150 to 157, depending upon which of several 
available lists of atomic scattering factors is 
adopted. These values are clearly in excess of 
anything which the curve of Fig. 2 could attain to. 

If we are to succeed in separating the scattering 
factors of barium from those of the rest of the 
crystal it is necessary to consider the cause of this 
discrepancy. It is not practicable to correct our 
scattering factors for the effects of thermal 
movements in this crystal but it is known that 
the correction would affect all ordinates of the 
curve of Fig. 2 proportionally and therefore leave 
the shape of the curve unchanged. Primary ex- 
tinction would have a like effect, but secondary 
extinction would alter the curve shape. Strong 
secondary extinction effects are judged to be im- 
probable, however, in a crystal with such strong 
ordinary absorption, since the extinction coeffi- 
cient affects the integrated reflection only as an 
additive term to an absorption coefficient al- 
ready large. James and Wood studied the similar 
crystalline sulphates of barium, strontium and 
lead and found discrepancies between theoretical 
and observed integrated reflections of the wave- 
length 0.710A of just the kind encountered here. 
Had this difference been due to secondary extinc- 


ROSS 


tion the discrepancy would presumably have been 
least in the strongly absorbent lead sulphate, 
whereas it was in fact greatest in this case. It 
seems correct to conclude that the shape of the 
curve representing the scattering factor of 
barium as a function of wave-length is as shown 
in Fig. 2, though the vertical scale of the plot is 
somewhat contracted and the constant scattering 
of the sulphur and oxygen atoms (about 55) has 
not been removed. 

Variations of atomic scattering factors with 
wave-length, in the neighborhood of the K limits 
with the ratio sin 6/X held constant, have previ- 
ously been measured by A. H. Armstrong" and 
by Wyckoft'® for nickel and copper, and by 
Glocker and Schafer,"* and Bradley and Hope" 
for iron. The present measurements, though 
performed on an element with twice as many 
scattering electrons, are in good agreement with 
the earlier work in disclosing a minimum in 
scattering power at or near the wave-length of 
the limit. Our results however do not show any 
maximum value of the scattering factor on the 
long wave-length side of the limit, such as 
Wyckoff'® described. In comparison with the 
predictions of the dispersion theory of Glocker 
and Shafer our curve (Fig. 2) has the predicted 
general shape on the long wave-length side of the 
limit, but, like the experimental curve of Glocker 
and Shafer, it rises too slowly on the side of 
shorter wave-lengths. 


@ A. H. Armstrong, Phys. Rev. 34, 931 (1929). 

3 R. W. G. Wyckoff, Phys. Rev. 36, 1116 (1930). 

4 R. Glocker and K. Schafer, Zeits. f. Physik 73, 289 
(1931). 

6 Bradley and Hope, Proc. Roy. Soc. A136, 272 (1932). 

% The evidence upon which this maximum is claimed is 
not considered strong. The curve of Fig. 3 (Phys. Rev. 36, 
1120 (1930)), showing the maximum, is not in agreement 
with the accompanying data of Table IV, sloping down at 
the right of the maximum three or four times as rapidly as 
the tabulated data would permit. This curve shows the 
values of the (200) reflection of nickel. Values of the 
(111) reflection given in Table IV and not plotted give no 
evidence of a maximum. 
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The cathodic disintegration of platinum by mercury ions 
has been measured for ion velocities up to 2000 volts. 
The mercury ions were isolated into a beam and allowed to 
strike a platinum surface under such conditions that the 
secondary electron emission could be measured as well as 
the disintegration produced. The disintegration per positive 


ion was found to increase linearly with the energy of the ion 
up to 2000 volts. Approximately two atoms of platinum 
are ejected per incident ion of 1000 volts energy. Curves 
are shown which give the characteristics of the secondary 
electron emission as a function of the energy of the ion. 





INTRODUCTION 


HE phenomenon of cathodic disintegration, 

or sputtering as it is commonly designated, 
has been utilized by many laboratories for 
deposition of thin films of metal. Despite its wide 
use, the facts concerning the process are not well 
established. It is generally agreed that the 
phenomenon results from the stopping of fast 
moving positive ions by the cathode. There is 
considerable disagreement among investigators 
even as to the ease with which different ions and 
different materials sputter. Most of this disa- 
greement arises from the manner in which the 
experiments are conducted. Nearly all experi- 
menters! have used the material to be sputtered 
as the cathode of a glow discharge. Our lack of 
complete understanding of this type of discharge 
has therefore prevented an understanding of 
cathodic sputtering. The stopping of a positive 
ion is in general accompanied by the emission of 
electrons from the material struck. These elec- 
trons will constitute a portion of the current 
measured if the material to be sputtered is in the 
field accelerating the positive ions. The total 
number of positive ions taking part in such a 
sputtering process is therefore in doubt. Cor- 
rection factors have been used? but, as the 
secondary electron emission is very much affected 
by surface conditions and history of the sample,* 


1Compton and Langmuir, Rev. Mod. Phys. 2, 186 
(1930). 

? Meyer and Guntherschulze, Zeits. f. Physik 71, 279 
(1931). 


’ Oliphant, Proc. Roy. Soc. A127, 373 (1930) 


these estimates may be very much in error for the 
particular material sputtered. The purpose of the 
present experiment is to isolate a beam of positive 
ions and cause it to enter a chamber in which the 
secondary electron emission may be measured 
upon the particular sample being sputtered, or 
may be eliminated from the quantities measured. 
This permits a determination of the amount of 
material sputtered per positive ion and its 
variation with the velocity of the ion. Such a 
method encounters a considerable sacrifice in the 
magnitude of the quantities measured. 


APPARATUS 


A diagram of the apparatus will be found in 
Fig. 1. The tube is in two sections. The adjoining 
ends of the sections fit into grooves on opposite 
sides of the aluminum plate A and are sealed in 
position with Dennison’s No. 2 hard sealing wax. 
A 14 mm opening through the aluminum plate 
joins the two sections. Across this opening is 
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Fic. 1. Diagram of apparatus 
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placed a grid of aluminum wire consisting of 
seven equally spaced parallel wires of 0.010 inch 
diameter. One section of the tube serves as the 
source of positive ions. The ions are obtained by 
ionizing mercury vapor between the filament F 
and the plate P which is of aluminum. The 
source of mercury vapor is the mercury in the 
side arm U. This mercury can be maintained at 
any desired temperature by a water-bath control. 
The mercury ions thus produced are accelerated 
through the grid across the opening in the 
aluminum plate into the section containing the 
plate C to be sputtered. The glass shield O serves 
to give some collimation to the positive ion 
beam. 

This beam of positive ions is incident upon a 
platinum plate C, the supports of this plate being 
entirely of aluminum. Aluminum is used here and 
elsewhere because of its known low rate of 
sputtering. The entire inside surface of this 
section of the tube is covered with an opaque 
conducting layer of platinum R obtained by 
cathodic sputtering previous to assemblage. 
Contact to this layer is made by the platinum 
spring 7. The platinum plate can be removed 
from the apparatus by means of the ground joint 
N. The arm B of the tube is immersed in liquid 
air to maintain the pressure in this portion of the 
tube as low as possible. Both sections of the 
apparatus are connected through liquid air traps 
to a vacuum system which includes a two-stage 
mercury diffusion pump. 


OPERATION 


The effective potential difference between the 
plate P and the filament F is between 75 and 80 
volts. The plate A is maintained at arbitrary 
values from 400 to 2000 volts negative with 
respect to F by means of a direct-current motor 
generator. The platinum plate C, which is to be 
sputtered, is maintained at a potential of 100 
volts positive with respect to A and R, to retain 
any electrons which might otherwise be emitted 
from C and thus be added to the reading of the 
positive ion current. Evidence that 100 volts is 
sufficient will be found in Fig. 2. This potential 
may be removed or reversed when it is desired to 
read the positive ion and secondary emission 
currents combined. 


AND EDITH GIDEON 


The water W, which controls the mercury 
temperature, is in general kept at about 57°C 
during the operation of the tube. In order to 
determine the order of magnitude of the effect 
due to the wax vapor, a run was made with the 
mercury U at room temperature. The pressure 
of the wax remained the same under these 
conditions, as the lower end of the plate A was 
still kept in cold water. The pressure of the 
mercury, however, decreased by a factor of 
fifteen. All currents were adjusted under these 
conditions to as nearly as possible the conditions 
which exist with the mercury at 57°C. The 
disintegration of the cathode then decreased to 
but 1/7 of the disintegration produced with the 
mercury pressure fifteen times as great, all other 
conditions remaining the same. We, therefore, 
conclude that the disintegration produced by 
positive ions of whatever constitutes the wax 
vapor is small compared to the effect produced by 
mercury ions. This is particularly true when it is 
considered that under the normal conditions of 
discharge the proportion of the positive ions of 
mercury to positive ions of wax vapor would be 
large. The effect of the wax vapor was to 
contribute to the magnitude of the positive ion 
current but not appreciably to the amount of 
material sputtered. Our values for the disinte- 
gration per positive ion may therefore be lower 
by a small amount than the true values. 

The magnitude of the positive ion current 
striking the plate is kept approximately the same 
for the different accelerating voltages. The 
platinum plate, being mounted through a ground 
joint, is removed before and after each run, for 
weighing. The loss in weight (determined by 
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Fic. 2. Showing the fraction of the secondary electrons 
retained by the plate C for various retarding potentials. 
Curve I, 600 volt Hg ions; curve II, 2000 volt Hg ions. 
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means of a Becker micro-balance) is taken as the 
amount of material sputtered. The length of 
time the platinum plate undergoes bombardment 
is noted, and this quantity multiplied by the 
current gives the amount of.charge and therefore 
the number of ions which have struck the plate. 
The ions are assumed to be singly charged.‘ 


DISCUSSION OF RESULTS 


Since the secondary electrons escaping from 
the plate C contribute to the current measured by 
the galvanometer Z, the secondary electron 
emission could be determined by measuring the 
combined positive ion current and secondary 
electron current obtained with no potential or 
small accelerating potential for electrons applied 
between plate C and the walls of the surrounding 
chamber, and subtracting from this the current 
to plate C when sufficient retarding potential is 
applied to hold the electrons on the plate. Such 
a method neglects the emission from the walls 
due to reflected ions and metastable atoms 
impinging there, or to other causes. Under the 
conditions of normal incidence prevailing in this 
experiment the number of reflected ions should be 
small. Experimental evidence, that electron 
emission from the walls due to any cause what- 
soever is negligible, was found in the fact that 
changing from no potential to an accelerating 
potential of as much as one hundred volts did not 
result in a measurable change in current to the 
plate C. In Fig. 2 are shown curves giving the 
fraction of the secondary electrons retained by 
plate C for different retarding potentials applied 
between the platinum plate C and the inner 
surface of the tube. It will be noted that ap- 
proximately 100 volts are necessary to hold the 
secondary emitted electrons on the platinum 
plate C for the highest velocity ions used. Two 
curves are shown, viz., for 600 volts and 2000 
volts. The intermediate voltage curves lie in the 
region between these two curves. The higher 
velocity ions result in a larger proportion of high 
velocity secondary electrons. Throughout the 
experiments on disintegration, 100 volts retard- 
ing potential was continuously applied so that 


‘ Bleakney, Phys. Rev. 35, 139 (1930). 
* Oliphant, Proc. Roy. Soc. A132, 631 (1931). 


the readings taken would represent a true 
positive ion current. 

Fig. 3 shows the number of secondary electrons 
emitted from the platinum per incident mercury 
ion over the range of 600 to 2000 volts. The 
values range from 0.3 to 0.5 for the platinum 
surface utilized in this experiment. The plate C 
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Fic. 3. The secondary electron emission from Pt per 
incident Hg ion for different energy ions. The platinum 
was removed from the apparatus between successive 
points. 


was removed from the apparatus for weighing 
between the observations plotted in this curve. 
When this same curve was taken without re- 
moving the plate from the apparatus, a smooth 
straight line was obtained, justifying the curve 
drawn in Fig. 3. The deviation of the observed 
points from this curve illustrates well the type of 
variations observed for secondary electron emis- 
sion, and shows the necessity of controlling the 
secondary electron emission in experiments of 
this character when extreme degassing of the 
parts is not possible. The characteristics of this 
curve are very similar to those observed by 
Oliphant* for helium ions incident upon cold 
molybdenum. The magnitude of the secondary 
electron emission observed shows the necessity 
for its consideration in determining the number 
of positive ions. The variation with energy of the 
ion makes essential the control of the secondary 
electron emission for experiments of the type 
under consideration. 
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Curve 1, Fig. 4, shows the disintegration 
expressed in atoms sputtered per impacting 
positive ion and in grams per coulomb as a 
function of the energy of the positive ions. A 
straight line is considered by us to be the most 
representative of the trend of the points. The 
deviation of the points from this line we believe 
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Fic. 4. The circles represent the number of atoms of 
platinum sputtered per impacting Hg ion. The crosses are 
the same data but neglecting secondary electron emission. 


to be due to errors which resulted mainly from 
handling and weighing of the plate. This may be 
realized when one considers that the greatest loss 
of weight measured was of the order of 0.005 
milligram. Disintegration by neutralized atoms 
which still possess velocity comparable to that of 
the impinging ions is a possibility which has been 
neglected in this experiment. Oliphant® gives 
evidence that disintegration by such atoms is far 
less intense than by positive ions with similar 
kinetic energy. Curve II is the one which would 
be obtained if the secondary electrons were 
permitted to contribute to the current measured. 
This curve is of the form obtained by other 


* Oliphant, Proc. Roy. Soc. Al24, 228 (1929). 
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observers when the secondary electron emission 
is not controlled.” ? 

Although the materials used prevent a direct 
comparison with the work of other observers, the 
similarities found by them for different materials 
and different ions is strong evidence that the 
same laws govern for all materials. The results of 
these experiments show that up to 2000 volts, the 
disintegration is a linear function of the energy 
of the impacting ion. This strongly suggests that 
the sputtering process is not a direct transfer of 
momentum, but is more likely a process involving 
the energy directly such as might be expected 
from evaporation resulting from a local high 
temperature condition. 

It is of interest to compute the efficiency of 
sputtering as a means for the evaporation of 
platinum. Taking the value of the disintegration 
at 1000 volts, the energy required to sputter one 
gram of platinum is 2.6X10° joules/gram. The 
energy required to heat and evaporate one gram 
of platinum at 4300°C is 3.3 10* joules/gram, 
giving as a means of evaporation of platinum an 
efficiency of a little more than one percent. The 
portion of the energy of the positive ion which is 
transformed into kinetic energy of the secondary 
electrons may be estimated from Fig. 2 and Fig. 
3. For 1000 volt positive ions this is approxi- 
mately one percent. This value for the secondary 
electron energy is undoubtedly higher than would 
have been obtained had it been possible thor- 
oughly to degas the platinum. The sputtering 
and secondary electrons can thus account for 
from two to three percent of the energy of the 
positive ions. Of the remaining energy of the 
positive ion incident upon the metal at normal 
incidence the major portion is undoubtedly given 
up as heat. Evidence on this point is given by 
Oliphant,‘ who reports that, with helium ions 
incident upon nickel at normal incidence, the 
energy given up as heat is that possessed by the 
ion to within five percent over the range from 40 
to 6000 volts. 


7 Von Hippel, Ann. d. Physik 81, 1043 (1926). 
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The Polarizabilities of Ions from Spectra 


JosePpH E. MAYER AND MARIA GOEPPERT MAYER, Chemical and Physical Laboratories, Johns Hopkins University 
(Received January 14, 1933) 


The polarizabilities of the gaseous alkali ions, and several 
other ions of rare gas structure, have been calculated by 
the Born-Heisenberg method from the corresponding 
spectra. Correction has been made for the effect due to 
“penetration” and for higher order distortion of the ion. 
The quantum mechanical derivation of the method has 
been given. The values of the polarizabilities of the alkali 


ions so obtained are found to be somewhat lower than 
those previously assumed. The polarizabilities of the 
gaseous negative ions must be correspondingly higher. The 
decrease in polarizability of the negative ion and increase of 
the positive ion on inclusion in a crystal is in agreement 
with known frequency shifts. 





INTRODUCTION 


HE sum of the polarizabilities of a negative 

and a positive ion in a crystal! or in solution 
may be obtained by measurements of the disper- 
sion. These sums indicate that the polarizability 
of a given ion is by no means independent of the 
state of the ion. It is different in solution and in 
crystal, and varies with the crystal. Fajans and 
Joos' have calculated values of the polarizabilities 
of the ions in the gaseous state, using reasonable 
but by no means certain assumptions. Pauling? 
by a semi-theoretical, semi-empirical method has 
calculated the polarizabilities of gaseous ions 
using the known polarizabilities of the corre- 
sponding rare gas and the assumed value of the 
negative ion. Since just the negative ions show 
polarizabilities extremely susceptible to the state, 
there is considerable doubt concerning their 
polarizabilities in the gaseous phase, and Paul- 
ing’s method cannot be considered unim- 
peachable. 

The only method of obtaining directly the 
polarizability in the gaseous state which may be 
called experimental, is that of Born and Heisen- 
berg,? which uses the values of spectral terms. 
The actual calculations made by them (1924) are, 
however, by means of the Bohr model of the 
atom, which introduces considerable error. Great 
improvements have also been made since in the 





1K. Fajans and G. Joos, Zeits. f. Physik 23, 1 (1924). 

*L. Pauling, Proc. Roy. Soc. All14, 181 (1927). 

*M. Born and W. Heisenberg, Zeits. f. Physik 23, 407 
(1924), 


values and assignment of, spectral terms. Al- 
though the method is frequently referred to in the 
literature‘ it appears that no comprehensive at- 
tempt has been made to recalculate polarizabili- 
ties by this method using the quantum mechan- 
ics. Further, no quantum mechanical derivation 
of the fundamental equation of this method, 
Eq. (1) of this paper, appears to have been pub- 
lished. Whitelaw and Van Vleck® have reported 
the conditions under which the proof can be 
made, but their complete work has not been 
published. It has been thought desirable to make 
the proof and calculate some polarizabilities by 
this method. Throughout, the term values given 
in the new book by Bacher and Goudsmit® have 
been used. 


METHOD 


The principle of the Born and Heisenberg 
method is simple. Consider the energy of a 
spectral term with high quantum numbers » and 
l, of an alkali metal. The state of the valence elec- 
tron may, to the zeroth approximation, be con- 
sidered to be given by a hydrogen-like function 
which is assumed not to have appreciable ampli- 
tude inside of the ionic kernel. That is, the elec- 
tron is assumed to be nonpenetrating. The poten- 


4 An excellent discussion is given by J. H. Van Vleck, 
The Theory of Electric and Magnetic Susceptibilities, page 
215. 

5 N. G. Whitelaw and J. H. Van Vleck, Phys. Rev. 41, 
389, Abstract of Report to the New Haven meeting (1932). 

6 R. F. Bacher and S. Goudsmit, Atomic Energy States, 
McGraw-Hill, New York (1932). 
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tial on the electron is that due to the singly 
charged undeformed ionic kernel, plus a dipole 
perturbation potential due to the deformed ion, 
polarized under the influence of the valence 
electron. The zeroth order electrostatic potential 
of the undeformed ion would account for exactly 
the energy of a hydrogen state with the same 
quantum numbers. The perturbation method 
allows us to calculate the decrease of energy due 
to the polarization of the ion, 


—AE=}ae(r-*), (1) 


In this equation r~ is the average value of r~‘ for 
the hydrogen state, that is, 


froin. u(r) | 2dr 
ri= (2) 
frirs, i(r) | °dr 





if R,, x(r) is the (unnormalized) radial hydrogen 
function with the quantum numbers 2, /, of the 
term in question. By use of experimental devia- 
tions of term energies from those of hydrogen 
terms the polarizability a of the ion may be cal- 
culated. This is, of course, the polarizability of 
the ion under the influence of a field with the 
frequency of the atomic term in question. In all 
cases treated in this paper the frequency is so far 
below the resonance line of the ion that this cor- 
rection is less than one percent. 

The justification and details of the calculation 
will be given under the heading “‘Calculations.”’ 
By this method the calculated polarizabilities for 
any single ion are found to vary considerably 
with the term used. In general the calculated 
polarizability is found to increase with the total 
quantum number » in a series of constant /, and 
to have lower values if a series of higher / is used 
in the calculation. 

However, the values of a so calculated from 
the P and D terms of Li I, Be II, B II] and CIV; 
the D and F terms of Na I, Mg II and Al III; 
and the F and G terms of K I, Rb I, Cs I, Ca IT, 
Sr II, are sufficiently close to the probable values 
of a for the corresponding ions, and show a suffi- 
cient regularity in the trend to give hope that by 
considering second order effects reliable values 
may be obtained. Two corrections have then 
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been made, a correction for the perturbation 
energy due to higher order distortion of the ion, 
and a correction due to penetration. 

The effect on the energy due to higher than 
dipole distortion of the ion can be shown to be 
due to a constant of the ion which may be termed 
its ‘“quadrupolarizability,” the ratio of which to 
the “‘dipolarizability” can be estimated. 

A further perturbation which is present even 
were the ion undistorted may be termed the 
penetration perturbation. If Z(r), the effective 
nuclear charge at the radius r from the center of 
the ion, be defined as the net charge, in units of 
the charge on an electron e, lying within the 
radius 7, then a perturbation potential at r exists 
Unpen(r) defined by the equation 


dU e*(Z(r)—Z.,) 


dr r? 





Although Z(r) has been calculated for many ions 
to obtain form factors which determine x-ray 
reflection intensities, values for only a few ions 
have been published. The values for Nat, Cl-, 
He and Rb* are found in an article by Hartree,’ 
for Lit, in an article by Hargreaves,* and a curve 
from which the values for K*+ may be obtained in 
an article by James and Brindley.’ Some other 
curves of electron density obtained from ex- 
perimental data have been published"® but are of 
little value for the range of large radii which are 
here important. The authors are indebted to 
Professor J. C. Slater for a copy of a table of Z(r) 
values for Cs+, obtained by him from Hartree. 

Graphical integration of the perturbation po- 
tential times ¥”, which gives the perturbation 
energy, indicates that the perturbation at large 
radii alone is important (r>1.5 in units of 
dy = 0.5285 X 10-8 cm). For this portion of the ion 
an empirical analytical expression for Z(r) may 
be obtained, the constants of which can be extra- 
polated for other ions. 

The calculations are given in more detail under 
that heading. In Table I the results of the calcu- 


7D. R. Hartree, Proc. Camb. Phil. Soc. 24, 89 (1927). 

8 J. Hargreaves, Proc. Camb. Phil. Soc. 25, 75 (1928). 

9R. W. James and G. W. Brindley, Proc. Roy. Soc. 
A121, 155 (1928). 

10 R. J. Havighurst, Phys. Rev, 29. 1 (1927). 
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TABLE I. 
Element, T/Z?—R/n* 
spectrum, TermvalueT = =AToxp/Z*  ATpen/Z AT, /Z* ao X 10% aX 10 
and charge Term cm" cm"! cm! em7! cm Cy cm 
Li 2P 28 582.5 1150.3 187.0 963.3 0.0310 ? 
Sp. Li I 3P 12 560.4 368.4 73.2 295.2 0.0289 ? 
Z=1 4P 7 018.2 160.1 32.2 127.9 0.0284 ? 
3D 12 203.1 11.1 0.5 10.6 0.0258 1.05 0.0246 
4D 6 863.5 5.4 0.2 5.2 0.0273 1.06 0.0258 
5D 4 389.6 0.5 0.0 0.5 0.005 1.07 
Bett 2P 114 959.6 1302.4 279.0 1023.4 0.00825 ? 
Sp. Be II 3P 50 384.4 402.4 73.0 329.4 0.00805 ? 
Z=2 4P 28 122 172 33 139 0.00785 ? 
3D 48 828.5 14.8 1.0 13.8 0.00841 1.11 0.00756 
4D 27 459.5 6.7 0.9 5.8 0.00746 1.14 0.00656 
5D 17 570.4 3.4 0.5 2.9 (0.00696 1.15 0.00605 
+++ 2P 257 555.5 1170 190 980 0.00340 ? 
Sp. B III 3P 112 976.8 356 62 294 0.00320 ? 
Z=3 4P 63 100.1 153 27 126 0.00312 ? 
3D 109 861.0 14.4 1.1 13.3 0.00361 1.16 0.00311 
4D 61 794.7 7.75 0.67 0.67 0.00405 1.20 0.00338 
5D 39 541.8 4.27 0.32 3.95 0.00421 1.21 0.00348 
4F 61 731.9 87 0.00 0.87 0.0035 1.01 0.0035 
Ctttt 2P 435 637.9 919 153 766 0.00154 ? 
Sp. CIV SP 200 110.4 314 48 266 0.00162 ? 
Z=4 4P 111 858.7 133 22 111 0.00155 ? 
3D 195 292 13.2 1.1 12.1 0.00185 1.20 0.00154 
4D 109 837.6 6.6 0.6 6.0 0.00193 1.25 0.00154 
5D 70 286.5 3.6 0.3 3.3 0.00198 1.26 0.00157 
4F 109 743.7 0.7 0.0 0.7 0.00157 1.02 0.00154 
Na* 3D 12 276.18 83.45 8.13 75.32 0.184 1.11 0.165 
Sp. Na I 4D 6 900.35 41.94 4.66 37.28 0.192 1.13 0.170 
Z=1 5D 4 412.47 23.09 2.82 20.27 0.194 1.14 0.170 
4F 6 860.37 1.96 0.02 1.95 0.070 1.01 0.069 
5F 4 390.37 0.99 0.01 0.98 0.062 1.01 0.061 
Mg*t 3D 49 777.0 251.5 31.5 220.0 0.134 1.27 0.105 
Sp. Mg II 4D 27 955.3 130.4 18.0 112.4 0.144 1.34 0.107 
Z=2 5D 17 846.3 72.2 11.2 61.0 0.146 1.36 0.107 
4F 27 467.4 8.43 0.17 8.26 0.074 1.03 0.072 
SF 17 577.2 4.92 0.14 4.78 0.076 1.03 0.074 
OF 12 204.8 3.02 0.12 2.90 0.074 1.03 0.072 
Al*+++ 3D 113 497.9 418 46.5 371.5 0.100 1.47 0.069 
Sp. Al III 4F 61 841.7 12.89 0.39 12.50 0.050 1.05 0.047 
Z=3 5F 39 578.6 8.24 0.31 7.93 0.056 1.06 0.053 
OF 27 484.5 5.64 0.23 5.41 0.061 1.06 0.058 
5G 39 526.2 2.41 0.01 2.40 0.070 1.00 0.070 
K* 4F 6 878.5 20.0 0.14 19.9 0.715 1.03 0.694 
Sp. KI SF 4 404.2 14.8 0.1 14.7 0.925 1.04 0.890 
Z=1 OF 3 056.5 8.3 0.1 8.2 0.840 1.04 0.807 
oo 4F 27 694.0 65.0 1.5 63.5 0.571 1.06 0.540 
Sp. Ca II SF 17 714.1 39.1 ia 37.9 0.600 1.08 0.555 
Z=2 OF 12 290.0 24.3 0.9 23.4 0.600 1.09 0.550 
6G 12 211.0 4.54 0.0 4.54 0.46 1.01 0.46 
Rbt* 4F 6 987.6 39.0 1.0 38.0 1.37 1.05 1.30 
Sp. Rb I 5F 4418.2 28.7 0.8 27.9 1.77 1.06 1.66 
Z2=1 6F 3 068.0 19.7 0.6 19.1 1.95 1.07 1.83 
Sr++ 4F 27 960.4 131.5 8.6 122.9 1.10 1.15 0.96 
Sp. Sr II 5F 17 896 84 7.0 77 1.22 1.19 1.02 
Z=2 6F 12 412 55 5 50 1.28 1.21 1.06 
Cst 4F 6 932.80 74.27 2.73 71.54 2.58 1.07 2.38 
Sp. Cs I 5F 4 433.00 43.55 2.26 41.29 2.61 1.09 2.39 
Z=1 OF 3 074.77 26.54 1.58 24.96 2.56 1.10 2.33 
5G 4 393.5 4.05 0.01 4.06 1.1 1.02 1.1 
6G 3 057.0 8.77 0.01 8.76 3.6 1.02 3.5 
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lation are given. The method is to subtract the 
experimental term value 7 divided by Z? from 
the hydrogenic value, R/m*. This difference is 
given in column 4. The part of this, AT y.n/Z?, cal- 
culated as due to penetration, tabulated in 
column 5, is subtracted, leaving A7,/Z* (column 
6). This AT,/Z? is the term defect, divided by the 
ionic charge squared, due to distortion of the ion. 
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From this a, the polarizability, is calculated, 
(column 7) without consideration of the effect 
due to the quadropolarizability. The effect of the 
higher distortion is then taken into account by 
dividing by a factor C,, the expression for which 
is given in Eq. (14) and which is tabulated in 
column 8. In column 9 the final value of the 
polarizability a is given. 





CALCULATIONS 


(a). The deformation effect 


The complete zeroth order function W, of the atom, neglecting interaction between the ionic 
kernel and the valence electron is written as a product of the ionic function y,,' with the quantum 
numbers m, and a hydrogen-like function y,” for the valence electron. p is an abbreviation for the 
quantum numbers », /, and the orientation quantum number m;. We wish to calculate the perturba- 
tion of a state k=0(m=0, p= py). 


' .. 3.,.,2 
Ww, - Vm'Vp"s Wy - Yo'p,"s Ey = En ’ + E,’. (3) 


The perturbation potential /7 due to the distortion of the kernel may be written 


R-r 1 1 R:-r\” 
Hie= —e —, us=——e—|3(—) -R'| (4) 


r® 2 #F r 


H=H*+H? with 


R = > Ti, 


in which r and R designate the vector radii of the valence electron and the ion, respectively. Cor- 
rectly R is a sum over the vector radii r; of all of the electrons 7 in the kernel. 7 and R are the re- 
spective magnitudes. //* is the dipole, //° the quadrupole distortion. 

If we define Ho, by the equations, 


R=|R| 


Hox = Hox? +1 o?, Hy = f W//¢w dr, Iq? = f Wl wd, (5) 


all products Ho," X I/;0° as well as the first order energy JJ) are zero because of the form of //. Further, 
the degeneracy in the orientation quantum number m, of the hydrogen state is not lifted by the 
perturbation, and the states of different / are the correct zeroth order approximations to the problem. 
Applying the conventional method of the perturbation theory, one obtains 


0 0 0 
AE» = AE, -_ p ag (11*p,, nel 24 | IT»), nl ?) /(Eo— En'— E,”) ]. (6) 


0 0 0 
The energy term £,” which is small compared to Ey)—E,,' is neglected. (The term m=0 does not 
occur in this sum.) The sums over # are then evaluated by observing that 


r 2 1/71 R-r\? 
| Ron (-) +- | ~ E (—) -Ren| 

r3 ee 4llr3 r ” aig 
2 2 e! 


{(R ) 1 [3(~*) P "| IR 3 
=-—é lame +— “ —— hae * al re rites 7 he yt) Dy a ae f° n®| *(97-* VoPo T° 
0 3 4p , % 0 = 3 ( ( Pol 5 om” | (1) | 


The complete expressions for the change of energy due to deformation becomes 


DL | L1%0p,, mp | 4 | 1° op,, mp! t= —e) » 
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0 0 0 0 
—AE - (e*/3) Cue) OD OP I | Rom | ?/(Eo— En) + : (7),.5.20. | Rom? | 2/(Eo- En) ] } ° (8) 
The polarizability a of the ion is given by 
0 0 
q= rp we | Rom | */(Eo- En) }. (9) 


It is now desirable to obtain an approximate value for the second sum in Eq. (8) in terms of the 
polarizability. A customary approximation is to assume that all the significant terms in Eq. (9) occur 
at approximately the same energy, obtaining, 


LmL| Rom|2/(Eo—Em)] 2Roo?/Ex (10) 


in which E, is some average energy, approximately the ionization energy. A similar approximation 
leads to the equation, 

—— 0 Root (Roo?)? Roo? 

> ml | Rom |?/(Eo— on) |S os = = = 

E, NEw N 








0 0 
py | Rom | 2/(Eo _ En) }, 





(11) 


in which N is the number of outer electrons. The approximation (11) has been obtained in the follow- 
ing manner, both Roo? and Ry are the sums of the average of these functions of r for the individual 
electrons 7 of the ion in its normal state. That is, 


Roo? = >. i(7 00"), Root = ds(r 00°) x. (12) 


There are N outer electrons of the ion whose individual 7? and fo* will be about equal, and these 
functions of all other (inner) electrons will be negligible. root will be about equal to, and somewhat 
larger than, the square of 7°, making Root=(Ro*)?/N. Further the average energy E, will be some- 
what larger than E,. Numerically the values N =8 (and 2 for Li+) have been chosen in making the 
arithmetical calculations. 
Substituting (9) and (11) in (8) one obtains for the change in energy, 
2 (p-6 
—AE =; race b Po Ree nore] 
5 N ta) 





(13) 


0 
The unperturbed energy is that of a hydrogen like state, —E=Rhc/n,?, and the term value T is 
the negative of the energy divided by fc. One may then write the equation for a, 





a= 


2(R/no? —T,,)he 
— -=B,Z,-'*C"LR/no?—T>, |, 


eee a (14) 
e?(r~4) ppg 1 +3Ro0?(7-*) p,n, SN (7~*),.»,] 


1C 
B, Ze ee (r ) Pers .. 
e 


Cc. - [1 +3Ro?(r ®) .o./SN(9-) o.0, }- 





The values of (r~“*),,», and of (r~*),,», may be 
easily evaluated knowing the hydrogen radial 
functions. Waller" has given an expression for 
(7~*) pn: Ro® is connected with the specific mag- 
netic susceptibility x by the equation 

oo” = (6mgc?/ Ne*) x. (15) 


The values of x used are the experimental values 


1], Waller, Zeits. f. Physik 38, 635 (1926). 


tabulated by Pauling,” except for the first row 
of the periodic system, for which Pauling’s 
theoretical values are used. The bracket in the 
denominator of Eq. (15), which is abbreviated 
as C,, the correction for higher distortions is 
tabulated in Table I, column 8. (7~*),,», is di- 
vergent for P states and the correction cannot 
be made. This merely means that the P orbits 


2 L. Pauling, Phys. Rev. 34, 954 (1929), 
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penetrate the ion too strongly for this method of 
approximation. 


(b). The penetration effect 

Instead of using the difference between the 
experimental term value 7 and the hydrogen-like 
term value in the numerator of Eq. (15) the 
change in 7, AT,.n, due to penetration has been 
subtracted. This effect is due to the deviation in 
the true potential from the assumed zeroth ap- 
proximation due to the undistorted ionic kernel 
and may be handled independently of the dis- 
tortion effects. 

If Z(r) is the effective nuclear charge as a func- 
tion of r and Z,, the charge on the ion, the per- 
turbation potential may be written. 


dU/dr=e?(Z(r) —Z,,)/r’, 


Z(r) —Z., © 1 dZ(r) 
= y-e|—— -{ ~ ar| 


r  &* ae 


(16) 








The first order effect on the term value alone is 
considered. 
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. AE e 7 
AT = ——=—— | ¥,°Uy,,"dr, (17) 
hc he 


Z(r) —Z,, is assumed to have the analytical form 
Bbr’e-8", in which case U(r) will have the form 
U(r) = — (e*b/ao)e-*", when r is measured in units 
of a, the radius of the first Bohr orbit of hy- 
drogen, (a) =0.5285 X10-§ cm.). This empirical 
equation for Z(r) can be made to fit all the given 
tables of Z(r) for large values of r by the proper 
choice of the constants 0} and 8. Since the whole 
effect is at large r values it is unimportant that 
the equation is absolutely false for small r. No 
tables were at hand for the multiply charged 
ions. It was assumed that the equation 8 = 2Z,/n 
was approximately true with ” the total quantum 
number of the outer electrons of the ion, and Z, 
the effective charge on the electrons. Z, was then 
assumed to increase by one as the charge on the 
ion increased by one, for the same configuration. 
In Table II are given the assumed values of } 
and £6 for the various ions. A few graphical inte- 
grations showed that the total penetration effect 





TABLE II. Empirical constants of the penetration potential. 




















Ion Li* Be** BY Crt Nat Mgtt Al*+t+ K* ce. Rbt or? Cst 
b 9.5 9.5 9.5 9.5 33 33 33 34 34 35 35 36 

B 5.43 7.43 9.43 11.43 4.05 5.05 6.05 3.09 3.76 2.42 2.92 2.11 
occurred in the region where the assumed form TABLE III. Polarizabilities of the gaseous ions. (a X10" 


of the function Z(r) was in agreement with the 
tabulated values. 


RESULTS 


In Table III are given the most probable values 
of a, as deduced only from the internal evidence 
of this calculation. Values previously published 
in the literature are included for comparison. The 
values for the rare gases are experimental, and 
quite accurate. They are copied from Pauling’s 
tables. 

Some considerations lead to misgivings about 
the accuracy of the values here computed. The 
consistently lower values of a obtained from a 
higher series in comparison with a lower series, 
F compared with D in sodium for instance, indi- 
cates some sort of systematic error. This may be 
a gross underestimation of the quadrupolariza- 


in cm*.) 








Born and 

















Fajans_ Heisenberg 

Ion This paper Pauling and Joos (k=/+}) 
He 0.204 
Lit 0.025 0.0295 0.079 0.075 
Bett 0.007 0.0080 0.039 
ptt 0.0033 0.00303 
ata a 0.0015 0.00135 
Ne 0.396 
Na* 0.17 0.182 0.197 0.21 
Mg*+ 0.10 0.092 0.110 0.12 
ae 0.53 0.054 0.065 
Ar 1.645 
Kt 0.80 0.844 0.879 0.87 
Catt 0.54 0.474 0.447 
Kr 2.49 
Rb* 1.5 (?) 1.42 1.485 1.81 
Se** 1.0 0.87 0.883 1.42 
Xe 4.05 
ie 2.35 2.45 2.577 2.79 
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bility correction, which is not unlikely since our 
approximation gives a lower limit to this correc- 
tion, or a gross underestimation of the penetra- 
tion effect which is less likely. Any such explana- 
tion would tend to give more weight to the 
values obtained from the higher series, indeed the 
true value would then be slightly less than that 
of the higher series, which seems unreasonable in 
the case of sodium at least. 

For some ions an error may be introduced by a 
false estimation of the series limit. That is, al- 
though the differences in term values are well 
known, the absolute values are uncertain. That 
such uncertainties exist is apparent. Bacher and 
Goudsmit and the International Critical Tables 
differ by 2.6 cm™ in the case of K I, and 98 cm 
in Sr II. However the term values of Na I are 
probably known with great accuracy. 
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The results of this paper agree quite well with 
the calculations of Pauling. It is particularly to be 
noted that the agreement for the ions of the 
elements of the first row of the Periodic Table 
lends great support to Pauling’s method, which 
for these ions is purely theoretical and assumes no 
empirical constants. 

For most of the singly charged positive ions, 
(alkalis), the polarizabilities appear to be slightly 
lower than Pauling’s values, and considerably 
lower than those of Fajans and Joos. Pauling has 
obtained the values which he gives by assuming a 
linear extrapolation of a screening constant, cal- 
culated from the polarizabilities of the corre- 
sponding rare gas atom and negative halide ion. 
If the deviation above noted is real, it would lead 
to the conclusion that the polarizabilities of the 
gaseous negative ions are greater, and indeed con- 





TABLE IV. 
Term 2P 3P 4P 3D 4D 5D 4F 5F 6F 5G 6G 6H 
Buy X 104 0.322 0.978 2.228 24.4 51.5 96.0 360 634 1025 2640 4080 12850 





siderably greater, than those assumed by Pauling. 
This means that the gaseous negative ions have 
considerably higher polarizabilities than the same 
ions in crystals, and the gaseous positive ions 
somewhat lower a’s than in the crystals. 

If the approximation of Eq. (10) is used, 
leading to 


a 2$e?(Roo?/Ea), (18) 
it might be assumed that Ry’? is approximately 
independent of the state of any given ion and that 
the polarizability @ varies inversely with the 
“main frequency.’ The absorption frequencies of 
negative ions in a crystal are more than double 


the probable absorption frequencies (near the 
electron affinity), in the gaseous state. The ab- 
sorption frequencies of the positive ions, on the 
other hand, are probably lowered by inclusion in 
a crystal. It seems then, qualitatively, probable 
that the polarizabilities of positive ions are in- 
creased somewhat in passing from the gaseous 
into the crystalline (or aqueous) state, while those 
of the negative ions are enormously (twofold?) 
decreased by the same change of state. 

For convenience the values of B,,,, of Eq. (14) 
are included in Table IV for different terms. 

In conclusion the authors with to thank Profes- 
sors M. Born and K. F. Herzfeld for helpful sug- 
gestions and criticism. 
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Effects of Chemical Combination on the X-Ray Emission Spectrum of Sulphur 


JosEPH VALASEK, University of Minnesota 
(Received February 28, 1933) 


With an improved method to eliminate the effects of 
chemical reactions with the target material, the wave- 
lengths of the sulphur Ka and @-lines were measured for 
the following substances: S, FeS, CoS, NiS, CueS and ZnS. 
The superiority of the new technique over that usually 
employed is particularly evident in the spectrum of ZnS in 


which the K@-line is observed as a triplet instead of the 
single broad line heretofore observed. This result is inde- 
pendent of the target material. The K@-lines of CoS and 
NiS are observed for the first time. Two crystal modi- 
fications of ZnS, namely, zincblende and wurtzite, gave 
spectra which showed no appreciable differences. 





HE first observations showing that the x-ray 

emission spectrum of an element depends 
on chemical combination were made in 1924 by 
Lindh and Lundquist. The effect is very small ex- 
cept in the case of the lighter elements, such as 
silicon, phosphorus, sulphur and chlorine. A 
summary of the work up to 1928 may be found 
in Lindh’s book on Roentgenspektroskopie. Re- 
cent studies of the subject are those due to 
Lundquist,' Faessler,2, Deodhar* and Svensson.* 
Deodhar, in particular, finds several new lines 
which are emitted by some compounds and not 
by others, indicating the creation of new energy 
levels in an atom when it is combined with one or 
more other atoms. In addition to this appearance 
of more or less characteristic lines, there are often 
shifts of a few X-units in the wave-lengths of 
some of the lines. 

Although sulphur compounds have been in- 
vestigated several times, in all but one of these 
studies,? and that one restricted to the intense 
a, and a»-lines, the spectrum was excited by 
direct electron impact on the powdered com- 
pound. This was rubbed into scratches on the 
face of the target. What happens to the com- 
pound under these conditions is only now begin- 
ning to be clearly understood.* The chemical re- 
action between the substance and the target is 
evident from the fact that the ‘‘same”’ substance 
often gives a different spectrum depending on 
the underlying target metal. 


1 Lundquist, Zeits. f. Physik 60, 642 (1930). 
? Faessler, Zeits. f. Physik 72, 734 (1931). 

3’ Deodhar, Proc. Roy. Soc. A131, 647 (1931). 
‘Svensson, Zeits. f. Physik 75, 120 (1932). 


In order to overcome such effects it has been 
often suggested that the spectrum of the com- 
pound should be excited by fluorescence instead 
of electron impact. This is the procedure followed 
by Faessler, who used one of the new type of 
x-ray tubes developed by Coster and Hevesy. In 
this form of tube, the powdered material is 
rubbed on a flattened portion of the electron 
focussing shield which surrounds the filament. 
However, Faessler’s work covers only the intense 
a, and a-lines of sulphur, presumably because 
the other lines could not be obtained with suffi- 
cient intensity. 

In order to study the less intense lines, par- 
ticularly the 8-lines which show the greatest de- 
pendence on chemical combination, the writer 
has used the target illustrated in Fig. 1. The sub- 


Ni Co 





Fic. 1. X-ray target. 


stance S is protected from the direct impact of 
the electron stream by the copper shield C. The 
focal spot is one to two millimeters from the front 
of the material. This does not give a pure 
fluorescence spectrum, for many of the electrons 
no doubt are scattered from the target to the 
substance. However, the opportunity for chem- 
ical reaction is greatly reduced as the results 
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show. The spectrum is of fair intensity. The usual 
exposure for the B-lines was 8 to 10 hours at 10 
to 15 milliamperes and 12 to 15 kilovolts. In the 
case of pure sulphur, the current was reduced to 
1 to 2 milliamperes to keep the sulphur from 
melting and running out onto the target. 

A vacuum spectrograph of the Siegbahn type 
was employed. The crystal was of calcite, the 
wave-lengths being measured with the aid of 
comparison spectra. For the sulphur Ka-lines, 
cobalt Ka, was used in the third order, while for 
the KB-lines nickel Ka; was employed. These 
lines have been measured with ample accuracy by 
Eriksson.® The Bragg angles for these lines, with 
Larsson’s value for the lattice constant of calcite 
in the third order,® are, at 18°C: 

Ni Ka; 55° 0’ 25” Co Ka; 62° 7’ 35”’. 
The angular separations between the reference 
lines and the sulphur lines were calculated from 
the distances along the plate as measured by a 
Gaertner comparator and the radius of the cam- 
era measured in the manner outlined by Lindh.’ 
Temperature differences between the sulphur 
and comparison exposures were corrected for by 
the equation A@ = 2.10’’ AT tan 6. From the glanc- 
ing angles so found, the wave-lengths of the sul- 
phur lines were calculated with the aid of Lars- 
son’s value of the lattice constant for the first 
order.® 

The wave-lengths of the various lines in the 
spectrum of sulphur and the various sulphides 
studied are given in Table I. 


TABLE I. Wave-lengths of the x-ray emission lines of sulphur 

















in X-untts. 

Line Ss FeS CoS NiS CueS ZnS 
a. 5363.50 5364.35 5364.16 5364.05 5363.97 5364.27 
a: 5360.39 5361.28 5361.34 5361.05 5360.94 5361.29 
e’ = 342.1 5343.0 5342.7 5342.1 5342.4 5342.7 
a3 —«-455330.21. 5330.75 = $330.77. = $330.49 $330.39 = 5330.61 
a3’ 5327.5 5328.2 5328.4 5328.1 5328.2 5328.1 
a 5323.27 $324.12 $324.04 5323.52 5323.45 5323.89 
as 5284.5 5285.2 5286.0 5284.2 5285.1 5285.4 
Bs — 5065.24 5065.26 5066.02 _ — 
Bs — —_ — —_ = 5030.4 
6: 5020.72 — $020.89 — 5021.84 5023.77 
Biz — 5017.3 — 5018.4 — 5016.95 
Bz 5013.35 _ 5012.62 — 5013.39 —_ 

5 Eriksson, Zeits. f. Physik 48, 360 (1928). 

®Siegbahn, Spektroskopie der Roentgenstrahlen, 2nd 


edition, p. 44. 
7 Lindh, Roentgenspektroskopie, p. 63. 
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The a, and ae-lines are probably correct within 
+0.1X, while the others may be in error by +0.1 
up to +0.4X according to their breadth and in- 
tensity. The a;-lines were particularly difficult to 
measure. 

Other faint lines which were occasionally ob- 
served and seem to be due to sulphur or sulphides 
are as follows: 5108, 5043, 4953, 4936, and 4900. 
These values may be in error by +2X. Of these, 
the 5108 and 4900 lines were most frequently ob- 
served. 

The a-group is seen to change very little from 
compound to compound. This is surprising in 
view of the fact that the K absorption edges shift 
by 2 to 4X-units.* Evidently these shifts are 
closely paralleled by the shifts in the L levels. 

On the other hand, the 8-lines show widely dif- 
ferent appearances in some cases. This is clearly 
illustrated in Fig. 2 which summarizes the re- 
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Fic. 2. The K§-lines of sulphur and sulphides. 
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sults for these lines. The spectra of CoS and NiS 
have not been previously observed. The CoS was 
prepared by precipitation from cobalt nitrate by 
ammonium sulphide. The NiS was the mineral, 
millerite. The spectrum of ZnS is entirely differ- 
ent from the single broad line heretofore de- 
scribed in the literature. The spectrum here ob- 
tained was the same when the copper target was 
replaced by iron. The spectrum of pure sulphur 
showed only a slight dependence on the target 
material. This result could not have been ob- 
tained by the usual procedure of electron impact 
directly on the powdered sulphur. 

However, extreme care must be exercised if the 
electrode influence is to be entirely eliminated, 
particularly if the substance is easily fused or 


8 Lindh, Roentgenspekiroskopie, p. 291. 
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volatilized. This is due to the almost unavoidable 
contamination of the target which usually gives a 
scarcely perceptible target sulphide spectrum 
superposed on that of the compound. This can be 
reduced to negligible proportions by applying 
proper care. The target must be kept smooth 
and clean. 

It is surprising that the frequently observed 
B;-line could only be seen on two of the plates, 
namely, two heavily exposed ZnS plates, and 
that on these it was scarcely visible. On the 
other hand, two of Deodhar’s recently discovered 
“‘weak”’ lines, namely §,; and 8;, were quite strong 
in some cases. The {,-line is a prominent line in 
the CoS and FeS spectra, while 8; appeared in 
the ZnS spectrum. 

The reasons for the observed differences in the 
B-lines of various sulphides are not clear. One 
may think of a readjustment of the outer energy 
levels of the combining elements to some com- 
mon mean value. It would be of interest to study 
this question by means of the spectra of the 
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metallic elements in the sulphides. A start has 
already been made on this problem, and it has 
been found that the K§,-line of iron is strongly 
shifted in the sulphide to 1758.46X, while its 
wave-length in the spectrum of pure iron js 
1753.01X. In energy units, this shift is approx- 
imately equal to the difference between the S Kg, 
and S KBix-lines and is in the same direction. A 
further study of these effects is in progress. 

The method of examination here described 
makes it a relatively simple matter to compare 
the spectra of various crystal modifications of the 
same compound. Insofar as this has been done, 
(ZnS as zincblende and wurtzite) there seems to 
be no appreciable effect of crystal modification 
on the x-ray spectrum. This study is being con- 
tinued. 

This investigation was aided by a grant from 
the Research Fund of the Graduate School of the 
University of Minnesota. The writer wishes to 
acknowledge the assistance of Mr. John W. 
Liska with a part of this work. 
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Projective Relativity and the Einstein-Mayer Unified Field Theory 
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The relationship between the Einstein-Mayer unified 
field theory and the projective relativity theory is de- 
veloped and it is shown that the field equations (47) of the 
former theory may be interpreted as conditions that the 


connection employed be reducible to a symmetric form by 
means of a suitable non-holonomic factor change. Grounds 
are found for questioning an assumption made by Schouten 
and Van Danzig. 





§1. INTRODUCTION 


N 1921 Kaluza put forward the idea that the 

addition of a specialized fifth dimension to 
the space-time of the general theory of relativity 
would provide a mechanism for the construction 
of a unified theory of gravitation and electro- 
magnetism and this was shown to be rigorously 
the case by Klein. The fact that no physical sig- 
nificance could be given to the fifth dimension 
was a considerable drawback to the Kaluza-Klein 
theory and this objection was removed in the 
projective relativity theory' the formalism of 
which closely parallels that of the five-dimen- 
sional theory. In 1931 Einstein and Mayer? intro- 
duced a new unified field theory that was also 
based upon the Kaluza model but, like the pro- 
jective theory, was four-dimensional; this theory 
has been further considered and modified by 
Schouten and Van Danzig,’ and also by Veblen 
in a forthcoming book. It is known that the 
various theories are closely related and in the 
present paper it is proposed to discuss certain 
formal relations between them.‘ We shall be con- 


1See O. Veblen, “A Generalisation of the Quadratic 
Differential Form,” Q. J. Math. (Oxford) I, 60-76 (1932), 
and O. Veblen and B. Hoffmann, “Projective Relativity,” 
Phys. Rev. 36, 810-822 (1930) where references are given. 

* Einstein and Mayer, Sitz. der Preuss. Akad. (Berlin) 
XXV, 541-557 (1931), and also XII, 130-137 (1932). 

*Schouten and Van Danzig, Proc. Roy. Acad. (Amster- 
dam) XXXIV, 1398-1407 (1931), and also Zeits. f. Physik 
78, 639-667 (1932). ¥ 

‘ Most of the work of the present paper was performed 
before the appearance of the second of the Schouten-Van 
Danzig papers and some degree of overlapping will be 
apparent in the present paper. 


cerned here solely with the first of the Einstein- 
Mayer papers cited; the second carries the 
theory beyond the stage it attained in the first of 
these papers but the physical significance of this 
advance has not yet been investigated. 


§2. ON THE ALGEBRAIC FORMALISM OF THE 
EINSTEIN-MAYER THEORY 


We begin by considering §1 of the first Ein- 
stein- Mayer paper in terms of a slightly modified 
formalism which depends upon the replacing of 
the singular transformation there used by a non- 
singular one. We shall refer to the Einstein- 
Mayer theory as E—M and to equations occurring 
in their first paper upon it as E-M(2), etc. The 
projective theory will likewise be referred to as 
P.R. and equations in Projective Relativity as 
P.R.(2), etc.; we shall deal with P.R. only for the 
case @=1. Greek suffixes will always take on the 
values 0, 1, 2, 3, 4 and Latin suffixes the values 
1, 2, 3, 4. We shall have to deal with two types of 
reference system; for the one we shall always re- 
serve the letters a, 8, y, ---; a, 6, c, «++, and for 
the other the letters \, yw, v, -++;1, m, nm, +++; to 
distinguish the null suffix in the \-system from 
that in the a-system we shall denote the former 
by @. 

To avoid complication we prefer to deal with 
the projective space in terms of its five-dimen- 
sional representation. Let x* be the coordinates 
of such a representation of which x’ are the ordin- 
ary coordinates of space-time; furthermore let all 
geometric objects be such that their components 
are independent of x®, and let us decide to keep 
the factor x® unchanged in the \-system. 
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Let g,, be such that 
2ru= Ory; (1) 


by means of the ‘‘mixed’’ tensor y,.* we define 
alternative components for all geometric objects, 
and thus for the metric, g,,, we have the a-com- 
ponents, gag, given by 


Lap YaVp"2ru= Val ¥a"Zim t+ Va’ Vs". (2) 


If we set 
vo =Aa (3) 


the restriction that x’ be invariant means that A, 
behaves like an a-vector and Eq. (2) becomes 


Lap= Val ¥p"Zim tA aA g (4) 


which is identical with E-M (13b). 
Since gy¢ = 5,9 we have from gyeg = 5" that 


go = 60: (5) 
then, defining y;* and A®@ as the result of raising 


and lowering indices of ya’ and Aq with appro- 
priate g’s as in E—M, we have 


yalAe= gehy.'A = g¥=0 (6) 
which is E-M (10). We have also, at once, that 
LapA*AF= go= 1 (7) 


which is E—M (11). 
Let b be a vector; then }), = 7%), or 


bi=yi%ba (8) 
and 

bo=A*b,; (9) 
Eq. (8) is E-M (14), (the second E—-M (14)!). 
Again, the inverse transformation is 

@ =e 7d* = yi%b'+ Ab? 
and using (9) and (5) this is 
be=¥,%b'+A%A bby 
=yi%b'+pA%, (p=Ag,b*) 

as in E-M. 


(10) 


§3. COVARIANT DERIVATION 


Denoting covariant derivation with regard to 
a connection [ by means of a solidus let us re- 
quire that 


(11) 
(12) 


2ru/» =0 
and also that 
LaB/y = Q; 
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it is a well-known result that (11) gives 


d 
r= | +£*(Syre-+Snet Son) (13) 
MV 


where are the Christoffel symbols formed 


from the g,,, and S,,, is an arbitrary tensor anti- 
symmetric in the second and third indices, 
Since gy¢ = dye and the gy, are independent of x® it 
is at once found that 


d l 
= 5)5,'"5,", 
MY mn 


where ‘| are the Christoffel symbols formed 


(14) 


from the gim. 
Instead of taking the general I given by (13) 
we shall impose the restriction that 


Srur= 50.5,,7"5," Fm ny 
where, necessarily, | 
(15) 
Fn _ Fen | 
and is an affine tensor; 


this is an invariant restriction since @ is a scalar 
suffix. 
As a result of (15) we now have 


ria 


= ’ Tn0' = Ton! = Fo." T im? = Fin; | 
| mn 


(16) 
the rest being — 


The conditions (12) will not suffice to give the 
law of transformation from T,,* to I's,*; we impose 
the following condition which is consistent with 
the condition that (12) be a result of (11); let d 
be any vector and )),, its covariant derivative, 
and let b,,, be the covariant derivative of the cor- 
responding b,; we require that 


Datu=VardrIu3 


the law of transformation from I,,* to T's,% given 
by this condition is then 


0 " 
Ye". (17) 





V7'T ap? = ve T,’+ 
Ox# 


By using (3) and (16) this becomes 
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2D ¥= ; are F,,°"+ Q o 
Vy tam'= Ya lm| alm ax” Ya» 
7) 
A Tam? = Ya! Fimt+ Aa, 
ox™ ] 
and 
V7"T a0? = Val Fr", Al ao?=0. (18b) 


The Eqs. (18b) serve to define, in terms of F and 
y, the quantities I'a46” which, like the I,«, do not 
enter the formalism of E—M explicitly; thus the 
present I'a6” and I’,»* do not involve any quanti- 
ties that are not contained in the other compo- 
nents of I. 

By using a semi-colon to denote E-M co- 
variant derivation and identifying the present 
T'gm* with that of E—-M, (18a) become 


Ye": m@Aal”.=, Aa: m— —¥a' Fim; (19) 


which are E-M (2) and (27), upon which the 
whole of the subsequent work in E-M is based. 
The Eqs. (19), together with (18b), are merely 
an expression of the relation 
Van =0 (20) 


which is thus g,°,,=0 in a different mixed refer- 
ence system. 


§4. TRAJECTORIES 


The E-M trajectories may be defined as the 
curves generated by the parallel displacement of 
the vector dx*/dt along the corresponding 
direction dx'/dt, the displacement being according 
to the connection 1; for, such curves are given by 


dx™ {dx 
(2) 
dt \ dt] ym 


d?x/dt? +-T ym(dx™/dt)(dx’/dt)=0 


(21) 


i.e., by 


or, using (16), by 


dx™ dx! d?x 
na QO —=(), 


dt dt — dt? 


—— 


dt? 


dt dt 


mn 


dx! | l |= dx" 


Thus from the second 


dx*/dt=p, (p=const.) 
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and the remaining equations become 


dx™ 
_ pF m9 


dt dt dt 


dt? 








d*x!' l \— dx" 


mn 


which are E—-M (28). 

The above treatment is essentially the ex- 
pression of the E-M treatment in terms of the 
\-system only and is formally the same as that of 
Schouten and Van Danzig. 


§5. THE Fre_p Equations E-M (47) 
The field Eqs. E-M (47), namely 


Fim, at Pun, it Fat, m=O, E-M (47) 


have no counterpart, as field equations, in the 
Kaluza-Klein theory, in P.R. or in the theory of 
Schouten and Van Danzig, though in the last 
they enter implicitly as a formal assumption; in 
all these theories the above field equations are 
identically satisfied. They imply that there 
exists a four-vector, say B,, arbitrary to within an 
additive gradient and such that 


Fim = 0B,/dx" —OB,,/dx'. (22) 


Let us specialize the y’s by imposing the re- 
quirements that 


Ayo=1, va'= 56,', (23) 


and let us consider the quantities g.s and 
I's47=V's,%yy7’ when this specialization has been 
made; the restrictions (23) imply that the a- and 
\-systems are now related by a non-holonomic 
change of factor, the space-time coordinates 
remaining unchanged. 

The antisymmetric part of a connection does 
not behave as a tensor under changes of reference 
system of the present type and therefore al- 
though I,,* be unsymmetric it will not necessarily 
follow that I's,* is unsymmetric also; the law of 
transformation of I is 


YaT'3,°= Puy s!¥ 7’ +7 7"(Ovs*/dx") (24) 


and the condition that I's,* be symmetric is 


therefore that 


a a 
(Tu—T a) ye'¥y’ +1" — v8\— ve" — yV=0, 
x Ox" 
(25) 


which, by (23) and (16), gives merely 
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2 Finndg”by"+64"(0A p/dx™) — 5a"(0A y/dx™) =0 
or 
2 Finn= a(A adm") /Ox" = (A aOn*)/Ox™. (26) 


This can only be satisfied if E-M (47) is satisfied; 
thus the equation E-M (47) may be interpreted 
as the necessary and sufficient condition that it 
be possible so to choose the Ag, the ya’ being 
restricted according to (23) and the quantities 
T'go*, Tye being given by (18b), (16), that the 
connection I'g,* be symmetric; and in this case 
Fran is given by (26). 
Moreover the gas, by (4), are now given by 


Zab= 5a 5o"Lmn tA aAb, £a=Aa, Z00= Ao= 1, (27) 


and, since the physical significance of the gn» is 
that of being the gravitational potential, and of A. 
that of being, by (26), the electromagnetic poten- 
tial, the present gas has both the form and the 
physical significance of the yas of P.R.® except for 
a difference in sign between the definitions of the 
mn of P.R. and the F,,, of E-M.® Now the intro- 
duction of the F,,, in the E-M theory does not 
fix the sign to be attached to it and the equation 
E-M II could as well have been written 


Ta; «™ —A'Fxq ’ 


since the field equations contain the F,,, either in 
the form of a homogeneous quadratic function, as 
in E-M (46a), or else, as in E-M (46b) and 
E-M (47), as 

f= 0, 


where f is a homogeneous linear function of the 
Fn; and therefore so far as the field equations 
are concerned the sign of F,,, is immaterial. 

Thus, if we consider the F of E—-M to be re- 
placed by its negative, and if we satisfy E-M (47) 
identically through the requirement that I's,* be 
symmetric when the y’s are given by (23), we 
have a one-to-one correspondence between the 
basic formulas and field equations of P.R. 
and E-M. 


§6. THE FIELD Equations E-M (47), 
(CONTINUED) 


In the preceding sections we exhibited a rela- 
tionship between the E—M theory and P.R. on 
the basis of a special choice of the quantities I’,»* 
that are not completely determined by the E-M 


theory, and this relationship gave to the vector A 
the physical significance of being the electro- 
magnetic potential, the tensor F being the curl 
of A;in the present section we wish to discuss the 
relationship between F and A more fully than 
was done in the previous section, where we were 
concerned with deriving the relationship between 
P.R. and E-M. 

The equations E~M (47) imply that there is a 
four-vector, which we shall write as B,6,,*, such 
that 


Fruin=90(B adm) /0x"—O(Babn*)/Ox™; (22a) 


it is evident that if we keep the special values 
that were assumed for the I',»* in the previous 
sections we can find a set of y’s, corresponding to 
a non-holonomic factor change, such that the 
new metric will contain the B as the electro- 
magnetic potential in the manner of P.R.; the 
factor change in question is merely 


ye" = 50”, Va’=Ba, (Bo=1), 


(23a) 


i.e., (23) with B written in the place of A; from 
the work of §5, if we continue to write B for A, 
we have at once that the new I,,* is symmetric. 
Thus we have proved a theorem closely related to 
that of §5; in §5 we started with A as given and 
satisfied E—-M (47) by its means through the con- 
dition that I's,* be symmetric. In the present sec- 
tion we started with any solution of E-M (47); 
this solution defines a vector B (to within a 
gradient) and it is shown that this vector B can 
always be made to play the part played by the 
vector A in §5, the I's,* at the same time being 
made symmetric automatically. 

The above was on the assumption that the I’, 
were chosen in a special manner not forced upon 
us by the formalism of E-M; let us now drop this 
restriction upon I’. The E—M theory is based es- 
sentially upon the Eqs. (19); from the first of 
these it is seen that we must necessarily have 


l 
Pan! = | > Tom'=Fm?; and Ton’=0 
mn 
and from the second of them we must add to this 
list 

Tim? = Fim. 


5 Compare P.R. (3.1). 


6 Compare P.R. p. 815 with the Eq. (26) of the present 
paper. 
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Thus the values we have used for the components 
I,» were forced upon us; it does not follow, 
however, that since the I,’ are not contained in 
(19) they may be taken as arbitrary; for their 
values depend on the values of the components 
of the tensor S,,, of (13) and this S must give the 
I,, correctly. It is not difficult to find to what 
extent the quantities I’, are arbitrary, and for 
this purpose it is unnecessary to consider ex- 
plicitly the alternative forms of S that may be 
used whilst still keeping the old values of I’,,.; 
for we need merely express the condition 


2ru/y a 0, 


which must be satisfied by I’,,, of which I’,, are 
already determined. The only equations in which 
Ie appear are 


£ru/0=0, 
which, since g is independent of x°, gives 


Loul no? + orl yo’ = 0; 
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and, remembering that gy» = 4,9, this at once gives 


P'mno’=T'oe°=0, gmel'se’+2rel'me*=0, (28) 


which limits the arbitrariness of I’,9. 

We may now answer the question, what re- 
strictions do the field equations E-M (47) place 
upon I; evidently the most general I that is 
possible if E-M (47) are satisfied is given by 


l 
I mn = Me ’ Ton! = F,.**, Tim? = Fim; T'mo'= g"Kmey 


where 


Fim= (0B,/dx" —OB,,/dx') and Km.s= —Kam, 


and all the components of I’,,* not included in the 
above are zero. 

The Schouten and Van Danzig interpretation 
of E—M implies that the antisymmetric tensor 
Kim is a multiple of the curl of the vector B, 
and the above would seem to show that such an 
assumption goes beyond the actual content of 
the E-M theory. 
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Neutrons in the Nucleus. I 
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It is supposed that the nucleus consists of a-particles, 
neutrons and zero or one proton, instead of the former 
structure scheme, a-particles, electrons and 0, 1, 2, 3 
protons. The calculated mass defects Am of the isotopes of 
a single element show an exact linear increase of Am with 
the mass number (instead of the fluctuating function in 
the former structure scheme). The Am-difference of the 
isotopes of each element is rather constant throughout 
the periodic system for even elements. The absolute values 
of the mass defects give a potential energy curve decreasing 
uniformly with increasing atomic number (instead of 
showing a minimum at the element No. 50 in the former 
structure scheme). In the region of the radioactive sub- 


stances the energy curve seems to turn up again. All of 
this helps the model that the isotopes of an element differ 
only by the number of neutrons incorporated by them, 
the binding energy per neutron being rather constant 
throughout the periodic system and about 0.009 in mass 
units for even elements. In several cases the isotopes give 
a picture of the construction of a complete outer shell of 
8 or 12 neutrons and indications that also in other cases 
there is a shell structure of the neutron arrangement, but 
none of the a-particles. The empirical rule that odd 
elements have no isotopes with even mass numbers can 
be explained, together with its exceptions, by help of the 
neutron structure scheme. 





OLLOWING Heisenberg,' who has explained 

several facts about the existence and sta- 
bility of isotopes by help of neutrons and 
protons, we assume that the building stones of 
the nucleus consist of a-particles, zero or one 
loose proton, and neutrons. It may be of interest 
to support this model from some other points of 
view throughout the periodical table, stressing 
more than does Heisenberg the rdle of the 
a-particles in the whole nuclear structure. 


I. IsoTOPES OF A SINGLE ELEMENT 


There are only a few elements with several 
isotopes whose atomic weights have been meas- 
ured to three decimals, so that one can calculate 
their mass defects with sufficient accuracy. 
These are Kr 36 (6 isotopes), Sn 50 (11 isotopes), 
Xe 54 (9 isotopes), Hg 80 (9 isotopes). Now the 
former theory interprets the structure of the 
Xe-isotopes Z = 54 in the following way: 


Xe"4=31la+8 el, 
Xe'*>, missing, 
Xe! = 31a+10 el+2 pr, 
Xe!’, missing, 
Xe§ = 32a+ 10 el, 
1W. Heisenberg, Zeits. f. Physik 78, 150 (1932). 


Xe”? = 32a+11 el+1 pr, 
etc. 


assuming four different numbers of a-particles 
from 31 to 34. This scheme gives the following 
values for the mass defects of the nine isotopes 
(compare the mass-spectrum of Xe in Table II): 
0.133, —, 0.150, —, 0.137, ?, 0.154, 0.162, 0.141, 
—, 0.158, —, 0.154, which are seen to fluctuate 
with increasing mass of the isotope (compare 
Gamow’s Table II of reduced weights with 
He = 4.000). 

Now since Aston’s weights? of the Xe-isotopes 
increase from 123.867 to 135.855 in. almost 
exactly equal steps (passing over the missing 
isotopes) it is rather unsatisfactory to account 
for this linear rise by a fluctuating function. 
One has rather to explain it by a linear sequence 
of equal binding processes. In the scheme of 
neutrons (v) the Xe-isotopes are to be described 
by 


Xe"™4=27a+16v» to Xe™§=27a+28». 


Here the mass defect of Xe'* becomes 27 
< 4.000 + 16 & 1.0075 + 16 XK 0.0005 — mMobserved 
=(.261, representing the energy required to 
split up the nucleus into 27 a-particles and 16 


2(. Gamow, Atomic Nuclei and Radioactivity, Oxford, 
1931. 
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neutrons and break up each neutron into its 
components. The last part could be separated if 
we knew the exact atomic weight of the neutron. 
If one calculates in the same way the mass 
defects of all Xe-isotopes, one gets the results: 
0.261, —, 0.279, —, 0.297, ?, 0.315, 0.324, 0.333, 
—, 0.351, —, 0.369, indicating a linear increase of 
the defect by equal steps of 0.009 each, parallel 
to the observed linear steps of the atomic 
weights. This result is based on the assumption 
that all of the Xe-isotopes consist of 27 a- 
particles and 16 to 28 neutrons, the energy 
required to split off each neutron from the 
nucleus of Xe and to break it up into a proton 
and an electron being exactly constant, equal to 
0.009 in mass units. 


II. IsoroPpEs oF VARIOUS ELEMENTS 


If one carries out the same calculations for 
the isotopes of other elements insofar as their 
atomic weights are sufficiently determined, one 
obtains in all cases the same exact linearity of 
the mass defects, parallel to the fact that the 
observed atomic weights of the isotopes form a 
linear set. We have then merely to report the 
constant intervals between subsequent mass 
defects for various elements, among them 0.009 
for Xe. These are given in Table I for those 


TABLE I. Mass defect differences. 


ORES eee te ae a > ae 0.009 
Ring aneaewaen Sere 0.008 
| re 0.009 ECS iiwaxss'a'eces 0.009 
irae i cs en 0.008 oe: 0.013 
SEE Ene Si > 0.011 
_ EE ae Se PG $e dea owns 0.010 


elements with more than one isotope, whose 
atomic weights are known up to 3 decimals. 

We interpret then the values of Table I as the 
energy in mass units to split off one neutron and 
to break it up. The atoms of even charge 
number display in general many more isotopes 
than the odd elements. For even elements the 
binding energy per neutron seems rather constant 
throughout the periodic system except for Ne. 
The binding energy per neutron for odd elements 
is more variable and on the average a little 
larger than for even elements. This may be due to 
the one extra proton which increases the attrac- 
tion of the neutrons and disturbs the symmetry 
of the neutron arrangement. 
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III. Nuctear STABILITY THROUGHOUT THE 
PERIODIC TABLE 


If one plots the mass defects according to the 
old theory (a-particles, loose electrons and 
0, 1, 2, 3 loose protons) one obtains for the 
potential energy of the nuclei a distribution 
having a minimum in the range of element 
number 50 (compare Gamow’s Fig. 16). This 
minimum is a serious difficulty, because it means 
that the elements above 50 should be unstable, 
while in fact they are unstable above 80. If 
however we plot the mass defects according to 
the new scheme (a-particles, neutrons and 0 or 1 
proton), then we obtain the mass defects given in 
Fig. 1. Each mass defect dot is marked with the 
number of neutrons contained in the isotope. 
The figure shows a general decrease of the po- 
tential energy with increasing atomic number. 
Among the radioactive substances the weights of 
Ra, Th and U are determined to only 2 decimals, 
in which the last one is not certain. So we have 
drawn dashes instead of dots in the region of 
their probable mass defects. Among the lighter 
elements the measured weights are rather in- 
complete, and sometimes only one of several 
isotopes could be used in our figure for that 
reason. The potential energy seems to have a 
minimum just in the region of the radioactive 
substances, where stable elements cease to exist. 
That our Am increases uniformly while the 
curve of Gamow’s book bends up after 50 is due 
to the fact that in the old scheme the number of 
a-particles is larger than in the neutron theory. 
For instance Hg 80 with m= 196 to 204 has in the 
old scheme 49 to 51 a-particles, in the neutron 
scheme only 40. 


IV. EXIsTENCE OF Opp CHARGED NUCLEI 


The rules of Harkins concerning the relative 
abundance of the elements and their isotopes 
throughout the periodic table read in our scheme: 
Elements with a loose proton (odd charge 
number) are much rarer than elements without a 
loose proton (even charge number). Furthermore 
if we consider first only even elements we find 
that isotopes with an even number of neutrons 
are more frequent than those with an odd number 
of neutrons. Missing isotopes in most cases belong 
to odd neutron arrangements. This fact cannot be 
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explained by energy balances, since we learned in 
Part I that the successive incorporation of one 
neutron after the other produces a linear increase 
of the mass defect without any break, at least 
for the neutrons added last. We have then to 
take for granted the preference of even arrange- 
ments of neutrons without an explanation by 
energy balances. 

As to the elements of odd charge number one 
may wonder why they exist at all, both with even 
and odd numbers of neutrons. For if processes of 
B-emission connected with a transformation of 
one of the neutrons into a proton are taken into 
account, then the original loose proton and the 
new proton should associate with two neutrons 
to form an a-particle, and the odd element 
should pass over into an even one. However, 
though the final product has one a-particle more, 
it has three neutrons less. Now we learn from 
Table I that the energy gain 0.032 of building up 
an a-particle is about the energy loss of sepa- 
rating three neutrons and one 6-particle from the 
nucleus; so these two effects may approximately 
balance each other, and we expect that minor 
secondary causes shall give the decision whether 


an element of odd charge number is or is not 
stable. 

Now experience shows as a rule that odd 
elements have isotopes with odd mass numbers 
only, or expressed in our scheme: Odd elements 
have no isotopes with an odd number of neutrons. 

We can explain this rule in the following way. 
The above-mentioned 8-emission transforms an 
odd arrangement of neutrons into an even one. 
If there is only a small preference of even 
neutron arrangements before odd ones, this may 
give the decision mentioned above that the 6- 
emission in this case always will happen, hence 
odd arrangements of neutrons in odd elements 
will be unstable. In the reverse case, if an even 
number of neutrons has only a small tendency to 
remain even instead of turning odd, this will 
give the decision that the B-emission is barred, 
and the isotope is stable. 

This consideration does not hold if the odd 
number of neutrons is just unity. For then a 
B-emission would leave no neutron to form a new 
a-particle. In fact, there are four exceptions to 
the rule, namely, the existence of the stable 
isotopes, H2, Li®, B, N'*. But these are just four 
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cases in which our structure scheme gives one 
neutron only, namely: 


goer Oa+1 pr+1p 
Mn seekies la+1pr+lyp 
ag ereeee 2a+1 pr+1y 
PP énestew 3a+1 pr+ip 


Thus from the point of view of the neutron 
scheme the four exceptions prove the rule. The 
same holds for the exceptions among the radio- 
active elements giving §-rays. 


V. NEUTRON SHELLS 


It is very tempting but precarious to conclude 
from the table of the isotopes known up to now 
details of the arrangement of the neutrons in the 
nucleus. The existence of closed shells of neutrons 
is suggested by the fact, so striking at the first 
view of an isotopic table, of similar isotope 
sequences repeated in different elements. Con- 
sider for instance the following cases of mass- 
spectra given in Table II. Missing isotopes are 


TABLE II. Mass-spectra. 














Ge Se Kr Mo Ru | Hg | Pb Po Sn Xe 
32 34 36 | 42 44 | 80 | 82 84 | 50 54 
a ae em. Pare ee 
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| — 
| | |= - 
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6» 6» 6v | 8» 8» | 36y | 42» 42» | 12» 16» 








indicated by dots. The second line means the 
charge number of the element, the last line of the 
table indicates the number of neutrons contained 
in the isotopes of lowest mass. We interpret this 
number as the total number in complete or 
incomplete inner neutron shells, to which neu- 
trons are then attached step by step until an 
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outer shell is completed. The examples of Table 
II we interpret as the incorporation of an outer 
shell of 8 or of 12 neutrons. As to the dots we 
observe that usually odd arrangements of neu- 
trons are missing. With Se, Kr, Mo, Ru and zlso 
Sn and Xe, the isotopes with one neutron in the 
outer shell and the isotope with one neutron 
lacking in the outer shell are not observed. 
Whether this is essential or merely due to small 
intensity may be judged from the case of Hg. 
Here Aston’s first experiments missed these two 
isotopes, but in a recent paper he reported them 
present in small intensity. So if one looks only 
for the possible population of neutron shells and 
for quantum rules controlling them, one should 
consider the missing isotopes more as an intensity 
problem. 

Until the table of isotopes is completed by a 
table of their exact atomic weights to 3 decimals, 
it is rather precarious to generalize the idea that 
the smallest isotope represents all neutrons in inner 
shells, and the largest isotope an additional com- 
plete outer shell. This idea may be justified for 
the group of elements in Table II, since for those 
of them which are measured the mass defects 
certainly increase in equal steps, with the step 
magnitude 0.008 and 0.009. From mere energetic 
considerations one should expect that the in- 
corporation of the neutrons takes place with 
decreasing energy steps, until the last neutron is 
bound with nearly no energy and bars the 
incorporation of more neutrons. This helps the 
conclusion that the neutrons are arranged in 
shells like the outer electrons, subject to Pauli’s 
exclusion principle and Fermi statistics.* At the 
same time one expects for the a-particles, since 
they consist of an even number of components, 
Bose statistics and no shell structure at all. 


Note added in -proof: After this paper was sent in, I 
received an article of E. N. Gapon‘ being rather identical 
with the contents of Sections I, II and III of the present 
paper. Compare also with D. Iwanenko® who first suggested 
the nuclear model used here. One may consider, then, 
these parts as an introduction into Part II. 


3J. H. 
(1932). 

4E. N. Gapon, Zeits. f. Physik 79, 676 (1932). 

5 D. Iwanenko, Comptes Rendus 195, 439 (1932). 


Bartlett, Phys. Rev. 41, 370 (1932); 42, 145 
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Neutrons in the Nucleus. II 


ALFRED LANDE, Mendenhall Laboratory, Ohio State University 
(Received February 10, 1933) 


In Part I we came to the conclusion, that in several cases 
the lightest isotope of an element represents an arrangement 
of neutrons in more or less closed shells, that the heavier 
isotopes indicate the attachment of one neutron after the 
other, until with the heaviest isotope a new outer shell is 
completed. In the present paper we account for the number 
of neutrons in the nuclei by the following model. Up to 
Z=16 the neutrons populate a first shell with at most two 
neutrons. From Z=17 to 28 an additional shell of four 
neutrons is allowed. From Z =32 to 36 we see the presence 


of neutrons in a third shell of eight. It is necessary to 
assume that between Z=36 and 40 the inner shells are 
filled in with two more neutrons, and the same between 
Z=50 and 54. Evidence for this filling in is seen in certain 
significant features of the binding energy. It is pointed 
out that the intensity distribution of the isotopes points 
to their stability with respect to a and 8 transformations. 
This can be used to give a more rational basis to the art 
of predicting missing isotopes. 





I, 


SING the scheme ‘“nucleus=a-particles 
+neutrons+zero or one free proton”’ in the 
case of an even or odd charge number, we have 
plotted in Fig. 1 the number of neutrons in the 
isotopes against the charge number Z=1 to 92 
of the elements as far as known. The result can 
be described as follows (compare with the sche- 
matic diagram of Fig. 2). Up to Z = 16 the number 
n of neutrons varies between 0 and 2. From Z=17 
to 28 the number ” runs from 0 to 2 to 6. From 
32 to 36 we have m between 6 and 14. From 40 
to 42 m is found between 8 and 16. The element 
Z=50 has from 12 to 24, and in Z=54 » runs 
from 16 to 28. The significant feature of this dis- 
tribution is found in the fact, that the upper and 
lower limit of the m-range remain constant 
through a series of Z-values, and then jump up 
suddenly by several units. In some cases the 
former m-maximum becomes now ”-minimum. 
We can account for this numerical arrange- 
ment by the model of Fig. 2. Up to Z=16 the 
neutrons populate a first shell with at most 2 
neutrons. From Z=17 to 28 an additional shell 
of 4 neutrons is allowed. In Z = 32 to 36 we see the 
presence of neutrons in a third shell of 8. It is not 


1 The odd elements give mass defects a little above our 
smooth curves »=const., if we calculate Am by the scheme 
nucleus = a-particles+neutrons-+ proton, and a little below 
the curves in the scheme ‘“‘neutrons+ protons.”’ This shows 
that one of the a-particles has a looser structure inside 
the nucleus than when isolated. 


meant that these shells are already complete. On 
the contrary we have to assume (verification in 
Section II), that between Z = 36 and 40 the inner 
shells are filled in with 2 more neutrons, so that 
the third shell of 8 neutrons rests now on 8 
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Fic. 1. Number of neutrons in the isotopes as a function of 
atomic number. 
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Fic. 2. Idealized schematic diagram corresponding to Fig. 1. 
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instead of 2+4=6 inner neutrons. With Z=50 
we meet a fourth shell of 12 neutrons built on the 
former three shells diminished by 4 neutrons. 
But with Z=54 these 4 neutrons are filled in 
again (see Section II), and the fourth shell is 
based now on 16 inner neutrons. The observa- 
tions of isotopes for higher Z are too incomplete 
to draw conclusions from them. 

Altogether we obtain for the elements up to 
Z=54 a table of their capacity for neutrons, 
(Table I) not made use of in the case of ‘“‘miss- 
ing” isotopes. 


TABLE I. Neutron shells. 














P 4 Neutron shells 
3-16 
20-28 24 
32-36 248 
40-44 268 
50 264 12 
54 268 12 








Our taking recourse to the assumption of inner 
shells being filled in seems to render the result 
very hypothetical indeed. It is with the help of 
energy considerations that this structure scheme 
can be verified in an independent way. 


Il. 


In Part I we calculated the mass defects Am 
of the isotopes derived from their measured 
atomic weights according to the scheme nucleus 
=a-particles+neutrons+zero or one proton. 
If we join by a curve the mass defect dots in 
Fig. 1 of Part I belonging to a certain constant 
number n, these curves will represent the 
energy of binding a constant number m of neu- 
trons to a varying number Z/2 of a-particles, 
considering only the even! elements. Although 
there are not many dots on each particular curve 
n=const., that are obtained from direct measure- 
ments of the exact atomic weights, we know from 
Part I, that these curves run nearly parallel at 
constant distances from each other throughout 
the periodic system.? This enables us to extra- 
polate these curves in a rather wide range with 
good reason. The result is shown for the even 


? This fact was first found by E. N. Gapon, Zeits. f. 
Physik 79, 676 (1932). 
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Fic. 3. Mass defect for constant m as a function of Z. 


elements in Fig. 3. In general the mass defect 
Am increases with increasing number of a-par- 
ticles attached to a constant number » of neu- 
trons. Only between Z =36 and 40 and between 
Z=50 and 54 the curves are descending. These 
are just the ranges, where the numerical distri- 
bution of the isotopes suggested the assumption 
of neutrons being filled into inner shells. We have 
to interpret the descending parts of the curves 
n=const. as a sign, that the uniform process of 
attaching a-particles to a fixed number of neu- 
trons is interrupted here. Namely a-particles are 
apparently bound in these ranges looser than in 
the preceding ones, as if they were attached more 
toward the outside of the nucleus, and some of the 
neutrons were shifted to inner shells. One could 
try to explain the stability against a-emission 
on a descending part of the Am-curve by the 
assumption that an emission of an outside a- 
particle would leave the nucleus in a state of 
excitation with a Am smaller than the measured 
Am of the stable product after such an emission. 

The extrapolation of our mass defect curves’ 
n=const. into the unexplored range Z = 54 to 80 
shows that corresponding descending parts must 
occur here, too, in order to join the curves onto 
the measured mass defect dots of Z=80. 


Ill. 


The isotopes of one element, though differing 
only by the number of neutrons in the outer 
shell, cannot be supposed to be in a state of 
mutual equilibrium. Else we should expect to 
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find a maximum abundance for a certain isotope 
and a uniform decrease of the intensity to both 
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Instead we find that the relative abundance var- 


ies in large steps up and down, as in the instance 





heavier and lighter isotopes of the same element. of Z=54: 
m=112 113 114 115 116 117 118 119 120 121 #122 «123 =©«©124 
Intens. 0.1 0 0.1 0 2.3 27.1 4.2 20.7 26.4 0 10.3 0 8.8 
Only if we consider the 4 families m =4N, 4N+1, 4N+2, 4N+3 separately: 
112 116 120 124/113 117 121/114 118 122)115 119 123 
0.1 2.3 264 88] 0 27.1 0/02 4.2 10.3) 0 20.7 O 





then we obtain a uniform decrease on both sides 
of the maximum abundance. This indicates of 
course, that the intensity distribution can only 
be understood from the historical point of view of 
a- and £@-disintegration or synthesis. For this 
reason we think that speculations about the 
formation of nuclei by successive addition of pro- 
tons and neutrons have not much significance. 
* For every charge number Z there should be a 
certain number of neutrons possessing the 
largest stability, and the stability maximum 
should run along a uniformly ascending curve in 
the (nZ)-diagram, smoothing out the corners 
of Fig. 2. From this point of view one can under- 
stand the abundance distribution of the isotopes 
of each radioactive family separately. In partic- 
ular the isotopes represented by convex corners 
of the (wZ)-diagram are extremely weak or even 
missing, indeed. 

Although the members of one family have mass 
differences of 4 units for a given charge number, 
the even families 4N and 4N—2 contribute in 
many cases 2 or even 3 members to the isotopes 
of one element. These belong to different branches 
of the family because of forking processes, 
familiar in the case of radioactive elements. The 
‘odd families 4N—1 and 4N—3 are, with a few 


exceptions, represented by only one isotope in 
each element. Forking seems to be much rarer 
here. Consequently there is no great chance of 
predicting missing isotopes in the odd families. 
Furthermore, according to the rule obtained in 
Section IV, Part I, there is no chance to find 
stable isotopes of odd charge number with even 
mass numbers, except for the first elements. 

A method of prediction based on the descend- 
ence theory has to consider two separate cases. 
Either the present abundance of the “‘stable”’ 
elements and their isotopes is the effect of a 
decay process. Then one will predict missing iso- 
topes on the heavier end of the mass spectrum of 
some elements. For instance, one will predict the 
missing isotopes 7 = 6 in Z = 26 and 28, since n =6 
in Z = 24 can be reached by a-disintegration only 
through the two missing isotopes, if one wants to 
avoid too large detours from the intensity max- 
imum. Or the present abundance is the result of 
a- and 6-synthesis. Then one will predict missing 
isotopes on the lighter end of the mass spectro- 
grams of some elements. It would be interesting 
to tell from a more completed table of isotopes, 
whether the present state appears to be the re- 
sult of decay or synthesis. 
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The Vector Model and the Pauli Principle 


M. H. Jounson, Jr.,* New York University 
(Received February 21, 1933) 


The vector model is examined in the light of the quantum 
theory of wave fields. First the distributions of electrons 
among the electronic states are limited to correspond to 
the atomic states arising from one electronic configuration. 
A set of independent angular momentum vectors 1, l:, 
***Si, Se*** is found, in which one orbital and one spin 
vector are associated with each group of electronic states, 
rather than individual electrons, having given values of the 
quantum numbers m and /. Any other vector T may be 
decomposed into a sum of independent vectors ti, one 


associated with each of the same group of electronic states. 
The usual commutation relations obtain among the 
vectors t;, 1; and s;. On further limiting the distributions to 
correspond to atomic states arising from electronic con- 
figurations containing no equivalent electrons, it is shown 
that the matrix of 2,f(l°, sx°) (an operator in coordinate 
space in which f is any algebraic function of the components 
of 1 and s) can be constructed from the vectors 1], le, «+ +s:, 
So, ++ according to the usual laws of matrix mechanics. 





N previous papers! a method of calculating the 

matrix H,=)>°,a,);-s; for any electronic con- 
figuration has been described. Here it was as- 
sumed that; (1) l- - 1,81: - +S, constitute a set of 
independent angular momentum vectors: satis- 
fying the usual commutation relations. (2) the 
total orbital angular momentum L is given by 


L= DY elk (1) 
and the total spin S by 
S= > Sk ° (2) 


This method enables one to obtain the matrix H, 
in any representation described by giving definite 
values to J? (J=L+S), L* and S& and by as- 
signing a definite structure to the 1 vectors and 
to the s vectors. It could easily be generalized to 
give this matrix in a representation described by 
any definite structure of the vectors 1,---l,s; 

**Sn. 

In the above references the results of applying 
this method were directly used to obtain the 
matrix of the spin-orbit interaction for various 
electronic configurations. Obviously some proof 
is needed to justify this step for it is not at all 
self-evident that any such set of vectors as 
1\-+-1,8::-+S, can be found when the Pauli 


* National Research Fellow. 
1M. H. Johnson, Jr., Phys. Rev. 38, 1628 (1931); 
ibid. 39, 197 (1932). 
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principle is correctly taken into account. We are 
naturally led to a more general question. In view 
of the Pauli principle which requires that all 
electrons be identical in nature, what meaning 
can be ascribed to the vector model of an atomic 
configuration? The point has already been con- 
sidered by von Neuman and Wigner’ who, in 
their discussions of vector coupling, take account 
of the symmetrization of the wave functions 
required by the Pauli principle. More recently 
Shortley® has justified the specific applications 
made in reference 1 in a similar way. Nevertheless 
we propose an examination of this question from 
a different point of view which is more directly 
concerned with the vector model itself. Whereas 
these authors use wave functions derived from an 
analysis of vector coupling and then sym- 
metrized for the Pauli principle, we shall employ 
a formulation of the problem embodying the 
Pauli principle at the outset and shall find therein 
the constituents of a vector model. 

For this purpose the quantum theory of wave 
fields‘ is convenient. Then if 


X= DieX%(xx) (3) 


2 J. v. Neumann and E. Wigner, Zeits. f. Physik 49, 73 
(1928). 

3G. H. Shortley, Phys. Rev. 40, 185 (1932). 

4P. A. M. Dirac, Proc. Roy. Soc. (London) All4, 243 
(1927); P. Jordan u. O. Klein, Zeits. f. Physik 45, 751 (1927); 
P. Jordan, ibid. 766 and 44, 473 (1927); Jordan u. Wigner, 
ibid. 47, 631 (1928); Fock, ibid. 75, 622 (1932). 
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is an Hermitean operator in coordinate space,° it 
will be represented by 


X =D pra p*a,(p| X°| 1) (4) 


in the space of the N,’s (VV, =the number of elec- 
trons in the state 7). Here 


(p|X°|r)= f U p*(x)X(x)a-(x)dx= (r| X°| p)*. (5S) 


JOHNSON, 


JR. 


Similarly if 
W= D> W(x, x1), W (xx, x1) = W(x, Xx) (6) 


is an Hermitean operator in coordinate space, it 
will be represented in the space of the N,’s by 


W= pps pars pa @'a.a,(pq | w® | rs), (7) 


where 





(pq| W®|rs)= fueron at (X2) W(x1, 2) Ur (%1) (Xe) dx dx2= (qp| W®| sr)=(rs|W®| pq)*. (8) 





Here the u’s form a complete set of solutions of a 
central field problem yielding definite values for 
the projections of the orbital and spin momentum 
along the z axis. The summation indices range 
over all the electronic states of this central field 
problem. The a’s are the noncommutative 
amplitudes which operate on the variables N, 
and satisfy the well-known exchange relations 


AptA-+O,0pt=5pr, @pta;++a,tapt=0, 
9) 


a pA, +a,dp,= 0. 
They have the following operational formulation.® 


a,t= N,A;v,= = N,v,A,-= A,(1 — N,) Vry 
(10) 
a,= v,A,N — - A,v,N,= (1 ~— N,) v,A,, 


where A, operating on a function of N, changes 
N, to 1—N, and », is the 
j.1(1—2N;). It is clear that 


sign function 


a;ta,=N,, a,a,+=1—N,. 


(12) 


As a consequence of the exchange relations 
satisfied by the a’s, the N’s commute among 
themselves and have the eigenvalues 0 and 1. We 
shall also need the following commutation re- 
lation which is easily obtained from Eq. (9). If Y 
is an operator similar to X 


[X, YJ=Ypa,*a-(p|LX°, Y°]/r), 


where 


(13) 


5 The superscript ° is used to distinguish operators in 
coordinate space from matrices and operators in the space 
of the N,’s. 

*W. Heisenberg, Ann. d. Physik 10, 888 (1931). 


(p|LX®, ¥°]|7) 
= L.[(p| X*|s)(s| ¥°|r) —(p| ¥|s)(s| X*|r)}. 


We shall use this relation in special cases in 
which at least one of the matrices X or Y is 
diagonal in some quantum number. For example 
X may be diagonal in the inner quantum number 
j but not in the magnetic quantum number. 
Then it is convenient to write X =} ;X; and 
Y=) ;Y;+>;; Yj;. For each X; and Y; the 
indices p, r cover the range of different magnetic 
quantum numbers belonging to the same 7. We 
can calculate [X,;, Y;] in the same way as Eq. 
(13) and the summation over s will then be 
limited to the subgroup of states having the 
same j. As long as at least one of the matrices 
X or Y is diagonal in j, this sum will still be equal 
to (p [X°®, Y°l'r). This fact will be used in 
obtaining Eqs. (18), (19) and (23). 

As the N’s commute among themselves, they 
may all be simultaneously diagonal. A physical 
state of the atom can be described by assigning 
the value 0 or 1 to each N, that is by giving the 
distribution of electrons among the electronic 
states. Thus we might have one electron in the 
electronic state 1, one electron in the electronic 
state 3 and no electrons in all other electronic 
states. The eigenfunction in the space of the 
N,’s for this state will be just the product of a 
number of 6 functions 6(1, N,)6(0, N2)d(1, Ns) 
5(0, N4)---, where each 6 function is analo- 
gous to the familiar spin function (because N,, 
like the spin variable, has only two possible 
values). With a set of such eigenfunctions, the 
matrix of any operator in the space of the N,’s 
may be constructed in the usual manner. On the 
other hand we know that the above atomic state 
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can be equally well described in configuration 
space by the antisymmetric eigenfunction 


Uy(x1) U3(X1) 








ui(X2) u3(Xe) 


The matrix of any operator in configuration 
space can readily be constructed from such a set 
of antisymmetric eigenfunctions. Since the for- 
mulation of a problem in terms of quantized 
amplitudes is mathematically equivalent to its 
formulation in configuration space with anti- 
symmetric eigenfunctions, the matrix of any 
operator constructed from eigenfunctions in the 
space of the N,’s is exactly equivalent to the 
matrix of the same operator constructed from 
antisymmetric eigenfunctions in configuration 
space. 

The problem of obtaining the eigenvalues of 
H, the Hamiltonian operator expressed in terms 
of quantized amplitudes, is equivalent to an 
exact solution of the Schrédinger equation with 
antisymmetric wave functions. Just as it is 
convenient to solve the Schrédinger equation 
considering a limited group of atomic states, so it 
is here convenient to consider together a corre- 
spondingly limited number of distributions of 
electrons among the electronic states. With this 
in mind we divide all the electronic states r into a 
number of subgroups 7 according to the different 
values of the quantum numbers ,/;. In each 
subgroup there will be 2(2/;+1) electronic states 
corresponding to the projections of /; and the 
spin along the z axis. Let 1; be the number of 
electrons in the subgroup 7. 


Mi=¥ (N= pra,ta(p| Mi |r), 


) (14) 
(p| M | 7) = 65,6;", 


where >" ,(7) means a summation limited to the 
states of the subgroup 7 and 46;" is a function 
which vanishes unless 7 is in the subgroup 7 in 
which case it has the value 1. We will consider 
together distributions of electrons among the 
electronic states for which M,, M2, etc., have 
fixed values. These distributions correspond to 
the states of an electronic configuration in which 
there are 1/7, equivalent electrons with ml =ml,, 
M, equivalent electrons with nl = nls, etc. 

Let us see what restrictions the constancy of 
M,, Ms, +++ places on an operator X which may 
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occur in the Hamiltonian. For M; to have a 
definite value it is necessary that 


[X, M.J=0. (15) 


On employing Eqs. (13) and (14) the condition 
becomes 


(p| X°| r)(5i" — 5,7) =0. 


To insure that this equation is fulfilled for all 7 
we must omit all terms from }>°,,a,+a,(p| X°|r) 
except those for which p and 7 lie in the same 
subgroup. Hence 


x= DiXi, X= > pr(1)a p'a,(p | X?| r). 


To see with what matrix components the omitted 
terms in Eq. (16) are connected, consider a,*a, 
when # and 7 are in different subgroups i and j, 
respectively. From Eq. (10) it is clear that 
a,*a, operating on an eigenfunction consisting 
of a product of 6 functions, gives rise to a new 
eigenfunction obtained from the first by re- 
moving an electron from the state 7 of the sub- 
group j and placing it in the state p of the sub- 
group 7. Thus such terms are connected with the 
matrix components between different electronic 
configurations. The omission of these terms is 
precisely equivalent to a first order perturbation 
theory for the Schrédinger equation which 
considers together only those states arising from 
a single electronic configuration. 

We now fix our attention on the angular 
momentum operators. By Eq. (16) we have 


L=S: di, L=D>,(t)a,ta(p|l\r). (17) 


We note that every component of 1; commutes 
with every component of 1;(i#7) as a result of 
the exchange relations satisfied by the a’s. 
From Eq. (13) and the remark following it, we 
have 


(16) 


1, X1,= > p(i)a,ta-(p| 2 XI*\r). 
But 1° 1° =71°; so we find 


1, X1= 11k. (18) 
In a similar way it may be shown 
S=D isi, si=b>-(i)a,ta,(p|s*|r), 
(19) 


s;Xsi= iS;. 


Also we find that every component of s; com- 








630 M. H. 


mutes with s;(i#j) and with all components of 
1. (where k may be equal to 7). 

By this decomposition of the total angular 
momenta we have found a set of independent 
angular momentum vectors ];, ls, «++, $1, Se, **°*, 
one orbital and one spin vector being associated 
with each subgroup 7, which have the correct 
commutation properties. Before the formulation 
of the Pauli principle with antisymmetric wave 
functions one could speak of the m’th electron, as 
distinguished by the coordinates x,, as having a 
certain angular momentum. When antisymmetric 
wave functions are introduced, the possibility of 
distinguishing an electron by its coordinates is 
lost so this manner of speaking is no longer 
correct. Nevertheless vector model consider- 
ations are still valid if we understand by such an 
expression as “the angular momentum of a d 
electron,” an angular momentum vector as- 
sociated with the d electronic states rather than 
with some electronic coordinates. It should be 
observed that we have succeeded in finding 
only two vectors 1; and s; belonging to the sub- 
group 7. If 1/;=3 (say), meaning there are three 
equivalent electrons with ml = njl;, we are justified 
in speaking of ‘‘the angular momentum of the 
shell of three electrons’ but we have not found 
anything that corresponds to the angular mo- 
menta of the individual electrons in this shell. 

It is of interest to find the commutation 
relations of other vectors with the angular 
momenta. Suppose 


T°= rm zt?(x;) 


is a vector in coordinate space which satisfies the 
usual commutation relations 


(20) 


[2.°, t2°]=0, [1,°, t,°]=[#.°, 1,°]=it2. (21) 
Then by Eq. (16) 
T=DViti, ti= DVp-(i)ap*a,(p|t’|r). (22) - 


Here t; is an independent vector in the sense 
that it commutes with all components of t;, 1; and 


s;(i#~j). Eq. (13) yields 
Lei, tyi J= DL pr(t)apta,(p| LI", ty’ ]| 7). 
Then because of Eq. (21) 


cm ty: j= it 2. (23) 


JOHNSON, 
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We therefore see that t;, 1; and s; satisfy the 
same commutation relations as t°, 1° and s°. 

Finally let F°= dor f(k°®, sx°), where f(1°, s°) is 
any function built up by algebraic processes from 
the components of 1° and s°. We decompose F by 
Eq. (16). F=}>0:F;. We now limit the values of 
M; to one or zero which means that we are 
dealing exclusively with electronic configurations 
containing no equivalent electrons. We propose to 
show that with this restriction 


F,=f(li, Si). (24) 


Suppose Eq. (24) is true for two functions G 
and H. Then it is true for their sum (for brevity 
we write g° for g(1°, s®) and h® for A(1°, s®)) for we 
have 


(G+H)=Di(G+H);, 
(G+H):= > p(i)a,ta(p| +h |r) =Gi+H; 


from which the desired result follows immedi- 
ately. Next consider the product 


(GH) = L(G), 


, (25) 
(GH) =D pr(i)apta,(p | gh? |r). 


Because of the limitation of 1/7; to the values 
zero or one, we are dealing with at most one 
electron in the subgroup 7. For this reason any 
operator representing the interaction of two 
electrons in the subgroup 7 must vanish identi- 
cally. Thus, in particular, 

LD vars(i)ap*a gta,a,(p|g®|r)(q|h°|s)=0. (26) 
As this is an essential point in the argument, let 
us consider it in a little more detail. From the 
fact that a,a,=0 we find a,(1— N,) =0, N,a,=0 
or if 7 is in the subgroup 7 a,M;=(1+ M,)a,. 
Now let M; be a variable so that 6(.M/;’, M,)f 
represents a state in which there are 1/;’ electrons 
in the subgroup 7. Here f is some function of the 
N’s in the subgroup 7 which gives the distribution 
of electrons among the states of the subgroup. 


Then 
a,6(M, M;) = (My, M,;+1)a, 
= 6(M,/—-1, V;)a,. 
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Thus the effect of commuting a, through the 6 
function is to decrease the number of electrons in 
the subgroup 7 by 1. If M;=1, the result of 
commuting a,a, through the 6 function is plainly 
zero so that in this case such an operator as 
Y pars(t)apta,ta.a,(pg|A\|rs) always yields zero 
when applied to an eigenfunction in the space of 
the V’s. Returning to Eq. (26), if the order of the 
a’s is changed, it is readily found that 


GilTi= LX pars(t)ap*ara gta,(p| g°|r)(q| h®| s) 
= LD pre(t)aptar(p| h° | s)(s | g° | = (GH). 


It immediately follows if Eq. (24) is true for G 
and H, it is true for their product (GH), Eq. (24) 
is surely satisfied if g° and h® equal any of the 
components of 1° and s°. Consequently it is true 
for any function built up by algebraic processes 
from these components. 
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In particular let 
P= >: Vir) °°, 
(27) 


H;= > pr(t)apta-(p| V(r)I-s°| r). 


Here V(r) is independent of the projections of the 
momentum vectors along the z axis so 


(P| V(P)P-s°|r)=ai(p|?-s°\7r), 


where a; is the integral of V(r) over the radial 
part of the wave functions u for the subgroup i. 
Hence 


IT;= aid pr(t)ap*a,(p| °-s8°|r)=aili-s;. (28) 


Thus we see in this case of an electronic con- 
figuration containing no equivalent electrons we 
are completely justified in building up the matrix 
of the spin-orbit interaction from the vectors 
1,, ls, «++ Si, Se, «+» according to the usual laws of 
matrix calculus. 
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(Received February 21, 1933) 


Heisenberg’s treatment of an almost closed shell is 
extended to electronic configurations consisting of several 
almost closed shells plus other electrons not in the shells. 
In this way Shortley’s results are obtained with additional 


information about the nondiagonal matrix elements of the 
electrostatic energy. The results of calculating the elec. 
trostatic energies by this method are presented for the 
configurations p°p, d°d and d°p. 





EISENBERG! has treated an almost closed 

shell by means of the quantum theory of 
wave fields. He considers those terms in the 
energy operator which are connected with atomic 
states arising from a configuration of equivalent 
electrons. He then changes from the variables N, 
(the number of electrons in the state 7) to the 
variables N,=1—WN, (the number of “holes” in 
the state 7). Corresponding to the transformed 
energy operator there is a Schrédinger equation 
in a configuration space whose coordinates are 
connected with the electrons missing from the 
shell. The results he obtains for atomic spectra 
and the anomalous Hall effect follow from this 
equation. 

More recently Shortley? has considered in 
atomic spectra the more general problem of an 
electronic configuration consisting of an almost 
complete shell plus other electrons not belonging 
in the shell. In calculating diagonal elements of 
the electrostatic energy matrix he showed how 
the sums which occur over the electrons in the 
shell may be replaced by sums over the electrons 
missing from the shell. His formulas may be in- 
terpreted by saying that so far as the diagonal 
electrostatic energy matrix elements go, the 
missing electrons interact with each other and 
with other electrons as though they were posi- 
tively charged. The formulas which he obtains 
for the matrix of the spin-orbit interaction may 
be interpreted similarly, for if the missing elec- 
trons act as though positively charged, the por- 
tion of the spin-orbit interaction connected with 


* National Research Fellow. 
1 W. Heisenberg, Ann. d. Physik 10, 888 (1931). 
2G. H. Shortley, Phys. Rev. 40, 185 (1932). 


them is reversed in sign. He also showed that in 
working out the J/; and Ws values belonging to 
the different atomic states (eigenvalues of L. and 
S., respectively) the m, and m, associated with 
each missing electron must be taken with a 
negative sign. 

It is the purpose of this note to point out that 
these results can be obtained rather nicely from 
the method Heisenberg used for an almost closed 
shell alone. Our procedure will differ slightly 
from that of Heisenberg in that we will not pass 
from the variables N, to a Schrédinger equation. 
Rather we will interpret operators in terms of 
these variables through their matrices, a method 
which naturally suggests itself as a result of the 
discussion on pages 628 and 629 of the preceding 
paper. 

As before, by assigning definite values to all 
the ./; (the number of electrons in a shell 7) we 
limit the distributions of electrons among the 
electronic states to correspond to all the atomic 
states arising from one electronic configuration. 
Our Hamiltonian is of the form 


H=I10'+]77, 
H?=3 poray*a g*a,a,(q| H*| rs), 


I'=)>",,a,+a,(p| II" r), 
(1) 


where we take //' to be the spin-orbit interaction 
and /7? the electrostatic repulsion of the electrons. 
In order that all the 7; have definite values, by 
Eq. (19) of the preceding paper we have 
M=>)H;, He=Dp(i)a,ta,(p|H"\r). (2) 


We must see what is required of H7? in order that 
all the 7; have definite values. Now 
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[H?, M;] = $>  peret|_@p*a q*QsQr, N. ]6i*(9 | H? | rs) . 


From the commutation relations among the a’s we find 


a,N;, — Nay = 5(p, tar, 
Hence 


a,*N:—N,a,* = —4(p, Aart. 


CH, Mi ]= —43>D porsdpta gta,a,(pq| H?| rs) (6;?+6;2—5 —5;*). (3) 


Therefore in order to have [ 7,24; ]=0 we retain in H? only terms for which r or s is in the subgroup i 
if p is in 7 or if g is in 7, or for which all four indices are in i. Hence we may write 


IP? =>) HitDdDsjAi;, 


y= 22 pors(t)Apta q*@.0,( pg | HT? | rs), 


His= Ved as(j)ap*ag*a,a,[ (pg| H?| rs) — (pg | H?|sr) ]. 


Now suppose in subgroups 1 and 2 we have 
almost closed shells so that it is advantageous to 
transform to the missing electrons in these sub- 
groups. We use the variables V,=1—N, if ris in 
subgroup 1 or 2; N,=N, if r is not in subgroup 1 
or 2. Then according to Heisenberg 


a;*=a, and d,=a,* if risin subgroup 1 or 2 
d,+ =a,* and d,=a, (5) 
if r is not in subgroup 1 or 2. 


The a@’s obey exactly the same commutation rela- 
tions as the a’s. 

Consider first the spin-orbit interaction. If 
1=1 or 2, H; transforms just as in Heisenberg’s 


(4) 
unchanged by the transformation. Thus 
iD = YH, 
H;= —Lopr(i)dp*a,*(p| H"|r)*, i=1, 2, (6) 


H;=Lpr(i)a,*G(p| |r), i>2. 
Here we have omitted a term which is the same 
for all our distributions and represents the mag- 
netic energy of the closed shells 1 and 2. Also we 
have used the fact (r| H"| p) =(p| H"|r)*. 
Consider next the electrostatic interaction. 
The terms ,,+He2+ He, transform just as in 
Heisenberg’s case. The terms H;; 1>2 and H;; 
i>j>2 are unchanged. There remain the terms 








case. If 7 refers to some other subgroup, H; is H;;i>2, 7=1,2 which on using Eq. (5) become 


Hij= ZX prli)>d os(j)Gp*G gas*G,[ (pq | I? | rs) — (pq| H?| sr) ] 
= — DL orlt)D os(J)Gp*d g*.4,[ (ps | H*| rq) — (ps | H*| qr) ]. 


A constant term representing the interaction of the electrons in the subgroup 7 with the closed shell j 
has been omitted. The first integral in the above expression can be brought in the convenient form 
(pq H?\rs) by using the complex conjugate wave functions for the “‘holes.”” However, this does not 
help with the second integral so we leave it as above. Dropping constant terms we have 


P?=> Hit Dollis, 





Hs=4% pros(i)Gp*d +G,4,(pq| H?|rs)*, i=1, 2, 
Hii =4D pors(i)G,tG gtG.a,(pq| H?|rs), i>2, i” 
Hey=> pr(2)¥ gs(1) Gp -*4.4,[ (pq | H?| rs)* — (pq| H?| sr)*], 
11; = —Xpr (i) gs(7) Gp (*G.4,[ (ps | H?| rq) —(ps|H?\gr)], i>2, 7=1, 2, 
Hij= Cor()Z oo(A)Gp*G G4, [(pg| H®| rs) —(pg| H*|sr)], i>j>2. 
We should also examine the angular momentum operators. We have 
L.=> Fete Li= Dd p(ia,ta,(p|l.|r) => p(t) N pm (8) 


and 
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S,=DiSsis S2i= DL p(t) Npmep. (9) 
Now 
La= —D>pr(i)d,+a,(p|1.|)*=—X (i) Npmy», i=1,2, Ia=X)p(i)Np»my, i>2. (10) | 
| 
Similarly | 
Su= —>,(i)N pm, i=1, ae S.:=>,(i)N mp 7>2. (11) 


With the set of eigenfunctions in the space of 
the N’s considered on p. 628 of the previous paper, 
we could construct the matrix of H' and of H? 
before transforming to the holes. According to the 
discussion of the preceding paper these matrices 
will be essentially the same as the matrix of H! 
and H? constructed with antisymmetric wave 
functions (Dirac-Slater determinants) in co- 
ordinate space. After the transformation the 
same set of eigenfunctions referred to the variable 
N, may be used to construct the matrix of H! 
and H?. Now the @’s bear the same relation to the 
N’s as the a’s to the N’s. Consequently if G has 
the same form in the @’s as G has in the a’s, its 
matrix may be obtained by comparison® with the 
matrix of G. For example 


G= LD verety q' s0,(pg | Go | rs), 
G=D porsdp*G gtG,4,(pq|Gi|rs). 


Then the matrix of G is obtained from the matrix 
of G by substituting in the latter the elements 
(pq|G.|rs) for (pg | Go| rs). 

On comparing Eq. (6) and Eq. (2) and remem- 
bering /7' is a real operator the matrix of the 
spin-orbit interaction is correctly obtained if we 
suppose the missing electrons in shells 1 and 2 act 
as though positively charged so that the part of 
the spin-orbit interaction associated with the 
holes is negative in sign. On comparing Eq. (7) 
with Eq. (4) and remembering H” is a real opera- 
tor the matrix of the electrostatic interaction is 
correctly obtained on the same assumption if in 
addition we replace the elements 


(pq| H*| rs) — (pq| H?| sr) 


which occur in the interaction terms between a 
hole and an electron, by 


(ps | H*|1rq) — (ps | H?| qr). 


8The comparison is between the matrix of G for a 
certain electronic configuration and the matrix of G for 
another electronic configuration obtained from the first by 
replacing the holes with electrons. 


For the diagonal elements (=r, g=s) the latter 
is equal to the former but for the nondiagonal 
elements it is otherwise. Finally on comparing 
Eq. (10) with Eq. (8) and Eq. (11) with Eq. (9) 
we see that the M, and Ms values of the atomic 
states are correctly obtained if the m; and m, as- 
sociated with the missing electrons are taken 
with a negative sign. Thus in the configuration 
d°p if m,;= —2, m,= —}3 for the d hole and m,=1, 
m,=% for the p electron, M,=3 and Ms=1. 

Appended are the electrostatic energies‘ for the 
configurations p°p, d°d and d°p as found by mak- 
ing the necessary changes in the calculations for 
the configurations p.p, d.d and d.p by Slater’s 
method. 


P’P 

‘D=—Fo— Fr, 1P= — Fot+5Fr, 

1D = — Fyo— F2+12G2, *S=—Fy—10F2, 

cP = = Fo+5 Fo, iS= = Fo- 10 F.+ 6G). 
d’d 


iG = — Fy)—4F,—F,, 

1G = — Fp—4 F2— F,+140G,, 
$F= — Fyt8F2+9F,, 

1F= —Fy+8Fo+9Fy, 

3D = — Fy + 3F2—36F,, 

1D = — Fo +3F2—36F,+28Ge, 
3P = — Fy—7F,+84Fy, 

1P = — Fy—7F 24 84F,, 

8S = — Fy—14F,—126F,, 

1§ = — Fy>—14F2—126Fy+10Go, 


4 The result for p*p has already been given by Shortley, 
reference 2. The notation here used for the F's and G’s is 
that of Condon and Shortley, Phys. Rev. 37, 1025 (1931). 

Added in proof: | am indebted to Professor Van Vleck 
for informing me that the electrostatic energies for d°p 
have been independently calculated by C. W. Ufford 
whose results agree with those given here. 
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d°p 
3F= — Fy—2Fo, 1D = — Fo+ 7Fr, 
1F= — Fy—2F2+90G;, *P=— Fy—7 Fo, 
sD = — Fot+7 Fo, 1P= — Fy —7F2.+20G,. 


Two things are at once apparent. The resonance 
separations of the singlets and triplets decrease 
with increasing L values; and starting with the 
next to the highest Z value, the resonance separa- 
tion for alternate L values is zero. These facts, 
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together with the larger value of the spin-orbit 
parameter for almost closed shells, explain why 
intermediate coupling sets in much sooner in the 
periodic table for closed shells than for the cor- 
responding configurations at the other side of 
the table. 

In conclusion the writer wishes to thank 
Professor Breit for many very helpful discussions 
and suggestions on the subject of these two 
papers. 
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Nuclear Magnetic Moments 
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This paper contains approximate formulas (Eqs. (5) and (6)) for the calculation of nuclear magnetic moments from 


observed hyperfine structure separations. 





HE hyperfine splitting of an energy level 
with total angular momentum Jh/27 is 
governed by the formula 


6v=ATIJ cos (IJ) 
=3A { F(F+1)—J(J+1)—1(1+1)}, 


in which J stands for the nuclear moment in 
units h/27, F is the fine structure quantum 
number, the resultant of J and J, and the cosine 
is supposed to be a “quantum” cosine. The 
proportionality factor A depends in a rather 
complicated way upon the quantum numbers of 
the energy state under consideration and is 
furthermore proportional to the nuclear g-value. 
In the hyperfine splitting A represents the 
proportionality factor obtained when one applies 
the interval rule. It is equal to the distance 
between two adjacent hyperfine levels divided by 
the largest of their F values and counted positive 
when the larger F value belongs to the higher 
energy. 

For a single electron one writes a instead of A 
and in a hydrogenic case the following formula 
holds: 

RaZ' g(1) 


a= cm7}, (1) 
n'(1+3)jG+1) 1838 





The symbols in this expression have the usual 
meaning, the nuclear g-value, g(J), is the ratio of 
the magnetic to the mechanical moment of the 
nucleus, the former expressed in “proton mag- 
netons” eh/4rMc, M being the mass of the 
proton. In analogy with the electron one expects 
for a single proton J= 3 and g(J)=2. 

For a single outer electron in a penetrating 
orbit one expects that an approximation similar 
to that derived by Landé for ordinary multiplet 


separations can be applied here. This gives!: 


ReZ;Z2 ~— (I) 
a= cm~!, (2) 
n§(I1+3)j(J+1) 1838 





The derivation of this formula assumes that the 
properties of the energy state can be considered 
as arising from two approximately hydrogenic 
regions of the atom, the outer one with an 
effective nuclear charge Z, and a principal 
quantum number »,, and the inner one with an 
average effective nuclear charge Z;. The value of 
Z, is 1 for a neutral atom, 2 for a state of a once 
ionized atom, etc., m is the Rydberg denomi- 
nator, which can be calculated if the absolute 
term value is known. The values of Z; are 
discussed below. 

For non-s electrons the ordinary spin doublets 
are given by the well-known formula 


Re®Z2Z2 
Av =————_ 
n,*l(}+1) 


cm", (3) 


Combining this with Eq. (2), assuming that 2; 
will have the same value in both formulas, one 
finds 
Av I(l+1) g(J) 
a= 


= . cm~!, (4) 
Zi(l+3)jU+1) 1838 





Earlier applications of these formulas gave 
unsatisfactory results, the nuclear magnetic 
moment calculated from the hyperfine structure 
of different levels of the same atom gave entirely 
different values. The discrepancy could be ex- 
plained only in part by the approximate nature 
of the formulas. Since then it has been shown by 


‘For this formula and the others compare Pauling and 
Goudsmit, Structure of Line Spectra, page 209, reference 1 
and page 60. 
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Breit? and by Racah* that large relativity 
corrections must be made to the above formulas. 
For non-s states a corrected equation (4) was 
used, which indeed seemed to remove the 
discrepancy between the values of the nuclear 
magnetic moment obtained from the two levels of 
the same doublet. For s states one did not use a 
corrected equation (2) but used numerically 
calculated approximate eigenfunctions because 
the factor in Eq. (2) is proportional to ¥(0)’, the 
square of the value of the eigenfunction at the 
origin. The calculations showed that there re- 
mained a large discrepancy between the values of 
the magnetic moment obtained from an s state 
and from other states. 

It is a remarkable fact that this discrepancy is 
much decreased if one uses Eq. (2) with relativity 
corrections instead of the approximate values of 
y(0)?. Already McLennan, McLay and Crawford‘ 
had shown that Eq. (2) is well in agreement with 
their observations on the hyperfine structure in 
the spectra of thallium. Furthermore Professor 
Breit® informed me that the value of ¥(0)? for the 
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6s state of Cs calculated by Nile by means of a 
more accurate method than the one used in 
previous applications comes out to be in very 
good agreement with Eq. (2). We are thus 
inclined to believe that the following equations 
(solved for g(J)) are suitable to give fair approxi- 
mations for the nuclear magnetic moment. 
For s electrons: 





3a 06Clom® =1838 
g(I)= , (5) 
8Ro? Z;Z27 (4, Z;) 
for non-s electrons: 
aZ; j(j+1)(1+4) ACI, 23) 
g(I) = 1838. (6) 








Av U(l+1) — x(j, Z:) 

In these expressions «(j, Z;) is the relativity 
correction by which Eq. (2) for the hyperfine 
structure has to be multiplied, A(/, Z;) is the 
same for the multiplet separation Eq. (3). These 
corrections are given by Eqs. (7) and (8). Table I 
contains these quantities for 7= 4 and 1} and for 
l= 1, for various Z; values.*: ® 





K(j, Zi) =479 +2) G4+1)/(40?—1)p, ep? = 4+3)?—(aZ;)’, (7) 
(I, Z;) =[2UI+1)/(aZ,)2]{ [(1+3)2— (aZ,)?}!—1—[P—(aZ,)*}}}. (8) 





In the application of the above formulas the 
following remarks are of practical importance. 
In many-electron spectra the same electron 
causes often part or all of the hyperfine splitting 
of several levels. Expressions are known’ to 
obtain the value of a for this electron from the 
various hyperfine structures. If these values turn 
out to be consistent with each other it may not 
be considered as a support for the formulas given 
in the present paper. Only when the values of 
g(Z) calculated from different electrons of the 
same atom agree can one take this as supporting 
the validity of the approximations involved. 

In a one-electron spectrum the choice of the 
quantities ”,, Z, and Z; is well defined. For 
s-electrons it will be correct to put Z;=Z, the 


? Breit, Phys. Rev. 38, 463 (1931). 

* Racah, Zeits. f. Physik 71, 431 (1931). 

‘McLennan, McLay and Crawford, Proc. Roy. Soc. 
A133, 652 (1931). 

5 Breit, Phys. Rev. 42, 348 (1932). 


actual nuclear charge, as the contribution to the 
interactions comes mainly from the region near 
the nucleus. For p-electrons Z; as obtained from 
the doublet formula (3) is throughout the 
periodic table about equal to (Z—4). It is not 

















TABLE I. 
K r 
Zi j=} j=} l=1 
10 1.01 1.00 1.00 
20 1.04 1.01 1.00 
30 1.09 1.02 1.01 
40 1.12 1.03 1.03 
50 1.30 1.05 1.05 
60 1.49 1.07 1.08 
70 1.78 1.10 1.12 
80 2.25 1.15 1.17 
8&5 2.61 1.17 1.20 
90 3.10 1.20 1.24 
92 3.36 1.21 1.27 








® Back-Landé, Zeemaneffekt und Multiplettstruktur, page 
96. 
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TABLE II.* Nuclear magnetic moments and g-values. 











Atomic Magnetic 
Z Element weight I moment g(J) Remarks 
3 Li 7 1} 3.29 2.19 
13 Al 27 3 2.1 4.2 B 
29 Cu 63 13 2.5 1.7 B 
65 13 2.5 1.7 
31 Ga 69 13 2.01 1.34 A g(71)/g(69) =1.27 
71 1} 2.55 1.70 
33 As 75 1} 0.9 0.6 C 
37 Rb 85 2} 3 a3 ¢ Rb87 perhaps J=2}, magnetic moment 
87 1} 2.7 1.8 3.0, g(Z) =1.2 
48 Cd 111 3 —0.67 —1.33 B It is not certain whether 111, 113 or both 
113 $ —0.67 — 1.33 cause the observed hyperfine structure 
49 In 115 4} 5.4 1.2 B 
51 Sb 121 24 2.7 1.1 B g(121)/g(123) =1.80 
123 33 2.1 0.6 
80 Hg 199 } 0.55 1.1 B g(199)/g(201) = —2.70 
201 13 —0.62 —0.41 
81 Tl 203 1 1.8 3.6 A g(205) about 1% or 2% larger than 
205 3 1.8 3.6 g(203) 
82 Pb 207 } 0.60 1.20 A 
83 Bi 209 4} 4.0 0.89 A 








* References: Li: reference 9. 
Al: Ritschl, Nature 31, 58 (1933). 
Cu: Ritschl, Zeits. f. Physik 79, 1 (1932). 
Ga: Campbell, Nature, in print. 


As: Tolansky, Proc. Roy. Soc. A137, 541 (1932). 


true that this is the same Z; as the one occurring 
in the hyperfine structure formula, the former is 
the average fourth power, the latter the average 
third power of the effective nuclear charge. One 
thus expects for the hyperfine structure Z; to be 
somewhat smaller than (Z—4), but how much 
smaller is quite uncertain. For heavy elements 
this uncertainty will be within the limit of 
accuracy of the method employed, but for lighter 
elements and also for larger /-values, d- and 
f-orbits, the error will become appreciable, 
causing the calculated g(J) to be too small. 

For a many-electron configuration one must 
first apply the formulas given in the above- 
mentioned papers®:7 to derive the a-values for 
the individual electrons. For non-s electrons 


7Goudsmit and Bacher, Phys. Rev. 34, 1501 (1929); 
Goudsmit, Phys. Rev. 37, 663 (1931). 





Rb: Kopfermann, Naturwiss. 21, 24 (1933). 

Cd: see Bacher and Goudsmit, Atomic Energy States. 

In: Jackson, Zeits. f. Physik 80, 59 (1933). 

Sb: Badami, Zeits. f. Physik 79, 206 and 224 (1932). 

Hg, Tl, Pb, Bi, see Bacher and Goudsmit, Atomic 
Energy States. 


similar formulas exist’ which give the Av separa- 
tions for the individual electrons, for p? and p* for 
instance, one has to use the total *P separation. In 
the case of p°, where the multiplet separations are 
zero in first order, all one can do is to take the *P 
separation of the next ion and decrease it 
slightly to allow for the screening by the third 
p-electron. 

When a single inner s-electron occurs in several 
configurations it is found that its a-value is 
practically the same for all of them. This means 
that there is hardly any screening effect caused 
by the other electrons. If there is a slight change 
the largest value of a should be the one where one 
expects the screening to be least and one may 
consider this a-value as the one the s-electron has 


8 Goudsmit, Phys. Rev. 31, 946 (1928); Humphreys 
and Goudsmit, Phys. Rev. 31, 960 (1928); Pauling and 
Goudsmit, reference 1, p. 157. 
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in that stage of ionization of the atom in which 
it is the only outer electron. One should use in 
Eq. (5) the Z, and , for that particular ion. For 
instance if one obtains a for the 5s electron from 
the hyperfine structure of the 5s 8s level in 
neutral cadmium one must use Eq. (5) for once 
ionized cadmium with this a-value. 

More complicated configurations cannot be 
used in general, unless one is able to obtain some 
detailed information from the multiplet structure 
about the coupling scheme of the various angular 
momentum vectors. The results obtained will not 
be very reliable. 

Table II contains the nuclear magnetic mo- 
ments and g-values for elements for which the 
hyperfine structure allows the application of the 
method described in this paper. The value for 
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lithium is not obtained by this method, but taken 
from the very accurate calculations made by 
Breit and Doermann.® All other results should be 
considered as preliminary only. In the table the 
results have been graded A, B, C to denote 
whether they have been obtained from a few 
independent hyperfine separations or not. Be- 
cause of the lack of sufficient data for several 
elements some personal judgment is needed when 
Eqs. (5) and (6) have to be applied to a more 
complicated electron configuration. I consider it 
therefore better not to discuss in detail the 
derivations of the results given in Table II, but 
rather wait for more extensive data to eliminate 
as much as possible the subjective factor. 


® Breit and Doermann, Phys. Rev. 36, 1732 (1930); 
Granath, Phys. Rev. 42, 44 (1932). 
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Paschen-Back Effect of Hyperfine Structure and Polarization of Resonance Radiation. 
Cadmium (6'P,—5'S,) 


N. P. HEYDENBURG, Department of Physics, State University of Iowa 
(Received February 3, 1933) 


Because of the Paschen-Back effect of the hyperfine 
structure of the odd isotopes, the polarization of cadmium 
\2288A resonance radiation excited by plane polarized 
incident light changes from 76.7 percent in zero field to 100 
percent in a strong field parallel to the electric vector of the 
exciting light. The experimental results agree very closely 
with computations based on Goudsmit’s modification of 
Darwin’s theory of the Paschen-Back effect and give for 
the separation of the two levels into which the 6'P, level of 


Cd (odd-isotopes) is split a value of 12.610"? cm™, 
From this result. and Schiiler’s measurements of the 
hyperfine separations of the triplet lines 6°Po, ;, »—6°S, the 
constants a and 6b giving the interactions of the p and s 
electrons, respectively, from Goudsmit’s equation 


61(2 —g) —2(g—1) L(L+1) 
(2L—1)(2L+3) 


are found to be a=8.4X 10-3 and 6=225.7 X 1073. 





A(J) =a(2—g)—a +6(g—1) 





HE polarization of the cadmium resonance 

lines AA2288A and 3261A has been dis- 
cussed by Ellett and Larrick! on the basis of 
Schiiler’s® * explanation for the hyperfine split- 
ting of the Cd energy levels. It will be recalled 
that Schiiler has assigned nuclear moments of 0 
and 3 (in units 4/27) to the even and odd isotopes 
of Cd, respectively. For the \2288A line, using 
nonpolarized exciting light, Ellett and Larrick 
found for the polarization in zero field 


P= (3+9y) /(11+9y) 100, 


where 7 is the ratio of the abundance of the even 
and odd isotopes which they found from experi- 
mental data to have the value 2.53. From obser- 
vations of Schrammen‘ and the narrowness of the 
Cd red line \6438.47A (7'D.—6'P)), it is known 
that the line \2288A is very narrow therefore the 
splitting of the upper level caused by the inter- 
action of the nucleus with the valence electrons 
must be small. Hence it is reasonable to suppose 
that a Paschen-Back effect should occur for these 
levels in rather small fields. 

Up to the present time the components of the 
\2288A line have not been resolved by spectro- 
scopic methods. However measurements of the 


1 Ellett and Larrick, Phys. Rev. 39, 294 (1932). 

? Schiiler and Briick, Zeits. f. Physik 56, 291 (1929). 
’Schiiler and Keyston, Zeits. f. Physik 67, 433 (1931). 
*Schrammen, Ann. d. Physik 83, 1161 (1927). 


polarization for different fields should lead to 
information on the separation of these compon- 
ents, for it will be shown later that the theoretical 
equation of the polarization as a function of the 
field strength depends on a parameter 8 which 
determines the separation of the hyperfine levels. 
In order to derive the expression for the polariza- 
tion as a function of the field strength it is neces- 
sary to obtain the various Zeeman transition 
probabilities as functions of the field. As the 
Zeeman pattern of the \2288A line for the even 
isotopes with J =0 is a normal pattern excitation 
with either perpendicular or parallel components 
alone will give complete polarization for all field 
strengths. Therefore these components need no 
further consideration. The Zeeman pattern for 
the odd isotopes with 7= 4 given in Fig. 1 is 
similar to that of the sodium D lines. Darwin® and 
also Heisenberg and Jordan® have published 
theoretical discussions of the Zeeman effect in 
intermediate fields for ordinary multiplets. The 
latter authors have worked out the explicit in- 
tensity formulas for the sodium D type doublets. 
However, these formulas must be altered some- 
what for our case as we are dealing with the inter- 
action of the nuclear spin and the angular mo- 
mentum of the valence electrons so that different 


* Darwin, Proc. Roy. Soc. A115, 1 (1927). 
6 Heisenberg and Jordan. Zeits. f. Physik 37, 263 (1926). 
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g values will be involved. Goudsmit’ has pointed 
out that we can use the formulas derived for 
ordinary multiplets by Darwin after replacing 
m, and m, by gim; and g;m;, respectively. After 
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making these changes and associating with a 
given Zeeman level the quantum numbers m; and 
m; appropriate to strong fields we find the intensi- 
ties to be given by the following formulas: 





|2’|?(1, 3, —23 0, 2, 2)=2{1+[6—w(g—g,) J/A_}, 
—})= 


|2’|*(1, —3, —3$;0, — ,, 





| x’ +iy’ 

|x’ +éy’|*(1, —9, —4; 0, — 
| x’ —ay’|2(1, 3, 330, — 

|x’ +ay’|?(1, —2, 2; 0, 2, 2) 


{ (8 —w(g—gi))? +867 }}; 





where A_= 


(1, —}, ite’ 


2{1+[8+(g—g:) ]/A4}, 
$)=2{1—[8—w(g—g;) ]/A_}, 
—}3)=2{1—[6+o(g—g;) ]/A,}. 
2) 
-)={1-b- w(g—gi) ]/A_}, 
t)= {1—[B+o(g—gi) J/A4}, 
3)=2, 
—$)= {1+[6—o(g—g;) J/A_}, 
= {1+[6+(g—g:) 1/44}, 


A,= {(8+w(g—g;))? +867}! 





and 8 is the separation constant of the hyperfine 
levels. The brackets on the left give the quantum 
numbers of the upper and lower levels involved 
in the transitions in the order (f, my, mi; f’, 
m,’, m;'). 

In calculating the polarization of the resonance 
radiation two additional factors must be taken 
into consideration, the intensity distribution in 
the source and the separation of the hyperfine 
levels in zero field. In order to obtain an approx- 
imate value of the separation of these levels the 
polarization will first be calculated for a source 
with uniform intensity distribution across the 
absorption region. It seems reasonable to sup- 
pose that this will give nearly correct values of 
the polarization as the source was operated at a 
fairly high temperature but very low Cd vapor 
pressure insuring a broad line free from self- 
reversal. The polarization is then obtained by 
the usual method of finding the population in 
resonance of the upper levels, then the intensity 
in resonance of each component, the polarization 
being given by P=(J,—J),)/(1,+J)), where J, 
and J), are the sums of the intensities of all the 
perpendicular and parallel components, _re- 


7 Pauling and Goudsmit, Structure of Line Spectra, p. 
219 


spectively. Completing these computations we 
find for the case where the incident light is plane 
polarized with its electric vector parallel to the 
field H 
16x*+32x°+27 
P\\= 100, (1) 
16x4+482"+ 63 





where x=gw/28=gehH/8xmcB. In plotting P), 
against H we obtain a family of curves, a curve for 
each value of the parameter £. 

The value of 8, giving that curve which most 
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Fic. 1. Zeeman pattern for the odd isotopes with J =} 
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Fic. 2. Intensity curves for the three components in the source. 


nearly coincides with the experimental curve, can 
then be used in calculating the polarization when 
the theoretical distribution of intensity in the 
source is taken into account. This intensity dis- 
tribution can be computed from the known for- 
mula for the Doppler broadening of spectral 
lines, 


I=I(c/vo) exp [ —a2c?(v/v9—1)?], (2) 
where a?=M/2RT, M being the molecular 
weight, R the molecular gas constant and ¢ the 
velocity of light. In Fig. 2 are plotted the inten- 
sity curves of the three components in the source 
taking its temperature to be 727°C. Curve (c) 
is the component due to the even isotopes. Curves 
(a) and (b) represent the two hyperfine com- 
ponents of the odd isotopes when the value of 8 is 
taken to be 4.2 X 10-* cm which was obtained by 
comparing the theoretical curves using Eq. (1) 
and the experimental curve of P); against / 
(Fig. 5). These add to give the resultant curve 
(d). As this curve has very nearly the same half 
value breadth as curve (c) it may be represented 
without appreciable error by Eq. (2) without 
changing the value of the temperature. Thus the 
energy radiated from the source at a given fre- 
quency range v to v+dy is given by 


d{[= I,(c/v4) exp [ —apcr(v/v1— 1)? |dv. (3) 


The energy absorbed in the resonance bulb by a 

component of frequency v2 for this same fre- 

quency range is 

dI’= {I2(c/ve) exp [ —ae2c?(v/ve—1)? ]} 
{Ii(c/v1) exp [ —avc?(v/y,—1)* ]}dv. (4) 


Integrating this value over all frequencies we 
obtain for the total energy absorbed 


B a;7a2c?(ve— 1)? . 
1'=— exp| - | (5) 


3...3 3...3 
V\V2 a;°V2" —a2°V71} 





where B is a constant depending upon the relative 
intensity of the absorbing line. The frequency 
shift of the Zeeman components with changing 
field strengths is readily obtained from the equa- 
tions for the energy of the Zeeman levels. The 
relative population in resonance of the upper 
levels can be calculated from Eq. (5), then the 
polarization is obtained as before. 

After comparing the curve obtained by plotting 
these new values of the polarization against the 
field strength with the experimental curve, it may 
be necessary to change the value of the separa- 
tion constant 8 to give better agreement. This 
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new value assigned to 8 must then be used to re- 
calculate the polarization. Hence it is seen that by 
a series of approximations we may approach very 
closely the theoretical polarization based on the 
theoretical intensity distribution in the source. 


EXPERIMENTAL 


The experimental set-up was the same as that 
described by Ellett and Larrick.' A monochro- 
mator was used toeliminate the \3261A resonance 
line. The vapor pressure of Cd in the source was 
maintained sufficiently low so that the red line of 
Cd was faint compared to the hydrogen //a line. 
Under these conditions no appreciable self-rever- 
sal of the A2288A line should have occurred. The 
resonance bulb was operated at temperatures 
between 100° and 105°C well below the tempera- 
ture at which depolarization due to secondary 
radiation occurs. A calibrated Helmholtz coil 
wound with copper tubing through which water 
flowed for cooling was used to obtain the magnetic 
field at the resonance bulb, its axis being parallel 
to the direction of the incident beam. 

Polarization measurements were made photo- 
graphically by the Cornu method with Wollaston 
prisms as in the work of Olson.’ In this method 
one of the Wollastons is rotated about their com- 
mon axis until two of the images on the photo- 
graphic plate are of equal intensity. The polari- 
zation is then given by the cosine of twice the 
angle of rotation measured from the position of 
the fixed prism. In the actual determination of 
this angle three or four photographs were taken 
at various angles within the region of the prob- 
able match position. A microphotometer curve of 
the cross section of each image was taken. The 
area under these curves as measured by a planim- 
eter was taken to be a measure of the intensity 
of the image. This is justified in the present case 
as the images compared were nearly of the same 
intensity and were exposed and developed to- 
gether. The function x =(C—D)/(C+D), where 
C and D are the intensities of the two images on 
a plate, was then plotted against the angle made 
by the Wollaston prism. Fig. 3 shows a typical 
curve obtained in this manner. The angle at 
which the line crosses the axis is taken as the 


* Olson, Phys. Rev. 32, 443 (1928). 
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angle where the two images would have exactly 
the same intensity. 

For the Cd line \2288A an exposure time of 

5 to 3 hours was necessary to obtain photo- 
graphs of sufficient intensity for measuring pur- 
poses hence after the general outline of the ex- 
perimental curve of polarization against field 
strength was obtained, the other points were 
found by estimating the angle at which the im- 
ages should be matched, then a photograph was 
taken at a small angle either side of this position 
and the actual match position determined from 
these. All of the points on the experimental curve 
with the exception of one were obtained before 
the form of the theoretical curve had been deter- 
mined. Several photographs were taken through- 
out the course of the experimental work with the 
heat removed from the resonance bulb and as no 
images could be detected on these films it was 
certain that stray light had been eliminated. 

It is necessary to correct the experimental 
values of the polarization for the depolarization 
caused by the slight convergence of the incident 
beam after passing through the lens system. The 
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depolarization due to the divergence of the beam 
which passed through the Wollaston prisms is 
negligible as the angle of maximum divergence of 
the rays was less than 3°. Gaviola® discussed this 
depolarization effect for the case of plane polar- 
ized incident light. A similar correction for the 
case of nonpolarized light is easily calculated. 
Fig. 4 illustrates how the incident intensity is 
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split up into components parallel and perpendicu- 
lar to the field. The intensity of a ray A may be 
represented by two mutually perpendicular 
equal components J; and J,’ one of which is 
chosen perpendicular to the field H. The other 
can then be split into components J». perpendicu- 
lar to the field and J; parallel to the field. Inte- 
grating over the aperture of the lens we obtain 
for the intensity of the perpendicular com- 
ponent 


over lens R r+2 L? 
I= f (11+J]2)ds= Tor f (——) rr 
0 \r+L? 


and for the parallel component 


over lens rR rdr 
n= f Ids= Ton f : 
0 #+L? 


where J) is the intensity of the radiation passing 
through a unit area of the lens. When the aper- 
ture is small compared to the focal length of the 
lens, Jo may be considered as constant over the 
lens aperture. The amount of depolarization 
resulting from the presence of the parallel com- 
ponent can be readily calculated for each field 
strength by determining the difference between 
the theoretical value of the polarization with and 
without this added parallel component. For zero 
field this correction becomes 





*Gaviola and P. Pringsheim, Zeits. f. Physik 34, 1 
(1925). 
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= 105(2—x) 
1149, 274-8x J’ 


where x= R?/L? while for very large fields it re- 


duces to 
Y=4x/(4+ x). 


The value of x for the experimental set-up used 
was 0.010. 


DISCUSSION 


Obviously we cannot compare the measured 
polarization values directly with those given by 
Eq. (1) as the measurements are for the case of 
excitation by perpendicular components and are 
the resultant polarizations from components of 
both even and odd isotopes. It can easily be 
proved that the polarization P, of the compon- 
ents from the odd isotopes alone in terms of the 
polarization Py of the components from both 
even and odd isotopes is given by the following 
formula 


Pi=(3y—3yPo—4Po)/(—y+yPo—4). (6) 
Then from the equation 
Py =2Pi/(1—P.) (7) 


relating the polarization P, for excitation with 
perpendicular components to the polarization P), 
for excitation with parallel components we ob- 
tain values which can be compared with those 
given by Eq. (1). 

In the first four columns of Table I are given 
respectively, the experimental values as meas- 
ured, the values corrected for the depolarization 
caused by the convergence of the incident beam, 
the values P, obtained from Eq. (6), and the 
values P,; obtained from Eq. (7). It will be 
noticed that the experimental value P» for zero 


TABLE I. Experimental values of the polarization for different 











magnetic fields. 

P Po Py P,, H 
observed corrected odd isotopes odd isotopes (gauss) 
76.3% 76.7% 27.3% 42.7% 0 
76.8 77.3 28.8 44.7 75 
79.3 79.9 36.0 53.0 144 
81.9 82.5 43.4 60.5 200 
85.5 86.2 54.4 70.5 255 
88.3 89.1 63.4 77.6 315 
91.3 92.1 73.0 83.7 375 
94.0 95.0 82.6 90.5 563 
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field differs slightly from that obtained pre- 
viously by Ellett and Larrick,' their value not 
having been corrected for the convergence of the 
incident beam. This gives 2.6 as the value of the 
constant y instead of 2.53. 

In Table II are recorded the theoretical values 
of the polarization for different field strengths. 


TABLE II. Theoretical values of the polarization as a function 
of field strength. 











, tl T H 
(calc. from 1) (calc.) (gauss) 
42.7% 42.7% 0 
48.6 48.9 100 
54.5 54.8 150 
62.0 62.4 200 
69.1 69.6 250 
75.0 75.6 300 
79.8 80.5 350 
83.3 84.0 400 
88.8 89.2 500 
91.8 92.1 600 








In the first column are the values computed from 
Eq. (1) by using 8=4.210- cm™, and in the 
second column are the values obtained on taking 
into account the theoretical intensity distribu- 
tion in the source. These differ very little from 
those in the first column. 

In Fig. 5 the theoretical values of the polariza- 
tion calculated from Eq. (1) are plotted against 
the field 7 and compared with the extrapolated 
experimental values givin in column 4 of Table I. 
Fig. 6 gives a comparison of the theoretical curve 
by using the values of the polarization given in 
column 2 of Table II with the experimental curve. 
As the general agreement of the curves in Fig. 6 
has not changed appreciably from that in Fig. 5 
it may be said that the consideration of the in- 
tensity distribution in the source has not effected 
a noticeable correction to the value of 8 previ- 
ously assigned. 

Fig. 7 gives a comparison of the experimental 
curve with two theoretical curves having values 
of 6 relatively 10 percent larger and smaller than 
that in Fig. 6. These curves lie outside the range 
of experimental error hence we may safely say 
that the value 8=4.2X10-* cm™ is correct to 
within 10 percent error. The separation of the 
two hyperfine levels of the odd isotopes is then 
38 = 12.6 10-* cm~. The accuracy of this deter- 
mination represents a resolving power of 3.5 X 10”. 
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As a resolving power of this magnitude is not 
available by use of present spectroscopic methods 
this is the only means of determining the separa- 
tion of these levels at present. 

In order to resolve spectroscopically the hyper- 
fine components of the line \6438.47A (7’D2 
— 6’P;) resulting from the separation of the 6’P, 
level a resolving power of 1.210° would be re- 
quired if the components had zero breadth. Since 
the line breadth cannot be entirely eliminated 
and the much stronger component due to the 
even isotopes is superimposed on the hyperfine 
pattern, it is not feasible at present to detect any 
variation from the intensity distribution given by 
a single line. 

Goudsmit"® has derived the following equation 
for the hyperfine splitting constant A, where A is 
just twice the constant 8 used above, 


A(J)=a(2—g) 
6r(2—g) —2(g—1)L(L+1) 
(2L —1)(2L+3) 


—da 





+b(g—1) (8) 


as usual [=[J(J+1)—L(L+1) —S(S+1)]/2, 
while a and } determine the absolute value of the 
splitting due to the and s electrons, respec- 
tively. The A’s for the triplet levels *P;, 2, *S; 
and the singlet level 'P; have been calculated 
using Eq. (8) giving 


A(®P;)=3a/2+6/2, =148.6X10%cm—, (9) 
A(®P2) = 3a/10+6/2, =120.0K10- cm, (10) 
A(®S,;) =), = 264.0 10-* cm~, (11) 


A('P:)=a, = 84X10 cm. (12) 


The numerical values given for the triplet levels 
were, obtained from Schiiler’s data®* on the 
separation of the hyperfine components of the 
triplet lines 6°Po, ‘. 2— 6° S; (AA4678, 4800, 5086A). 
The value of A('P,) is that found by the author. 
Schiiler’s results provide a unique solution as 
well as a check for the values of a and 0b. It will 
be noticed however that his results are not suffi- 
ciently accurate for a determination of the smaller 
constant a, for if the value of } from Eq. (11) 
is substituted into Eqs. (9) and (10) positive and 
negative values for a are obtained respectively. 


1 Goudsmit, Phys. Rev. 37, 663 (1931). 
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Fic. 5. Theoretical values of polarization from Eq. (1) and experimental values (column 4, 
Table I) plotted against magnetic field. Circle with dot, experimental curve; square with dot, 
theoretical curve. 8 =4.2 X10™ cm™. 
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Fic. 6. Plot of theoretical values of polarization (column 2, Table II) with the experimental values. 
Circle with dot, experimental curve; square with dot, theoretical curve. 
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Fic. 7. Comparison of experimental curve for polarization with two theoretical curves. Circle with dot, 
experimental curve; square with dot, A =9.25 X10™ cm™; triangle with dot, A =7.57 X10™* cm“. 


If we substitute the mean value of 6 from Schiil- from those obtained directly from Schiiler’s re- 


er’s results and the value from a of Eq. (12) into sults they are probably within the range of his 
Eqs. (9) and (10) we obtain 


A(®P1)= 135.0X10- cm—, 
A (®Ps)=124.9X10- cm—. 


experimental error. 

In conclusion I wish to express my thanks to 
Dr. Ellett for his many helpful suggestions 
throughout the progress of the work. 


As these values differ by less than 10 percent 
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The theory of the Stark effect in atomic spectra is 
discussed in general terms and it is shown that large 
effects arise when two terms which are related to each other 
in such a way as to satisfy the optical combining rules come 
close together in the spectrum. Fairly good values of the 


Stark displacements in complicated atoms may be obtained 
simply by using the hydrogenic values of the matrix 
components involved in the theoretical formulas. The 
ideas are illustrated by discussion of the existing data on 
the spectra of nickel, lithium, carbon spark, and argon. 





HERE exists in the literature of spectro- 

scopy a large amount of experimental data 
referring to the Stark effect in atomic spectra 
which has not been hitherto brought into any 
kind of relation to the modern theory of atomic 
structure. Of course the effect, and its theoretical 
interpretation, for atomic hydrogen is very well 
known and has played an important réle in de- 
velopment of the theory. Likewise the work of 
Foster! and of Dewey’ on the effect in helium has 
brought the theory into satisfactory relation with 
experiment there. Some work has also been done 
on the small quadratic effect on the resonance 
lines of sodium and potassium, but in addition to 
these contributions a glance for example at the 
Stark effect chapter in the H/andbuch der Experi- 
mental Physik (Volume 21) will show a consider- 
able amount of uncorrelated experimental data. 
In this paper the aim is to show how the theory 
may be used to throw light on some of this ma- 
terial and to point the way for further work in 
this field. 


§1. THEORY OF THE STARK EFFECT 


If Ho is the Hamiltonian of the unperturbed 
atom and if P is the electric moment of the atom 
and E the applied electric field, then the Hamil- 
tonian for the atom in the field is 


H=H,)—E-P. (1) 


The alteration of energy levels and proper states 
is therefore governed, in accordance with usual 


1 Foster, Proc. Roy. Soc. All4, 47 (1927) and A117, 137 
(1927). 
? Dewey, Phys. Rev. 28, 1108 (1926) and 30, 770 (1927). 


theory, by the matrix components of P. More- 
over the only component of P that is effective is 
that along E so we may choose axes so this is the 
z-component. Now P, is the same quantity which 
governs that part of the dipole radiation of atoms 
which is polarized with its electric vector parallel 
to the z-axis. Therefore the selection rules and 
other calculations of its matrix components which 
have been made for the theory of line intensities 
are applicable here. Writing (A|P|B) for the 
matrix component of P, which connects two un- 
perturbed energy states A and B we know that 
(A|P}B) vanishes unless: 

(a) A and B are of opposite parity (Laporte 
rule), (b) Je=Ja or Ja+1 where Jy and J, are 
the resultant angular momentum quantum num- 
bers of B and A respectively, 

(c) Mya=Myp where My, and My,» are the 
z-components of resultant angular momentum of 
the atom in the two states measured with h/2r 
as unit. 

These three properties are rigorous. In addi- 
tion, insofar as it is accurate to assign electronic 
configuration labels to the terms we may say: 

(d) Configuration A may differ from configura- 
tion B at most in regard to one electronic n, | 
value. 

Likewise in case of Russell-Saunders coupling 
in which resultant spin angular momentum S and 
resultant orbital angular momentum L are quite 
accurately diagonal when [7p is diagonal we may 
say that (A|P!|B) vanishes unless 

(e) Sa ™ Sp. 

(f) Le = ) am La + 1. 

In addition for Russell-Saunders coupling the 
actual dependence on L, S, J and M, is the same 
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as that which in the intensity theory leads to the 
well-known Kronig-Russell-Hénl formulas. 

So far as orders of magnitudes are concerned 
we may expect that the nonvanishing components 
will be of the same order as if calculated with 
hydrogen-like eigenfunctions for the individual 
electrons of the atom. For this reason we note 
that Schriédinger’s calculations for hydrogen give 


(n, 1—1, m| EP |n, l, m) 


3 
=— eEndy | 
2 





(n?—1?)(12—m?)}! (2) 
4/?—1 | 


for the matrix component connecting two states 
of equal » value and of / value differing by unity. 
Here dp is the first Bohr radius. Expressing en- 
ergy in cm™ and taking 100 kv/cm as the unit of 
E the coefficient of the irrational expression in 
(2) becomes 


6.43n cm per 100 kv/cm. 


The maximum value assumed by the radical for 
fixed m and varying / and m is about /3' so far 
n not greater than 10 the matrix component is 
under 200 cm in value in all cases for fields up 
to 100 kv/cm. This field is about the largest which 
has been used in ordinary Stark effect investiga- 
tions. 

Suppose now we put such a perturbing term 
into the Hamiltonian of an atom and consider 
what may be said in general about the result. 
The perturbed states will become linear combina- 
tions of such of the unperturbed states as are 
joined by nonvanishing matrix components of 
the perturbation. Therefore quantities which 
were quantum numbers in the unperturbed prob- 
lem lose their significance after the perturbation. 
Thus it is no longer possible to speak of a state as 
being definitely odd or even, also the J value loses 
its precision, but /, retains its status as an ac- 
curate quantum number. 

The extent to which this tendency is realized 
the perturbation theory shows to be related to 
the nearness of the unperturbed terms. If W4 and 
Wz are the unperturbed energies of states A and 
Band (A|EP|B) is the perturbation component 
connecting them then the perturbed states will be 
appreciable linear combinations of A and B if 
(A | EP| B) is comparable with or large compared 


to (W4—W,). If (A| EP| B)<(Ws— Ws) then 
the alteration in the perturbed states is relatively 
small as is also the change in energy. In this case 
so far as the change in energy is concerned due 
simply to the matrix component connecting A 
and B it is such as to displace the upper state 
upward and the lower state downward each by 
the same amount, namely. 


| (A | EP| B)|?/|Wa—Wal. (3) 


The states A and B thus seem to “repel” each 
other. 

As to the amount of this “repulsion” of two 
terms it is generally rather small. For weaker 
fields a common value of (A|EP|B) would be 
15 cm whereas the terms (A and B) may be 
separated by 1000 cm~ so the perturbation of 
each would amount to 0.22 cm. The small 
quadratic Stark effect of any state A is the sum 
of the actions of this type of all other states B of 
the atom which have matrix components connect- 
ing it with A. Most of the theoretical work on the 
Stark effect of alkalis, like that of Unséld* and of 
Kirkwood,‘ is concerned with attempts to sum 
up, for particular states A, like the normal state, 
or the first resonance state, the whole effect of all 
other states which perturb it. 

But such effects are never comparable with the 
relatively large effect in atomic hydrogen. 
Nevertheless such large Stark effects are ob- 
served in other atomic spectra. Large Stark 
shifts may occur in other atoms if two terms con- 
nected by a matrix component (A|EP|B) are 
closer together in the unperturbed energy scheme 
than the energy value (A | £P| B). In particular 
if A and B have zero interval in the unperturbed 
scheme then they will give rise to two levels in 
the perturbed scheme, one of which is moved up 
by (A|EP|B) and the other down by the same 
amount. This makes it clear why large Stark 
effect displacements are relatively rare in spectra 
other than hydrogen. The matrix components 
(A|EP|B) are rather severely restricted by the 
selection rules. Even when nonvanishing they are 
quite small so that it is something of an excep- 
tional case when two states A and B which are 
connected by the perturbation matrix are close 


3 Unséld, Ann. d. Physik 82, 355 (1927). 
‘ Kirkwood, Phys. Zeits, 33, 521 (1932). 
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together relative to the magnitude of (A | EP|B). 

The condition for large Stark effect just set 
out is more general than the one given by Bohr® 
before the new quantum mechanics which is 
usually used by experimental workers in qualita- 
tive discussions of their data. Bohr pointed out 
that we may expect large effects for a spectral 
term that is nearly hydrogen-like, i.e., one which 
if equated to Rh/n*® leads to a nearly integral 
value of n*. This was regarded as indicating that 
the electronic orbit was effectively in a Coulomb 
field and hence the hydrogen-like behavior was 
to be expected. But it did not give a means of 
estimating the magnitude of the effect. The ac- 
count of the theory we have given shows clearly 
the justification of Bohr’s rule. If a particular 
term is almost hydrogen-like then we may be 
fairly sure that terms arising from an electron 
configuration in which one of the electronic / 
values is increased by unity will also be hydrogen- 
like. If both are nearly hydrogen-like terms they 
will be near each other and so produce a large 
effect. The essential thing is nearness of W, and 
Wz for nonvanishing (A | EP| B) and Bohr’s rule 
gives us a convenient way of noticing what is in 
point of fact a large class of the cases in which 
large Stark effects are observed. The rule is suffi- 
cient though not necessary so we may expect to 
find cases in which the effect is large even though 
the perturbed states are not hydrogen-like in 
their energy values. 

The rule was enunciated before the develop- 
ment of Pauli’s principle with its sound basis for 
getting true ” values. In those days if ”* was 
nearly an integer it was thought that the nearest 
integer was the true value of ». Now we know 
that n*—n may become as great as several units 
so that a term may be apparently hydrogen-like 
through having (n*—n) be almost exactly equal 
to an integer instead of zero. This is in fact the 
case in silver, where there is a good-sized Stark 
effect that comes about from the nearness of the 
6d*D to the 4f?F, both being very close to the 
hydrogenic value for »=4. 

In what follows attention will be paid to the 
application of the perturbation theory to the dis- 
cussion of the existing data for the spectra of 


®’ Bohr, Proc. Phys. Soc. London 35, 275 (1923) and 
Ann. d. Physik 71, 228 (1923). 
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several elements. The general standpoint will be 
that the effect has to be calculated by finding the 
allowed values of /J,;—E-P which can be done in 
view of the smallness of the matrix components 
of E-P by solving the finite secular equation 
which represents the interaction of just those 
interacting terms which are within a few hundred 
wave numbers of each other on the energy scale. 


§2. NICKEL 

To illustrate how the principles of the preced- 
ing section may be used to throw some light on 
the experimental data, even in a case that is too 
complicated for detailed discussion at present, let 
us consider the Stark effect in nickel. The only 
observations apparently are those of Takamine,® 
which were made a decade before the analysis of 
the spectrum by Russell.’ Of the fifty lines for 
which Takamine records Stark displacements, 
twenty-nine may be found in Russell's list of 
classified lines. These are all transitions from 
levels of the high even configurations d®s5s, d°4d, 
d*6s, d°5d and d’s4d to the intermediate odd con- 
figurations d°4p and d*s4p. The final states have 
term values between 25,000 and 33,000 cm 
above the normal state. As there are no even con- 
figurations in this region we conclude that the 
Stark displacement of these terms is very small. 

The initial states involved in the identified 
lines fall into two groups, one between 49,000 and 
51,000 cm~', and the other between 54,000 and 
57,000 cm. Taking the ionization limit to be 
61,579 cm”! we find that the hydrogenic terms 
are 49,400 cm~'! for »=3 and 54,730 cm~ for 
n=4. For the lower group the configurations are 
(3d)°4d and (3d)*4s5s, so these are not u=3 
states for valence electrons. But these lower even 
configurations occupy the same region in the 
energy diagram as does the odd d°5p which is 
known from Russell's analysis. So it is more in 
accord with our general principles to regard the 
Stark effect as an interaction of d°4d and d's5s 
with d°5p which is induced by the applied field 
and to regard the nearness to a hydrogenic value 
as accidental. Similarly the upper group of even 
terms is in the correct position to interact with 
d°6p but as this configuration has not been identi- 

6‘ Takamine, Astrophys. J. 50, 1 (1919). 

7 Russell, Phys. Rev. 34, 821 (1929). 
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fied in the spectrum as yet we must confine our 
attention to the lower group of initial states. 

This further restriction reduces the number of 
identified lines for which Takamine gives Stark 
displacements to twelve, listed in Table I. The 
first column gives Takamine’s wave-length, the 
second that of the identified line in Russell's 
paper which I take to be the same line (there is a 
systematic difference of about 0.2A). The third 
column gives the identification in terms of Rus- 
sell’s multiplet analysis, the fourth and fifth are 
the shifts in wave numbers of the parallel and 
perpendicular polarized components as observed 
by Takamine for the values of the field strength 
given in the sixth column of the table. The line 
5142.77 is given two alternative identifications by 
Russell both of which are listed in the table. 

It will be noticed that e*Fs occurs as initial 
state for the first, second, fifth and eleventh lines 
in the table. The nearness of the equality of their 
Stark displacements, when these are reduced to 
the same field strength assuming a variation with 
E?, is an indication of the accuracy to which the 
whole effect may be attributed to the upper state. 
In the twelve lines there are eight different initial 


TABLE I. Stark effect in nickel I. 











(Taka- Av (cm~) Field 
mine) (Russell) Identification || 1 (kv/cm) 
4410.66 4410.50 25D,°—-eF, +6.1 +4.1 39 
4937.45 4937.33 2F,°-—eF, +50 +5.5 38.5 
5018.48 5018.30 2°D,°—e’F, -—2.8 0.0 38.5 
5082.55 5082.38 2°P,°—eP,; —-1.7 -0.97 38.5 


5084.20 5084.07 2°D;°—e'Fy +50 +54 38.5 


5142.91 5142.77 \aRe_PD. +0.34 40.26 21.8 


5146.61 5146.48 2°D.°-eF, +044 +0.26 21.8 
5155.90 5155.76 2'D.°-e'F; +30 +30 21.8 
5176.72 5176.56 2'D.°—f'D, +0.56 +0.34 21.8 
5184.78 5184.59 2'D,°—e'P, —1.1 —1.1 21.8 
5462.69 5462.48 3c'F;°-e' Fy +2.9 41.9 21.8 

21.8 


5588.09 5587.85 a*P.,—y'D;° +048 +0.35 








states represented. All the terms show a displace- 
ment upward in energy except e*P; and e*F». 
According to the ideas of §1 we therefore expect 
that the nearest odd combining term to each of 
these terms will be below the perturbed term, 
except for these two for which the nearest com- 
bining term is above. This is in fact the case as 
Table II shows. In Table II are listed in the suc- 
cessive columns: the name of the term, its con- 
figuration, the name and configuration of the 
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TABLE IT. 

Perturbed Nearest Interval 
term Configuration perturbing term (cm~) 
eF, d°4d 3D3°d°Sp $3 
e*F d°4d 3D\°d°5p +16.7 
e'F; d°4d 3F°d*5p —42.5 
e'P, d°4d 3D.%d°5p +12.0 
eF, d®s5s 6°*J =3 —53.3 
f?Dz d°4d 3D °d*Sp —0.5 
; 3D d°Sp —86.5 
fds od 3D °d°5p +56.1 
f'D, d°4d 1D..d°5p — 65.0 


y'D;° d*s4p 8D;d°5s —15.11 








nearest term which could perturb it, and the in- 
terval between the perturbed term and this 
nearest term in cm counted negative if the per- 
turbing term is below the perturbed term and 
positive if it is above. Of course the designations 
perturbed term and perturbing term are simply 
relative to the particular lines under discussion; 
actually the perturbation is mutual and (if it 
existed) Stark effect data on lines involving the 
terms called “perturbing”’ here should show that 
they are perturbed by amounts equal and oppo- 
site to their perturbing action on the terms here 
called ‘‘perturbed.” 

The table shows that we should expect f*D2 to 
show a strong upward perturbation whereas f*D; 
would probably show a very small perturbation 
since the nearest terms are relatively far from it 
and there are two, one above and one below of 
roughly the same distance whose effects tend to 
cancel. Since the Stark effect for 5142.77 is 
actually quite small this tells us unambiguously 
that this line is to be identified with the second 
of the two alternatives offered by Russell's 
analysis. It is believed that this is the first time 
that a doubtful point in an analysis of a complex 
spectrum has been settled by reference to the 
expected Stark effect of the lines. 

In the absence of more complete experimental 
data it is not thought worthwhile to attempt more 
precise discussion of the relations involved in the 
theory. It is felt, however, that even these rough 
agreements suffice to show that the Takamine 
data are in accord with the general outlines of the 
theory of the Stark effect. 


§3. LitHIuM 


The experimental data for lithium are sum- 
marized in the Handbuch der Experimental 
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Physik (Volume 21) article by Stark, page 479. 
All of it refers to combinations in which 2s or 2p 
are the final states. Either of these should have 
extremely small effects owing to their great sep- 
aration from the nearest terms of opposite parity. 
Therefore we expect the line shifts to be essen- 
tially those of the corresponding upper states. 
Evidence for this is found in comparing the 
recorded shifts of 2s-5p and 2-5), assuming 
these to be quadratic effects. The former is 
—5.5 cm at 0.26 unit (1 unit=10° volt/cm) 
while the latter is —45.5 cm™ at 0.8 unit. Writ- 
ing Av=kE* corresponding values of k are 81. and 
72. respectively. 

In this spectrum the doublet separations are 
negligibly small. Also the excited states are fairly 
hydrogen-like. Let us first consider the effect of 
the field on the states with m = 4. The unperturbed 
energy values relative to the 4d term are: 


4s —1611.7 cm" 4d 0.0 
4p — 154.7 4f +7.4 


The matrix components of electric interaction 
we may expect to be close to the values given by 
(2) for hydrogen. In any case the relative magni- 
tudes will be given quite accurately by (2) and 
this is all that matters for a description of the 
pattern. We may expect a fairly strong interac- 
tion between 4d and 4f, a smaller effect on 4p and 
a quite small effect on 4s. The secular equation 
for m= 3 is linear and tells us that this substate of 
the 4f term is not affected by the field. Owing to 
the selection rule on m this does not show up in 
combinations with 2p. For m=2 it is a quadratic 
equation connecting these substates of 4d and 4f. 
Using the matrix components for (2) and the 
empirical 4d—4f separation we find the roots to be 


=3.74[3.72+25.8°E?]} 


when measured from the 4d level. 

Likewise for m=1 the secular equation is a 
cubic connecting corresponding substates of the 
4p, 4d and 4f terms: 


—154.7-rX —39.8E 0 
— 39.8E —X —32.5E |=0. 
0 —32.5E 7.4—X 


For m=0 the secular equation is a quartic but 
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here the 4s term is so far away compared to the 
size of the perturbation that its effect may be 
neglected. The matrix component connecting 4s 
and 4p has the value 56.9 for unit field strength 
so the perturbation due to interaction of these 
terms is (56.9)?/1457.0 = 2.2 cm which is negli- 
gible compared to the main effect. Even the 4p—4d 
interaction is quite small so that all of the secular 
equations may with good accuracy be taken 
simply as quadratics connecting the 4d and 4f 
levels. If we do this we obtain as the levels 


Ae = 3.7+[ (3.7)? + 25.8? E? ]}} 
Ai =3.74[ (3.7)? +32.5°E? }! 
Ao =3.74[ (3.7)? +34.5°E? ]}. 


In Fig. 1 these levels are plotted against E. The 
experimental points are from the date of Snyder® 
on 4602 which is the combination with 2p. The 
agreement is quite good. This calculation makes 
no reference to the doublet character of the 
spectrum and so is not in accord with Snyder’s 
suggestion concerning the line’s behavior.® 
Similar remarks may be made concerning the 
group of levels for 2=5. In combination with 2p 
only the substates with m=0, 1 and 2 are effec- 








20F ‘4 SOF 


{0r 








0.1.2.5 





i 




















2 


4 


Fic. 1. 





2 
6 





4 


Fic. 2. 





Fic. 1. Stark displacements due to interaction of 4d and 
4f terms in lithium. Ordinates, displacement in cm; 
abscissas, field strength in 10* volt/cm. The points are 
observed displacements of the lines 4d—2p and 4f —2p. 

Fic. 2. Stark displacements due to interaction of 5d, 5f 
and 5g terms in lithium. Ordinates, displacement in 
cm~; abscissas, field strength in 10‘ volt/cm. The points 
are observed displacements of the combinations with the 
2p term. 


8 Snyder, Phys. Rev. 33, 354 (1929). 
® Reference 8, bottom of page 359. 
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tive. The energies relative to 5d are: 


5s —798.2 Sf +78 
Sp — 84.0 Sg. 7.8<?<10 
5d 0.0 


where probable limits for the unknown 5g are 
merely suggested. In the field the perturbed 
eigenstates become linear combinations of the 
almost coincident 5d, 5f, 5g terms. The 5 is not 
so far away but that its displacement is appre- 
ciable. With the second-order perturbation 
formula this term must be pushed up by 5s and 
down by 5d giving as its displacement in unit 
field —57 cm~ which is in fair agreement with 
the values —81 and —72 obtained experimen- 
tally. The displacement of the line 2p-5s corre- 
sponds to a downward displacement in unit field 
of the 5s term of 28 cm™. The theoretical value is 
12. cm™. 

The group 2p— (5d+5f+5g) may be discussed 
easily by omitting the 5/ interaction and neglect- 
ing the small and empirically uncertain intervals 
between these initial state terms. Doing this the 
energy levels are given by 


Wwe=+64.2E, A=ATTSE, A=+81.4E. 


In Fig. 2 the data of Snyder are plotted together 
with the straight lines given by the theory. 

Similarly the line 2p—6s shows a shift corre- 
sponding to a displacement of 6s of —32 cm™ 
in unit field whereas the theoretical push from 6p 
is —45 cm. The displacement of the 6p term 
in unit field as inferred from the 2p-6) line is 
—163 cm. The theoretical value is the differ- 
ence between the upward push from 6s and the 
downward push of 6d which comes out to be 
—200 cm-'. Since the experimental data are 
rough and the theoretical values depend on the 
square of matrix components which are taken 
from hydrogen it is felt that the agreement here 
presented is good enough to indicate that the 
ideas of §1 are adequate for a discussion of the 
main effect. 


$4. CARBON SPARK 


The data are some observations made by 
Ishida and Fukushima" and are believed to be 


© Ishida and Fukushima, Sci. Papers Inst. Phys. Chem. 
Research 14, 123 (1930). 
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the only data on a nonhydrogenic ion spectrum. 
The displacements, as is to be expected, corre- 
spond to smaller values of the matrix components 
than in arc spectra. The matrix components are 
in fact roughly half the hydrogenic values which 
is the factor of reduction to be expected for orbits 
in a Coulomb field with Z = 2. 

Their Fig. 2a gives the results of observation 
on \2747.31 and 2746.50 which is the transition, 
S3p?P—s*4d°D, in CII. The shift is to longer 
wave-lengths and corresponds to k= —0.004 in 
the equation Av=kE? where Ay is in cm=' and E 
is in 10‘ volt/cm. Likewise their Fig. 2b gives the 
data on 4267.27 and 4267.02 which is s°3d 
—s*4f. Here the k is about +0.002 in the same 
units. Attributing these shifts to the upper 
states we see that they are due to the interaction 
of the s*4d term and the s*4f term. The 4d term is 
855 cm below the 4f term and so the 4d term is 
pushed down in energy and the 4f term is pushed 
up, which is in accord with the facts. 

An interesting feature of their Fig. 2b is the 
fact that there is an unshifted component in the 
. polarization which is absent in the || polariza- 
tion. This is due to the fact that the m=3 sub- 
states in *F are unperturbed since there is no 
m=3 substate in *D to perturb them. In com- 
bining with 7D to form a radiative transition the 
m =3 substate of ?F must jump into an m= 2 sub- 
state of 2D which accounts for the presence of 
the unshifted component in the | but not in the 
|| polarization. 

The hydrogenic matrix component for m=0 
connecting the 4d and 4f states for a field of 10* 
volt/cm is 3.44 cm from Eq. (2). The value of 
the matrix component corresponding to empirical 
k of 0.002 and Av of 855 is 1.37 cm™ or roughly 
half of the hydrogenic value. This is largely due 
to the fact that we are dealing with an ion rather 
than a neutral atom, also to the fact that the 
pictures are blends of components for m=1 and 
m=2 for which the matrix components are 
smaller. . 

Their data on 2992.63 furnish an excellent 
example of the interaction of terms as sketched 
in $1. This is the transition s°3d—s*5f. The ex- 
perimental value of k is —0.018. At first sight 
this appears to be an exception to our general 
rules as the field makes it go down in energy and 
yet the 5d term, below it, should push it up. 
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This is due to the presence of a 5g term. This 
term is not known experimentally but we will not 
be far wrong if we place it at the exact hydrogenic 
value for n=5. The 5f term is 144 cm™ lower 
than this and the 5d is 462 cm™ lower than 5f. 
Assuming hydrogen-like matrix components (di- 
vided by 2 because we have an ion) and these 
empirical term separations we calculate for k the 
value 


k = —(2.43)2/144+ (3.45)2/462 = —0.015. 


Here the first term is the downward push due to 
5g and the second the upward push due to 5d; the 
final value is in good agreement with the experi- 
mental, —0.018. 

These are all the data on this spectrum which 
are available. Theory agrees with experiment as 
accurately as could be expected and, in particu- 
lar, explains the difference in the behavior of the 
4f and Sf terms. 


§5. ARGON 


The effect in argon has been investigated re- 
cently by Ryde." He measured the displacements 
of a large number of lines and has given in his 
Table VII the displacement of the energy levels 
produced by a field of 100 kv/cm. The configura- 
tion p°5d is given with some completeness and his 
Fig. 3 shows how the displacement of a term is 
greater the smaller its difference from the hydro- 
genic value for m=5. In argon the states may be 
separated into two sets, one built on the state 
p’*P, of the ion, the other on p**P, of the ion. 
Combinations between the two sets are weak. 

Looking at the data from the viewpoint of this 
paper, we see that the p*5d(?P;) terms are per- 
turbed downward because of a repulsion from 
the p°5f(?P;) terms which lie in a close group at 
just the hydrogenic value reckoned down from 
*P, of A II. Hence distance from hydrogenic 
value is here synonymous with distance from 
nearest perturbing term. At present we do not 
have very good knowledge of the coupling rela- 


l Ryde, Zeits. f. Physik 77, 516 (1932). 
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tions in argon so it is hard to give a detailed dis- 
cussion of the interaction of the 5d with the 5f 
terms. Neglecting the separation between the 5f 
terms, the perturbation of any particular 5d term 
due to the 5f terms is 


| (Sda| E+ P| 5fs) |? 
B  (Wsa*— W5/) 








where the sum extends over all the 5f terms in the 
group. If the numerator is about the same for all 
the d terms the Av values should lie on a hyper- 
bola whose axis is the location of the 5f group, i.e., 
about 4400 cm. Using the hydrogenic value with 
m =0 of the 5d—5f interaction from Eq. (2) for the 
numerator we have 5770. In Fig. 3 is plotted the 
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Fic. 3. Stark displacements of the several p°5d(?P3,2) 
terms in argon, showing relation of the displacement to 
nearness to the perturbing p*5f and p*°6p groups of terms. 
Ordinates, displacement in energy in field of 10° volt /cm; 
abscissas, term values of the terms in cm~'. The hyperbola 
is drawn, using hydrogenic values of the matrix component 
connecting 5d and 5f. 


curve (5770)/(Wsa—Ws,;) together with points 
showing Ryde’s experimental values for the cor- 
responding term displacements in a field of 100 
kv/cm. It is seen that the observed values agree 
with this curve in order of magnitude but drop to 
zero more rapidly than does the hyperbola repre- 
senting interaction with 5f. This might be due to 
a variation of the numerator with the 5d, but 
another cause which is certainly acting is the up- 
ward push on the 5d terms from the 6p group 
which is not so very far below the 5d terms. 
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On the Magnetization of Ferromagnetic Crystals 


FRANCIS BITTER, Westinghouse Research Laboratories 
(Received December 15, 1932) 


In a previous article a function Eg is discussed. This 
function gives the energy of a crystal as a function of the 
direction of magnetization. The present paper contains 
photographs of plaster models representing Eg for an 
undistorted crystal of nickel, for an undistorted crystal of 
iron in zero applied field, and in a field of 100 oersteds 
parallel to the [100], [110], and [111] axes respectively, 


and for an iron crystal distorted by compression and 
extension along the above mentioned axes. Assuming that a 
crystal is magnetized to saturation in the direction in 
which Ee has a minimum, magnetization curves for 
crystals distorted as mentioned above are calculated and 
illustrated. 





N a previous article! a procedure was outlined 

for describing some properties of strained 
ferromagnetic crystals. It is proposed to discuss 
the consequences of this procedure in greater 
detail. In the following pages certain aspects of 
magnetization are taken up. 


PROCEDURE 


Let i, 7, k be subscripts denoting measurement 
along three axes parallel to the tetragonal axes of 
a cubic ferromagnetic crystal whose spontaneous 
magnetization is J,. The orientation of J, is 
given by the direction cosines a@;, aj, a,. The crys- 
tal may be distorted and the distortion described 
by the components of a strain tensor, A;;. The 
energy for an arbitrary direction of magnetiza- 
tion is assumed to be:— 


Ey = cy'aa;?+K,>0A na?+Ke>-'A {jQjQ; 
—I,,>-a;}1;+const. (1) 
From this expression it is possible to calculate the 


magnetostriction in an arbitrary direction 8), 


B;, Bx. 
81/l= xo+x10;78?+ x°D'a,0/8B;. (2) 


The constants involved in Eq. (2) are related to 
those in Eq. (1). 


x1i= —Ki/(¢u—¢i);) x2= —Ke/2cs. (3) 


The quantities ¢);, C12, C4, are the elastic constants 
of the cubic crystal under discussion. It is further 
assumed that the actual direction of magnetiza- 





'F. Bitter, Phys. Rev. 42, 697 (1932). 


tion of the crystal for any given set of external 
conditions is that direction in which Ey has a 
minimum. If Ey has minima in two or more direc- 
tions ambiguities may arise. 


EVALUATION OF CONSTANTS 


For nickel? the constant ¢ is approximately 
— 2.5104 ergs/cc and J,,=500. For iron all the 
constants may be evaluated with sufficient ac- 
curacy to give some idea of the predictions of 
Eq. (1). From observations by Honda and 
Masiyama® we can evaluate xo, x1, and x2, as has 
been done in a previous article.' Certain con- 
stants, S)1, etc., related to the above elastic con- 
stants, €);, etc., by the following equations: 


C11 =Sirt$12/ (S11 +2512) (S11 — Si), 
Ci2= —$12/(Sir +2512) (Sin —Si2), 


C44 = 1/544, 


were measured for iron by Goens and Schmid.‘ 
From their measurements we obtain ¢,;=0.237 
X10"; ¢cyw=0.141K10"; and cy=0.116«10" 
ergs/cc. Knowing the magnetostriction and the 
elastic constants, we may with the help of Eq. 
(3) evaluate K, and Kg». The constant c for iron 
is taken from Akulov.® The values of the various 
constants of Eq. (1) for iron as used in this paper 


2 F. Bitter, Phys. Rev. 38, 546 (1931). 

3K. Honda and Y. Masiyama, Sci. Rep. Tohoku Univ. 
15, 755 (1926). 

4 E. Goens and E. Schmid, Naturwiss. 19, 520 (1931). 

5N. Akulov, Zeits. f. Physik 57, 249 (1929); 67, 794 
(1931). 
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are: ¢=2.15X10°; K,=—3.06X10'; 
X10’ ergs/cc; J» = 1718. 


2=2.85 


THE REPRESENTATION OF Es 


It is difficult to visualize the properties of the 
function Ey without a rather wearisome lot of 
computing. Since Eg promises to be of impor- 
tance in the study of ferromagnetism, and per- 
haps also in the study of certain non-ferromag- 
netic crystals,° I have thought it worth the 
trouble to make a few models for representative 
values of the parameters involved. Various 
methods of representation are, of course, pos- 
sible. Spherical polar coordinates were chosen. 
Thus, since EZ» defines an energy for any direc- 
tion, we may represent E, by a surface surround- 
ing a certain point in such a way that the length 
of the radius vector to the surface in any given 
direction is equal to Ey, for that direction. The 
disadvantage of this representation is that nega- 
tive radii are not desirable. The arbitrary con- 
stant in Eq. (1) is therefore so chosen that the 
figure has a reasonable size and does not become 
negative. In other words, the length of the radius 
of the model is not proportional to the corre- 
sponding energy, but the difference in length be- 
tween any two radii is proportional to the differ- 
ence between the two corresponding energies. 

The particular conditions illustrated below 
were chosen in such a way as to cover the most 
important aspects of the problem. The intention 
is to show what shapes the function E, may as- 
sume rather than to carry out the calculations 
for cases of experimental importance. This is 
especially apparent below when we speak of 
elastic compressions of 50 percent, etc., which 
are, of course, not experimentally realizable. 

The first term in Eq. (1) contains only one 
parameter c, giving rise to two different types of 
material, one, like iron, for which c>0, and 
another, like nickel, for which c <0. The remain- 
ing terms in Eq. (1) represent the influence of 
crystal distortion and of an externally applied 
magnetic field. We shall consider these last terms 
only for a crystal of iron, and only for the six 
special cases listed in Table I. 

It is to be observed that we have chosen dis- 
tortions which give rise to no volume change. 


* F, Bitter, Phys. Rev. 42, 731 (1932). 
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TABLE I. The parameters used in the illustrations of E». H is 
the magnetic field, and e is the elongation. 


H || [100] H;=H; H; =H. =0 
H || (110) Hi =H;=H/2'; H,=0 
H || (111) Hi=H; =H. =H/3! 


e \| [100] Ayw=e; A;;=Am= —e/2; Ay=Apn= eee =0 

e \| [110] Ay=Aj;=e/4; An= —e/2; Ai; =A i =3e/4; 
An= +++ = 

e || [111] Ai =Ajj=Aum=0; Aij=Ago=Ap=--- =e/2 


These distortions are not produced by the appli- 
cation of ordinary tension or compression. On 
the other hand the results obtained will be very 
similar to those resulting from true tension and 
compression, and in addition are particularly 
simple to calculate. 


THE FUNCTION E, 


Figs. 1 and 2 represent the function Eg for the 
special case A;;=A,;= +++ =0, Z=0, i.e., for the 
undisturbed ideal lattice. Since the constant c 
for nickel is not only negative, but also ten times 
smaller than for iron, Fig. 1, if constructed on the 
same scale as Fig. 2, would not deviate much 
from sphericity. The scale used for Fig. 1 is such 
that 1 cm represents 1718/500 times as many ergs 
as a like distance in Fig. 2. Since the ratio 1718/ 
500 is that of the saturation intensities, the scale 
chosen is such that a field 7 produces the same 
effect on the model for iron as on the model for 
nickel. In all the photographs here reproduced the 
plaster model stands next to a crystal model. In 
every case the axes of adjacent models are 
parallel. From Fig. 1 we see that in nickel E,» has 
maxima in the direction of the tetragonal [100] 
axes, and minima in the direction of the trigonal 
[111] axes, whereas in iron, as shown in Fig. 2, 
the maxima and minima are reversed. In Figs. 3, 
4, and 5 we have £, for iron in a field of 100 
oersteds parallel to the [100], [110], and [111] 
axes respectively. In Figs. 6, 7, 8 and 9 are il- 
lustrated the effect of compression and elongation 
on Ey in the absence of a magnetic field. Figs. 8 
and 9 are two views of the same set of figures. 
Referring to Table I, compression corresponds to 
e<0, elongation corresponds to e>0. In the 
photographs the model on the left represents 
e>0O, or elongation and the model on the right 
e<0, or compression. The stick resting on the 
crystal model in the center of each illustration 
































MAGNETIZATION OF FERROMAGNETIC CRYSTALS 657 





Fic. 1. E¢ for an undistorted nickel crystal. Fic. 2. Eg for an undistorted iron crystal. 





Fic. 3. Eg for an iron crystal in a field of 100 oersteds Fic. 4. E.for an iron crystal in a field of 100 oersteds 
parallel to a [100] axis. parallel to a [110] axis. 





Fic. 5. Ee for an iron crystal in a field of 100 Fic. 6. On the left, EZ» for an iron crystal elastically 
oersteds parallel to a [111] axis. stretched 1 percent along a [100] axis. On the right, EZ for 
an iron crystal elastically compressed 1 percent along a [100 ] 

axis. 





Fic. 7. On the left, EZ» for an iron crystal elastically Fic. 8. On the left, Eg for an iron crystal elastically 
stretched +1 percent along a [110] axis. On the right, stretched] 1 percent along a [111] axis. On the right, E¢ 
Fs for an iron crystal elastically compressed 1 percent for an iron crystal elastically compressed 1 percent along 
along a [110] axis. a [111] axis. 
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Fic. 9. A second view of the models shown in Fig. 8. 


represents the direction of elongation or con- 
traction. The models were all made for e= +0.01. 
All the models in Figs. 2—9 are made on precisely 
the same scale except the model on the left of 
Fig. 6. For this model a slightly larger additive 
constant was used to avoid negative values of Eo. 
The models in Fig. 6 depend only on ¢ and K,, 
while those in Figs. 8 and 9 depend only on ¢ and 
Ko, as may readily be verified by substituting the 
values listed in Table I into Eq. (1). Because K, 
and Kz in iron have opposite signs compression 
along a [100] axis has an effect on Ey similar to 
extension along a [111] axis. Further, as has al- 
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ready been pointed out,' E, for compression along 
a digonal [110] axis depends on Ke, K,+ Ko, 
and K,—Ks. 


MAGNETIZATION 


The magnetization curves obtained by apply- 
ing a field parallel to the elongation or compres- 
sion in Figs. 6, 7, 8 and 9 are shown in Figs. 10, 
11 and 12. The cross-hatched portions of the 
curves indicate that there is a transition from 
one minimum of Ee» to another, and that we 
therefore have a transition in some unspecified 
manner from one magnetization curve to an- 
other. Some of the magnetization curves are 
vertical for J7=0. These vertical portions also 
involve transitions from one minimum of E, to 
another. Experimental evidence already pre- 
sented® indicates that for e=0, the curves as 
drawn represent the behavior of actual iron 
crystals fairly well. As far as I know, no experi- 
mental data exist with which to check the other 
curves. 
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Fic. 10. Magnetization curves for an iron crystal elastically stretched or compressed along a [100] axis by an 
amount é, the applied magnetic field being parallel to the compression. 
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Fic, 11. Magnetization curve for an iron crystal elastically stretched or compressed along a [110] axis by an amount 
e, the applied magnetic field being parallel to the elongation. 
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Fic. 12. Magnetization curve for an iron crystal elastically stretched or compressed along a [111] axis by an amount 
e, the applied magnetic field being parallel to the compression. 
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The formulae derived from Eq. (1) represent- 
ing the curves drawn in Figs. 10-12 are: 


e|| [100 }—Fig. 10. 
H=al+bl*, 
b= 


a=(3eK,+4c)/J/,?, 

(4) 
—8c/I,'. 
e|| [110 ]—Fig. 11. 


For e>0 and for the upper portion of the curve 
for which e= —0.005 


BITTER 


H=al+bI’, a=(3eK2—4c)/I,?, 
b=8c/I,', ) 
for the lower part of the curves corresponding 
to e<0 
T=alI+6)*, a=[(3/2)(Ki+Ke)e+4c ]/T,,?, 
b= —6c/I,*. ” 
The shape of Es should be examined in some de- 


tail before applying this formula for small values 
of e. 





e|| [111 ]—Fig. 12. 


[l= 


3! 2c sin @ cos 6[2 cos? @—sin? 0]+eK2[sin 0 cos 0+2(cos? 6—sin? @) ], 





Iw 
1=I,,[(1/34) cos 0+(2/3)! sin 61. 


—sin 0+2! cos 0 


’ 


(7) 


Anyone interested in securing casts of the plaster models here illustrated should communicate 
with the Magnetic Division of the Westinghouse Research Laboratories, East Pittsburgh, Pennsyl- 


Vania. 
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New Techniques in the Cosmic-Ray Field and Some of the Results Obtained 
With Them* 


RoserT A. MILLIKAN, California Institute of Technology 
(Received February 16, 1933) 


ACH new field of scientific study goes through a period 

of development which may be likened to the develop- 
ment of a man child. Up to the age of six or seven a boy 
interests everybody intensely. In babyhood he is new and 
sweet and altogether charming because of the awakening 
qualities which bloom forth as he becomes a real person, 
and from then on up to seven or eight he is intensely 
interesting because he is revealing day by day what kind 
of a person he is likely to be. Then from say eight to 
thirteen he goes through a period in which he is in general 
best described as a whining brat who knows no law, is 
always against the government, and has no real friend 
except his mother. From about thirteen on he begins to 
find himself as a member of the human family, and by the 
time he is twenty-one is ready to don the toga and be 
admitted regularly into the ordered society of adults. 

It is now just over twenty-one years since Gockel' made 
the balloon ascent which revealed the new and quite 
unexpected fact that there were penetrating rays at an 
altitude of 4 kilometers of greater intensity than at the 
surface, i.e., at a height to which no rays could possibly 
reach, as everybody then knew, if they came from the 
earth. The cosmic rays ought therefore to be now of age. 
But the war robbed them of six or seven years of their life, 
and this may account for the fact that they are just now 
emerging from the lawless, friendless period into the period 
of orderly, recognized behavior. But until this emergence 
is more or less complete, the less public attention is focussed 
on the “ugly brat” the better it is both for him and for the 
public. In the lawless or controversial state, before the 
simple facts have been determined and sufficiently checked 
to become generally recognized by physicists, the child is 
of no interest or value to the public. Better for him to have 
no friend except his scientist mother until he learns to 
behave himself. 

I shall then have little to say today about the con- 
troversial aspects of the cosmic-ray field. These will not be 
controversial long, nor are they anything like as contro- 
versial now as the public, for some reason, seems to think. 
Indeed so far as all the major and really significant facts 
of observation are concerned, there is general agreement 
now. There are differences as regards accuracy of measure- 
ment and there undoubtedly are differences in interpre- 
tation, but so far as I myself am concerned, nothing has 
happened which alters in any essential or fundamental 
way the views which I expressed in my last comprehensive 


* Address delivered at Atlantic City on December 30, 
1932, as part of a symposium on cosmic rays. 
1 A. Gockel, Phys. Zeits. 10, 845 (1909); 12, 595 (1911). 


report written more than a year ago and printed in the 
Proceedings of the International Electrical Congress*® held 
last July in Paris. I hope very much that anybody who is 
really interested in the subject will take the pains to read 
that report before expressing or even forming an opinion, 
for within the past ten months the American newspapers 
have been quite consistently misunderstanding and mis- 
representing my findings, and I suspect that even some of 
my physicist friends, who have been too busy with their 
own work to look up the facts, have been pretty thoroughly 
misled by incorrect reports.* However, I shall not consume 
precious time by making the necessary corrections now, 
but content myself with the mere statement that the present 
report is to be regarded as a supplement to the aforementioned 
one, for that report, insofar as it goes, will actually stand 
today with scarcely a change beyond the deletion of a 
single line which represents merely a slip that should have 
been deleted anyway, even if no further work had been 
done (see below). I shall indeed present today new experi- 
ments, new techniques, and new results, but these have 
actually served in the main merely to confirm formerly 
expressed views, though adding some new and important 
experimental facts, especially in the fields of energy 
measurements and of geographical distribution. So far as 
possible, too, I shall confine myself to results which have 
been so fully checked and so often repeated that they will 
now meet, I think, with general acceptance; for this is the 
stage at which they become both of interest and of value to 
the public. 

The first important and non-controversial fact is that the 
cosmic rays ionize through the mechanism of the passage 


2? Etat Actuel de nos Connaissances sur le Lieu et le 
Mode de Production des Rayons Cosmiques, Congrés 
International d’Electricité, 1932. See also Annales de 
Institut H. Poincare. 

*Even A. H. Compton (Phys. Rev. 41, 682 (1932)) 
altogether erroneously reports Bowen and myself as finding 
but 6 ions at 16 km, and then draws conclusions from this 
fictitious number which have not the remotest relation 
to our findings. We actually have never made any esti- 
mate as to the ionization at 16 km. Our only statement 
has been that our flight ‘“‘seemed to give the mean value of 
the ionization between 5 km (where the value is 15 ions) 
and 15.5 km as 46 ions.” For three years we have repeatedly 
stated that up to 9 km Kolhérster and ourselves are in 
substantial agreement so that Dr. Compton's observations 
at about half that height have not the remotest bearing on 
16 km values. 
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with speeds approaching that of light of charged particles 
which shoot through electroscopes or other observing instru- 
ments and ionize in the well-known fashion of electron 
rays. In the term electron rays I shall include both positives 
and negatives, i.e., protons and negative electrons, What 
seemed to Bowen and myself unambiguous proof of this 
mechanism was first worked out in the spring of 1931, this 
proof being found in the observed pressure-ionization rela- 
tions in pressure electroscopes. We did not publish our 
findings until I had worked out in July and early August, 
1931, also the predicted influence of diurnal changes in 
temperature upon the observed ionization within a pressure 
electroscope. When these also checked we sent our findings 
to Nature, which published them in the October 3d issue, 
1931, but on September 14th I had presented them orally 
to Steinke and, a few days later, to Hoffmann in Europe, 
who at once agreed with our conclusions. Also, a short time 
after our first publication Bennett, Compton and Stearns, 
who had independently reached similar conclusions, pub- 
lished their findings, so that, so far as I know, this con- 
clusion, namely, that the cosmic-ray ionization that we 
observe at the earth’s surface is due to high-speed charged 
particles, is now considered as definitely proved. At least 
it seems to be universally accepted. 

However, the very significant fact that proton rays, as 
above stated, as well as negative electron rays, are the 
immediate ionizing agents and in something like the same 
proportions was first brought to light in September and 
October, 1931, by Carl D. Anderson’s measurements’ on 
the curvature of these charged-particle-rays in a very 
powerful magnetic field with the aid of a special vertical 
Wilson cloud chamber devised by himself and the writer 
more than two and a half years ago, and permitting 6-inch 
ray-tracks to be obtained in a uniform magnetic field of 
strength up to 18,000 or even 20,000 gauss. This is one of 
the new techniques mentioned in the title and one that 
extended the range of the measured kinetic energies of 
charged particles from about 15 million volt-electrons, the 
limit prior to 1931, up to the order of from two to four 
billion (10°) volt-electrons, where it now is. The possibilities 
of measurement go even higher but no higher energies have 
been found. Also much caution must be used in studying 
these nearly straight tracks because of minute sudden 
changes in direction which sometimes appear. For example, 
we have cases of nearly straight tracks which seem to show 
a very slight positive curvature before traversing a lead 
block a centimeter thick and a slight negative curvature 
after so doing—an obviously secondary influence of some 
kind, probably a minute deflection due to a close encounter 
and not a real curvature at all, and of course such an 
influence may produce either an apparent increase or 
decrease in a curvature. 

But that positive electrons or proton rays constitute 
something like half of the immediate ionizing rays, and 
that some 15 percent of these positive ray tracks are associ- 
ated with negative ray tracks starting from the same center 
constitutes the strongest sort of evidence that both tracks arise 


3 Millikan and Anderson, Phys. Rev. 40, 325 (1932); and 
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immediately from the disintegration of the nucleus of an atom 
which has been hit by a primary ray of some sort, in other 
words, that the tmmedtate tonizing agents in the cosmic rays 
are themselves secondaries released from the atoms of the 
atmosphere by non-ionizing primaries. I say non-ionizing 
primaries because we have placed strips of lead a centimeter 
thick in the middle of the cloud chamber and have re- 
peatedly directly observed ionizing tracks, some times 
negatives, some times positives, some times both, simul- 
taneously to spring out of the lead when no visible ionizing 
track enters the lead.* This is direct, visible proof of what 
Hoffmann, Geiger and others long ago showed by placing 
thin lead sheets above a cosmic-ray electroscope and noting 
that the ionization in the electroscope was increased 
instead of decreased by the interposition of the lead, a 
phenomenon long known with gamma-rays and meaning 
merely that new secondary charged particle rays are 
created in the lead by the absorption by the atoms of the 
lead of the incoming gamma-rays. 

Dr. Workman, one of the National Research Fellows in 
our laboratory, has recently completed an elaborate study 
of the ionization produced in an electroscope by the 
hardest gamma-rays as a function of the thickness of a 
given substance (expressed generally in terms of number 
of electrons per square centimeter) interposed between the 
source of the gamma-rays and the ionization chamber 
Each material, of course, has its own characteristic curve. 
If the aluminum curve, for example, is being plotted and 
successive sheets of aluminum have already been intro- 
duced, and if, from a given thickness on, lead sheets are 
used in place of aluminum, the new points at once leave 
the characteristic aluminum curve, rise above, then fall 
below it, and slowly move over to the characteristic lead 
curve, and reach it as soon as the secondaries released from 
the originally interposed aluminum sheets are no longer 
able to get through the added lead sheets. Now, Schindler* 
in Steinke’s laboratory has recently published exactly the 
same sort of curves taken with cosmic rays instead of with 
gamma-rays, and has obtained exactly the same relation 
that Workman has obtained with gamma-rays, thus 
demonstrating again and in a new way that precisely as in 
the case of gamma-rays the ionizing cosmic-ray-particles 
are indeed secondaries released by incoming non-ionizing 
primaries. 

But the third and perhaps the most complete demon- 
stration of this conclusion is furnished by Carl Anderson’s 
measurements of the energies of the actually observed 
cosmic-ray particles; for these measurements show that 
the majority of these energies lie below 600 million volts, 
Dr. Anderson and I having published the estimate that 
not more than a tenth of these tracks reach appreciably 
above the billion (10%) volt range. But now we know with 
great definiteness and great certainty, too, the maximum 
distance of penetration of billion volt-electron rays through 
the atmosphere; for such rays, be they positives or nega- 
tives, (for at these energies the differences are small) make 
at least 30 ions per cm of path (45 ions seems now a better 
mean), and each ion removed from the molecules of air 


4 Heinz Schindler, Zeits. f. Physik 72, 625 (1931). 
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represents an expenditure of 32 volt-electrons of energy; 
or such a billion volt-electron as we are considering, and in 
the case of electrons the same is true of all lower speeds so 
long as the energy is as high as a million volts (in the case 
of protons the penetration at low speeds is much less)— 
loses at least a thousand volts per cm of air path, which 
means that the maximum range of a billion volt-electron 
in air cannot be more than a million centimeters or 10 km, 
more accurately 9.6 km. Indeed Heisenberg® makes it but 
a third that much for billion volt particles. But the depth 
of a homogeneous atmosphere of the density existing at 
ordinary temperatures and pressures is 1033 divided by 
0.00117 and therefore equals 8.8 km. In other words, a 
billion volt-electron, if it happened to come in normally to 
the earth surface, might possibly just get down to sea level 
and then pass through a roof and an electroscope-wall 
having a combined thickness equal to a cm of brass. Since 
the rays come in from all directions, even if all of them were 
billion volt-electrons, a wholly negligible number of them 
could possibly reach such an observing chamber. Since 
something like nine-tenths of the actually observed ionizing 
rays have energies under this value, nearly all of the rays 
observed in such a chamber must then be secondaries, not 
primaries. This conclusion seems to me altogether inevi- 
table from the foregoing facts alone. The only possible 
escape is to say that our measurement of the energies is 
wrong. It might be slightly wrong in the higher volt range, 
but it cannot be wrong in the only significant range which 
is the lower volt range, as anyone can see if he will examine 
our photographs. Further, Kunze* in Fuchtbauer’s 
laboratory in Rostock, has duplicated our measurements 
in Germany and gets the same distribution of energies as 
we do. He lists a negligible number of energies up to 4 or 5 
billion volts, but, as heretofore indicated, this listing 
perhaps represents less caution about his high energy 
measurements than we use, though this is of course quite 
unimportant for our present considerations. 

So far as I can see then the only escape from this argu- 
ment that these ionizing particles observed at sea level 
originate practically wholly in our atmosphere is not only 
to deny the validity of these direct measurements of 
energies but to postulate, as some have attempted to do, 
incoming charged particle energies at least ten times, or 
even a hundred times, higher (10 or 10" volts) for the 
usual calculations require these energies if the rays are to 
get through the earth’s magnetic field into the equatorial 
regions at all. This postulate clashes head on and altogether 
disastrously both with the observed distribution of the 
cosmic rays in altitude and in azimuth (inclination with 
respect to the vertical). For 10" volt or 10" volt incoming 
rays of this sort will produce an atmospheric ionization 
very nearly independent of altitude, instead of increasing 
very rapidly as the facts require; and as to azimuth, such 
rays would produce a uniformity of distribution of the 
observed ionization at sea level over the whole celestial 
dome. But for at least seven years I have been testing this 
last point in mountain valleys and basins, and I have never 
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found any measurable decrease in the observed ionization 
when the horizon was screened by mountains up to 30 
degrees or even 40 degrees from the horizon. In other 
words, practically all the rays come in at sea level through 
a vertical cone of apex at the observer and of 100 degrees 
or at most 120 degrees opening. 

But I am perhaps spending time unnecessarily upon the 
point of the atmosphere origin of the observed ionizing 
cosmic-ray particles, for I do not see how it can be doubted 
for a moment by any one who is at all generally familiar 
with cosmic-ray facts. There is, however, one further closely 
related point which should be mentioned. For the main 
facts which have at any time led anybody to postulate 
incoming electrons of an energy of 10" or 10" volts are, 
first, so-called Geiger counter coincidence measurements 
on absorption coefficients, and second, the aforementioned 
geographical distribution of the rays. With respect to the 
first, I have been pointing out for two years in Pasadena 
seminars, in the Rome congress on nuclear physics in 
October, 1931, in New Orleans last Christmas at the 
A.A.A.,S. meeting, and in the report for the Paris Electrical 
Congress, that these counter experiments never in my 
judgment actually measure the absorption coefficients of 
anything. I shall presently show that no appreciable 
number of these observed ionizing particles ever go through 
more than 30 cm or at most 60 cm of lead, and yet both 
Regener and Cameron and I have proved that the cosmic 
rays penetrate through the equivalent of more than 20 
feet of lead. These figures cannot both be correct without 
carrying with them the conclusion that the primary rays 
at sea level and below are not charged particles. The 
mechanism of these coincidences is, I think, as follows: 

When two counters are separated only by a foot or two 
of interposed air the coincidences are due indeed wholly to 
these observed electrons (using now this word to cover both 
negatives and positives) which shoot through both 
counters. When, however, 15 cm of lead are interposed 
between the counters, as in Mott-Smith’s experiments, the 
electrons of energy five hundred million volts or more, of 
which there are actually a considerable number, do indeed 
shoot in succession through the upper counter the lead and 
the lower counter just as is usually assumed: but the number 
of coincidences thus produced is now augmented by new 
coincidences due to the practically simultaneous emission of a 
succession of secondary particles released along the path of a 
photon by Compton and other encounters, the chance for the 
emission of such secondaries being now enormously in- 
creased by the large amount of new matter interposed 
between the counters. The result is that the number of 
coincidences with the interposed lead is very much larger 
than that which corresponds to the absorption coefficient 
of the charged particles that get through. Messrs. Pickering 
and Sharp in our laboratory have just now obtained certain 
confirmations of the correctness of this explanation. It will 
be seen that if this view is correct these so-called absorption 
coefficients when thick layers of matter are interposed 
mean then nothing at all, and the conclusion about particle 
energies that have been drawn from them are quite invalid. 
The conclusions to be drawn from latitude effects will be 
treated below. 
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The foregoing interpretations gain further support from 
the new measurements made in the Norman Bridge 
Laboratory by Carl Anderson on the actual losses in energy 
experienced by cosmic-ray secondary electrons in going 
through lead sheets one or more centimeters thick. These 
measurements are upon the changes in curvature of the 
same electron tracks (plus and minus) above and below 
strips of lead a centimeter thick through which they have 
passed. The average loss in energy is 35 million volts per 
centimeter, which agrees within less than the uncertainty 
of measurement with Heisenberg’s computed value of 39 
million volts per centimeter. It will be seen from this that 
a billion volt ray has a range in lead of but 30 cm. The 
great majority, then, of all coincidences obtained through 
more than this thickness of lead are due if this theory is 
correct to successive secondaries released along the path 
of a photon. 

Measurements of the type I have been considering in 
the foregoing paragraph yield other important information; 
for even the straightest of these tracks always show measure- 
able changes in direction in traversing the lead. The statistical 
way in which these deviations vary with incident energy 
furnishes a new way of measuring energies, supplementary 
to the method of magnetic-deviability, and applicable to 
higher energies. This method alone seems to show that no 
tracks corresponding to energies of 10" or even 10" volts 
are actually present among those thus far worked with. 
Really quantitative conclusions, however, must await 
fuller statistical data upon such energy losses in traversing 
lead. 

With all these lines of approach, and still others to 
follow, the two conclusions that seem to me no longer 
controversial are thus far (1) that the observed ionization 
at sea level is due to charged particles, and (2) that 
certainly more than nine-tenths of these particles are 
secondaries formed in our atmosphere. A‘ third thus far 
non-controversial result has to do with the general nature 
of the curve representing the way in which the ionization 
varies with altitude. I have repeatedly pointed out for the 
past two years, both in speech and in print, that the curve 
obtained by Dr. Cameron and myself relating the observed 
ionization to depth beneath the top of the atmosphere 
and carried by us with great accuracy up nearly to 5 km 
runs sufficiently closely into Kolhérster’s balloon observa- 
tions, which he carried up to 9 km. It also runs into the 
more accurate airplane observation up close to 9 km 
recently made by Mott-Smith and Howell’ and Millikan 
and Neher (see below). In other words, Kolhérster and 
ourselves are in reasonably close agreement up that far. 
I have further based the conclusion that this exponential 
rise, observed by all of us, and reaching in its upper part 
a value close to » = 0.6 per meter of water could not possibly 
continue to the top wholly upon the integrated value of 
the ionization in the 15.5 km balloon flight with self- 
recording instruments which Bowen and I made in 1922, 
and I have repeatedly drawn the important conclusion 
from this and other phenomena that the incoming rays 
could not be in complete equilibrium with their secondaries 
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upon entering the atmosphere, since incoming rays in such 
equilibrium should yield such a continuous exponential 
rise. The absence of this exponential rise, and of the equilib- 
rium of the primary beam with its secondaries upon entering, 
is the sole conclusion which I have ever drawn from our 
1922 flight. If the entering rays are entirely unmixed with 
secondaries upon entering, i.e., if they constitute a pure 
photon-beam, there should of course be a maximum in 
the depth ionization curve somewhere on the way to the 
top, and a year or more ago when Piccard reported a lower 
observed ionization at 16 km than at 9 km, I thought that 
this maximum must be found within those limits and so 
stated, chiefly on the basis of Piccard’s briefly-reported 
findings, in my review for the Paris congress. This is the 
one line that should now be deleted from that review since 
Piccard has now in his final report* changed his judgment 
about these two readings. But I further reported in that 
review that current cosmological theories ‘‘seemed to lead 
to rays which should be in equilibrium with their second- 
aries (particle rays) upon entering the atmosphere; that 
one must then regard this question as not yet settled’’ and that 
Bowen and I were repeating these high balloon flights for 
the sake of get-more accurate data near the top of the 
atmosphere, upon the degree of purity, or freedom from 
secondary particle rays, of the incoming beam. 

With the second of our new techniques, for which Dr 
Neher’s skill is largely responsible, and the cooperation of 
the Weather Bureau, which is also very important, Bowen 
and I have been repeating then during the past summer 
our high balloon flights with self-recording instruments, 
and in the meantime Regener® has succeeded in getting 
one such flight, too, in which he reached an even higher 
altitude than ours, his altitude record being 28 km while 
ours is 26 km, though the returns from our last flight are 
not yet in, and we of course hope that it reached his altitude 
or better. However, the two flights that are already in and 
computed check reasonably Regener’s up to 20 km, where 
our electroscope records stop, and all three of these flights 
seem to confirm the conclusion of lack of equilibrium which 
we drew from our 1922 flight, though the temperature 
effect on the electroscope used in 1922 now turns out to 
have been considerable. How beautiful is the new electro- 
scope technique in the matter of freedom from temperature 
effects is shown by the fact that a change in temperature 
from minus 77 degrees C to plus 23 degrees C causes a 
change of not over 2 percent in the observed readings in 
our test chamber. 

This lack of equilibrium of the incoming rays with their 
secondaries,—a subject discussed more fully in the report 
to the Paris Electrical Congress—is attested first by the 
shape of the depth ionization curve as found now both by 
Regener and by Bowen and Millikan, for this curve is not 
even concave upward as the upper part of the atmosphere 
is traversed as it should be, and very strongly so for rays 
coming in from all directions, in case these rays are in 
equilibrium with their secondaries. Instead of this the 
curve is actually found to have bent over so that it is 


8 A. Piccard, Comptes Rendus 195, 71 (1932). 
*E. Regener, Naturwiss. 20, 695 (1932). 
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somewhat concave downward. It is attested in the second 
place by the fact that the ionization at the highest point 
reached (computed by assuming that the coefficient of 
absorption actually found by all of us, Kolhérster, Regener, 
Bowen and Millikan, Mott-Smith and Howell, at about 
9.0 km, namely about 0.6 per meter of water, maintains its 
value unchanged to the top, as it should for incoming rays 
in equilibrium with their secondaries) is actually not half 
the computed value, the computation of course being made 
for the actual case of rays which come in from all directions. 
Further, if the observed coefficient 0.6 is the composite 
result of rays of different hardnesses, like helium rays, 
oxygen rays, etc., as Cameron and I have suggested, the 
discrepancy with the observed ionization becomes much 
greater. Regener himself, then, draws from his curve 
precisely the conclusions which we drew from our 1922 
data, as well as from our present data, so that the whole 
field of ionization at very high altitudes is one in which 
there is as yet no disagreement either as to experimental 
findings or as to conclusions drawn therefrom by the 
observers themselves. It should be noted in this connection 
that the discovery that the secondaries are positives as well 
as negatives means that even with a quite small admixture 
of incoming particle-secondaries—though no approach to 
equilibrium—these positives, released from the nuclei of 
atoms by the relatively soft incoming primaries which 
seem to exist in the upper layers, say of 27 million volt 
energy, would ionize heavily and thus tend to push up the 
maximum in the ionization-altitude curve to close to the 
top of the atmosphere. The data close to the top are, how- 
ever, not yet reliable enough to make quantitative con- 
clusions possible. 

But why is it so important to determine whether the 
incoming rays are in equilibrium with their secondaries or 
not? Because, if the cosmic rays originate in the stars they 
should be in such equilibrium, the primaries being mixed 
with secondary electron rays just as they are at sea level, 
and if they are not in such equilibrium, then the conclusion 
that the primaries originate in interstellar space gains 
added support, though even so no definite proof, since even 
if they do so originate they might still be in equilibrium 
(see above). In any case the incoming rays should be mixed 
to some extent with their secondaries, as I pointed out in 
1931, since they must at least come through some of the 
matter, which small though it be, is distributed throughout 
space. I also pointed out in the paper written for the Paris 
Congress, before latitude effects had become in any way a 
controversial subject, that ‘On ne peut se reposer entiére- 
ment sur ces seules indications”; that new light should be 
thrown upon this important question by looking more 
carefully for such latitude effects, especially at high eleva- 
tions and high latitudes; for at sea level under no assump- 
tions which seemed either likely or natural could the 
secondary rays make their effects felt save as they had 
energies much greater than a billion volts, and no very 
considerable number of such energies were either then or 
now thought to exist. In other words, practically all such 
secondary incoming electron rays would in any case, no 
matter where they entered the atmosphere, be filtered out 
before the earth’s surface was reached. But, at an elevation 
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of seven kilometers, for example, or 23,000 feet, secondary 
entering rays of energies greater than 400 million volts, 
and many such were known to exist, might make their 
effects felt in electroscopes so far as penetrating power 
alone was concerned, and according to my own reasoning 
and Epstein’s calculations” the earth’s magnetic field 
should let these high altitude effects through in the high 
magnetic polar latitudes, but keep them out of the equa- 
torial magnetic latitudes. 

I did, indeed, suggest in the Paris paper that incoming 
secondary charged particles might produce latitude effects 
even in temperate latitudes through the possible agency of 
photons arising from the nuclear impacts of such high- 
speed electrons plunging through the upper part of our 
atmosphere, and since no such latitude effects were found 
between Pasadena, California and Churchill, Manitoba, 
the very region in which such effects were to be expected, 
I stated that there seemed to be as yet no evidence that 
the incoming rays were mixed in any measurable degree 
with their secondary electron rays. This line of approach 
involves the assumption that very penetrating gamma-rays 
can be produced by the passage of high-speed electrons 
through air, an assumption pretty certainly contrary to 
fact, for, first, although very soft gamma-rays are indeed 
produced when K electrons are knocked out of atoms by 
high-speed bombarding electrons, and although somewhat 
harder gamma-rays might be occasionally produced by 
similarly knocking out nuclear electrons, I think no gamma- 
rays capable of penetrating kilometers of atmosphere, as 
would be necessary if such secondary gamma-rays become 
observable at sea level, are to be expected from the plung- 
ing of high-speed electrons through the light atoms of the 
atmosphere. The energy available from the mere rearrange- 
ment of the parts of such a smashed atom is wholly inad- 
equate. Even if such secondary gamma-rays might thus be 
produced they could not account for Carl Anderson's 
observed effects no matter what the energy of the incoming 
particles, for his ray-tracks in the great majority of cases come 
in from above, while 10 volt electrons shooting through the 
atmosphere should produce at the earth's surface photons 
Shooting out equally in all directions from the bombarded 
nucleus in case it be assumed that they can produce such high 
energy photons at all. 

Nevertheless, if the incoming rays are actually mixed to 
any appreciable extent with their secondary electron-rays 
of the observed energies the effects of the latter should of 
course be observable at high altitudes in polar latitudes and 
absent in the equatorial latitudes, and accordingly more 
than a year ago, i.e., in 1931, I laid before the Carnegie 
Corporation of New York a definite and detailed program 
for the carrying out of this high altitude exploration 
program and also described it explicitly in the afore- 
mentioned report for the Paris Congress long before any 
latitude controversy was ever dreamed of. I was not 
primarily interested in further sea-level observations 
because, first, the foregoing reasoning showed that if any 
magnetic latitude effects could be found at sea level they 
would of necessity be small, and second, Cameron and I 


10 Paul S. Epstein, Proc. Nat. Acad. Sci. 16, 658 (1930). 
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had in 1927 carried sea-level observations from latitude 34 
north to 17 south and found no variation which was outside 
the limits of our experimental uncertainty, which we 
estimated at 6 percent. Further, the next year I had helped 
Dr. Gish to calibrate electroscopes which were taken on the 
good ship ‘‘Carnegie’’ and with which 500 observations 
were made in the equatorial latitudes from 20 degrees 
north to 20 degrees south, and an equal number from 20 
degrees to 60 degrees south, and also half that number from 
20 degrees to 60 degrees north, without bringing to light 
any sea-level latitude effects which were greater than the 
observational uncertainty which I estimated again, upon 
privately getting this data from Dr. Gish, to be of the 
order of 6 percent. However, these Carnegie observations 
showed, if anything, slightly lower equatorial values than 
higher latitude values, the mean of 507 readings between 
20 north and 20 south being 3.32 ions/cm*/sec., of 262 
readings between 20 north and 60 north being 3.39, while 
the mean of 509 readings between 20 south and 60 south 
was 3.44 ions/cm’/sec. The estimated zero was here about 
2 ions/cm*/sec. The readings in the equatorial belt are 
indeed 5 or 6 percent lower than those in the north tem- 
perate belt, but those in the south temperate belt differ 
from those in the north temperate by nearly as much. 
Further all of these differences might have been due to 
barometric differences, for which at that time dependable 
correction could not be made. It was for these reasons that 
in 1930 I made an elaborate and accurate comparison of 
sea-level intensities at Pasadena (lat. 34) and at Churchill 
within 700 miles of the North magnetic pole, for here I 
expected the maximum differences to appear, but I found 
no difference that exceeded my observational uncertainty 
which was stated in these measurements to be not over 1 
percent. 

For these 1930 observations from latitude 34 northward 
both at sea level and on mountains were taken with 
improved, very sensitive electroscopes and were very much 
more accurate than our South American ones. Bennett and 
Dunham," as well as Neher and myself, have this last 
summer completely confirmed these findings, not only at 
sea-level but on mountains from five to nine thousand feet 
high, so that there is now entire agreement upon the only 
point having to do with latitude effects upon which I have 
at any time been insistent or even claimed to have made 
accurate measurements (see also below) namely that from 
Pasadena northward, within the limits of our observational 
uncertainty, there is no effect of latitude so far as observations 
on terra firma, high or low, have gone. It is to be noted, too, 
that according to Epstein’s'® and Lemaitre’s and Vel- 
larta’s'® calculations, these high latitudes, north or south, 
are the only regions in which latitude effects on intensities 
can be expected as a result of the action of the earth’s 
magnetic field on incoming charged particles of the ob- 
served energies. There is no doubt at all about the abun- 
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also Nature 128, 704 (1931) for Lemaitre’s superradioactive 
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dance of secondary charged particles up to energies of 500 
million volts. If the effects of such secondary incoming 
rays do not appear about the magnetic poles, the alter- 
native inferences will be (1) that the incoming primaries 
are not appreciably mixed with their magnetically deflect- 
able secondaries, or (2) that we have not extended our 
observations to high enough altitudes in the atmosphere 
or to high enough latitudes to bring these effects to light, 
for the earth’s atmospheric skin ts so thin that no latitude 
effects of secondaries formed in our atmosphere can be 
expected. 

A third possible hypothesis might be made, namely, that 
the earth’s magnetic field, like the sun’s, does. not extend 
out as far into space as that of such a theoretical magnetic 
dipole should. This would shift the expected effects of the 
incoming secondaries to lower latitudes, where we should 
find them if only we carry our observations high enough. 

As to the mere observational facts at sea levels divorced 
from all theorizing, it should be reported that Clay" has 
two or three times during the past five years reported from 
10 to 20 percent of equatorial lowering on voyages from 
Batavia to Europe via the Suez Canal, lowerings which he 
has related to magnetic latitudes and that A. H. Compton" 
has reported similar results obtained in the western hemi- 
sphere. On the other hand, E. Oeser™ of the University of 
Géttingen has just repeated the sea-level observations 
between latitude 50 degrees north and latitude 7.1 degrees 
north (Hamburg to the Panama Canal) and reports no 
change whatever in his electroscope No. 1, that was air- 
tight, but a lowering of 17 percent in the equatorial regions 
for a second electroscope (No. 2) which was definitely 
proven to be leaking and which therefore contained more 
rarefied air in the hotter areas than in the colder ones. But 
the most significant and dependable measurements of all 
have just been made by Dr. Neher, who has carried one of 
our new sensitive recording electroscopes by sea from 
Pasadena (34 N.) to Mollendo (17 S.) and back without 
finding at sea level any changes appreciably larger than 
the limits of the uncertainties indicated above, namely, of 
the order of 6 or 7 percent. However, the mere uncertain 
indications previously found of about this amount of 
equatorial lowering are now confirmed by these quite exact 
and dependable readings, which show an increase in coming 
from Mollendo (17 S.) to New York (41 N.) of 6.8 percent 
correct to one or two percent. 

With respect to all my former sea-level observations, | 
have always been very careful to make no other statement 
than this: First, that “I myself have never been able to 
find any sea-level effect of latitude which was outside the 
limits of my observational uncertainty.” That statement | 
shall obviously never have to retract, and it cannot 
properly be made the subject of controversy. But, further, 
my estimated uncertainty of 6 percent between Pasadena 
(34 degrees north) and Mollendo (17 degrees south), and 


3]. Clay, Proc. Roy. Acad. (Amsterdam) 30, 1115 
(1927); 31, 1091 (1928); 33, 711 (1930); and J. Clay and 
H. P. Berlage, Naturwiss. 37, 687 (1932). 

14 A, H. Compton, Phys. Rev. 41, 111 (1932). 

16 E, Oeser, Zeits. f. Geophysik 8, 242 (1932). 
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of 1 percent between Pasadena and Churchill may of 
course some time have to change, but I think there will 
be general agreement that with the observations of Millikan 
and Cameron, the Carnegie Institution of Washington, 
Oeser and Millikan and Neher all in reasonably good 
agreement it would be quite premature to change them 
appreciably now. 

But with respect to high-altitude latitude effects, as 
already indicated, | have thought them so likely as to 
publicly state in several addresses in Europe, delivered in 
October and November of 1931, and in the aforementioned 
article written in 1931, that I was undertaking a rather 
elaborate high altitude latitude survey. The carrying out 
of such a survey required the development of the second 
new technique referred to in the title of this paper, for the 
observations had obviously to be made in airplanes which 
can hold the same level for long periods, and they had to 
be as accurate as those made in a laboratory. It is largely 
Dr. Neher’s extraordinary skill which has made the success 
of this program possible. The instruments are self-register- 
ing and work just as well in an automobile or a railway 
train running over rough roads, or even in a diving airplane, 
as in a laboratory. They yield the most accurate airplane 
observations of which I myself have any knowledge, though 
I have assisted in making many inaccurate airplane tests. 
They are very much more accurate than any balloon 
observations made by ourselves, and I think these compare 
well with the balloon observations made by others. 
Further, any one who henceforth wishes to check for him- 
self our results need only come to our laboratory and 
measure up our films, which of course are without prejudice 
or preconception of any kind. 

There are two main results of these high altitude airplane 
tests that I wish to report. With the aid of Colonel Arnold 
and the staff at March Field (latitude 34), we took up one 
of our new electroscopes inside a four-inch lead shield to 
an altitude of 21,000 feet, or about 7 km, and then we 
repeated the observations without the lead, so that we 
now know exactly how much of the rays at each altitude 
are cut out by the lead, not roughly but quite accurately. 
At 21,000 feet the ionization inside the lead is less than 
one-third what it is with the lead removed. This means 
that more than two-thirds of the ionizing secondary rays 
existing at that altitude are unable to pass through four 
inches or 10 cm of lead. This means in turn, since many 
new secondaries are of course produced within the lead, 
that very much more than 66 percent, probably as much as 
80 percent of the secondary ionizing rays existing at 21,000 
feet, have insufficient energies to carry them through 10 cm of 
lead, and that means energies below about 350 million volts. 
This checks very beautifully our conclusion that practic- 
ally all the ionizing rays are secondaries produced in the 
atmosphere, for 350 million volt rays cannot possibly 
penetrate one-third the thickness of the atmosphere. This 
way of finding the distribution of the energies of the ion- 
izing particles cannot, I think, be seriously in error. It 
shows that at least 35 percent of the ionizing particles 
existing at sea level, to take another example, have energies 
under 350 million volts and this checks well with Carl 
Anderson’s direct measurement of the distribution of 
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particle-energies. Again, Mott-Smith and Howell's high 
altitude measurements’ tell precisely the same sort of a 
tale, for at 25,000 feet these observers report that a one- 
inch lead screen cuts off 40 percent of the ionization. This 
means that certainly much more than half the particle- 
rays existing at 25,000 feet have an energy under 90 million 
volts. All these new airplane observations, then, confirm 
most emphatically the banded structure of the cosmic 
rays, upon which I have been insisting since 1925. Indeed, 
the rapidly progressive softening of the rays with altitude, here 
brought sharply to light, is I think inexplicable on any charged 
particle theory as to the nature of the primary rays. 

The other result of the airplane tests relates to the high 
altitude geographical findings. Exactly the same sort of 
flights in which the pilot flew for an hour at 10,000 feet, 
an hour at 14,000, an hour at 19,000, and an hour at 21,000, 
were made first at March Field (latitude 34, magnetic 
latitude 50 reckoned from pole position at 70°); then at 
Spokane (latitude 47, magnetic latitude 64), where the 
flights were most skillfully made by Major Breene, then at 
Cormorant Lake, Manitoba (latitude 55, magnetic latitude 
75), where Flight Commander Gordon of the Royal 
Canadian Air Force generously and skillfully flew for us. 
Also, a similar flight was made in Peru (latitude 17 south, 
magnetic latitude 2° N.), and still another, with the assist- 
ance of the U.S. Army at Panama, 7.1° N., magnetic lati- 
tude 24.1 N. In these flights the electroscope is automati- 
cally recharged once every 15 minutes, and it is so sensitive 
that it moves over the 100 divisions of the scale in 15 
minutes at sea level. The Cormorant Lake and the Spokane 
curves show no differences whatever, and yet this is where the 
magnetic effect of the earth's field on negative electrons of 
energies up to 10° volts should be most felt. Its failure to 
appear there may perhaps be considered as pointing again 
toward the absence of appreciable numbers of secondary 
deflectable rays of the above-mentioned, directly observed 
energies in the incoming beam, but an alternative explana- 
tion is that the earth's magnetic field is not as strong as is 
assumed in the Epstein, Lemaitre calculations. 

Quite unexpectedly, however, the results at March Field 
at 21,000 feet are 11 percent lower than those taken at the 
same altitude at Spokane and at Cormorant Lake. With 
respect to the high altitude tests in the equatorial belt, 
the Peru flight and the Panama flight are in perfect 
agreement and yield distinctly lower high altitude readings 
than any of the others, 40 percent lower at 21,000 feet than 
the Spokane and Cormorant Lake readings. It is interesting 
to observe, also, that this diminution in the high altitude 
intensities of the cosmic rays in the equatorial regions 
actually stands out clearly from the published readings 
which Cameron and I took in the High Andes in 1926,'® 
though we did not then call attention to it, nor ourselves 
assign to it special significance; for these readings were 
taken before we had developed sensitive high pressure 
electroscopes and we were unable at that time to lay much 
claim to precision. Further, we had taken but a dozen or 
so under-water readings at Lake Maguilla at an altitude 
of 15,000 feet when, as indicated in our published article, 


© Millikan and Cameron, Phys. Rev. 31, 165 (1928). 
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leaks developed in our electroscopes. We therefore relied 
mainly for our latitude findings upon our sea-level readings, 
of which we had a great many and which were quite free 
from suspicion as to possible effects of leaks. However, in 
view of the thinness of the earth’s atmospheric skin and the 
universally assumed character of the earth’s magnetic field 
no latitude effects due to the earth's magnetic field could 
possibly be expected at high equatorial altitudes if they failed 
to appear at sea level. This is why we based our conclusion 
as to a lack of a latitude effect chiefly upon our sea-level 
readings and drew one single ionization-depth curve 
through all our points taken up to that time at high alti- 
tudes both in the northern hemisphere and in Peru and 
Bolivia and assumed that the points that fell off this curve 
represented observational uncertainties, though, as any 
one can now see by looking at them, these points actually 
do show clearly the differences bet ween readings at a given 
altitude in Bolivia and in the United States. For our 
published composite curve embracing all readings taken up 
to the end of 1926" is actually considerably lower at its 
upper end, e.g., at 7 m beneath the top of the atmosphere, 
than is the corresponding curve taken in the northern 
hemisphere either in 1925'* or in 1931.7 Indeed, the mean 
values of the ionization at a depth of 7 m beneath the top 
of the atmosphere is 16 percent higher in the Muir Lake 
data than in the Lake Maguilla data, and when we compare 
the more accurate 19317 curve we find just as does Comp- 
ton, a difference of more than 30 percent at this same 
altitude. In other words, from this last curve the reading 
corresponding to a depth of 1 m beneath the surface of 
Lake Maguilla should have been more than 6 ions instead 
of the observed 4.1 ions. In other words, the published 
Bolivia readings are actually at least 30 percent lower than 
the readings taken at the same altitude in the central part of 
the United States. | have often discussed with others this 
difference between our curves taken in South and North 
America, and it was chiefly because of this discrepancy and 
of the aforementioned readings taken on board the Car- 
negie that I decided in 1931 to make what I regarded as a 
crucial test of this latitude question by comparing inten- 
sities at Pasadena and at Churchill, the only latitude 
stretch in which I estimated that the effect of the earth’s 
magnetic field could cause large differences at any altitudes. 
In a word, then, the lowering of the observed cosmic-ray inten- 
sities at high altitudes in the equatorial regions is consistent 
with all the high altitude observations which have yet been 
taken by anybody, so that there is no chance for controversy 
about this point. Only at sea level is the observational data 
not yet in quantitative agreement. In the matter of in- 
terpretations there are quite naturally differences of 
opinion. 

With respect to these interpretations it is quite clear that 
this new high altitude data in different geographical areas 
cannot easily be brought into consistency merely by 
invoking the earth’s magnetic field, for there are real 
difficulties which a magnetic interpretation encounters. 
The first is to explain the fact that the small sea-level dip 
in intensity occurs only in quite low, i.e., in equatorial 


17 Millikan and Cameron, Phys. Rev. 37, 244 (1931). 
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latitudes. According to a hypothesis which Epstein and | 
have considered and which he will presently publish in 
quantitative form, it is sought to explain this dip as follows: 
The photon hypothesis as advanced by Millikan and 
Cameron is adopted to account for the whole of the ob- 
served ionization in the equatorial belt, for the earth's 
magnetic field would in any case prevent the penetration 
of particles of the observed energies into this belt. But of 
the ionization observed at sea level between Pasadena 
and Churchill, only 93 percent would be due to incoming 
photons while to account for the 7 percent equatorial dip 
it is assumed that there is a small admixture of charged 
particles superimposed upon these photons. From Epstein’s 
point of view these particles might be looked upon as 
having the same relation to photons as beta-rays have to 
gamma-rays in radioactivity, or they might be secondaries 
produced by the incoming photons and therefore entering 
our atmosphere as a beam either in partial or total equi- 
librium with its secondaries. But to enable the earth’s 
magnetic field to remove them only in equatorial latitudes 
and leave the sea-level intensities from Pasadena to the 
pole constant, Epstein is obliged, following his own and 
Lemaitre’s computations, to postulate that these particles 
have energies as high as from three to eight billion volts. 
Such particles, however, should have so high a penetration 
that, in the first place, they should be observable at all 
angles from the vertical clear down to the horizon, which 
they do not seem to be, and, in the second place, they 
should appear in greater numbers than they seem to do in 
Dr. Anderson's direct measurement of energies. 

The third difficulty is as follows: Although the assumed 
Millikan-Cameron photon hypothesis can take care of the 
larger part of the enormous rise in the intensity of ioniza- 
tion with altitude, these high energy entering particles must 
take care of the whole difference observed at any and all 
altitudes between the ionization in the equatorial belt and 
the temperate belt. This is a difference which near the 
Equator increases from 7 percent at sea level to 40 percent 
at 21,000 feet, but since the actual ionization at 21,000 feet 
is seventeen times that at sea level and the percentage of 
particle-rays is, as above stated, five times that at sea level, 
the ionization at 21,000 feet due to the incoming particle- 
rays alone must be 5X17=85 times greater than the 
incoming-particle-ionization at sea level. But the particle- 
rays which must produce this sort of ionization are also 
required to get through the earth’s magnetic field and 
reach the earth’s surface between Balboa (latitude 7.1) 
and Pasadena (latitude 34), and it is this condition that 
limits their energies to the range from three to eight 
billion volts. If such high energy rays are to produce 
anything like the observed ionizing effects, then the 
particle energy necessary to just get through the earth’s 
atmosphere must be assumed to be some five or six billion 
bolts, instead of but about one billion as calculated above. 
This seems to be somewhat higher than the present rather 
incomplete experimental situation allows. If it should be 
found permissible the three difficulties mentioned might 
all disappear. 

But there is still a fourth difficulty. Since incoming 
particle-rays are continually losing energy along their 
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paths their energy or penetrating power should increase, 
not decrease, with altitude, but as was shown sharply by 
the foregoing measurements of Millikan and Neher in 
their airplane flights made with and without their 10 cm 
lead shields the mean penetrating power actually decreases 
rapidly with altitude, being about one-third as much at 
21,000 feet as at sea level. This behavior defies any natural 
interpretation on the basis of incoming charged particles 
alone, but fits nicely the theory of incoming photon bands 
made up of He, O, Si, Fe, etc., rays. Now, according to 
the foregoing hypothesis as to the effect of the earth's 
magnetic field, the ionization in the equatorial belt is due 
wholly to such bands of photons, but in the temperate 
belt, as is shown by the difference between the readings 
in temperate and equatorial latitudes, as much as 30 
percent of the ionization at 15,000 feet arises from the 
3 to 8 billion volt incoming particle rays. Hence, at 
15,000 feet a much larger percentage of the observed rays 
should get through the 10 cm lead shield in the temperate 
and polar belts than in the equatorial belt. Actually 
Neher and I found by quite accurate airplane measure- 
ments in the United States at 15,000 feet 38.7 percent of 
the ionization inside the 10 cm lead shield that we found 
outside it and in the equatorial belt Neher got the same 
way 39.2 percent, i.e., the rays at 15,000 feet are of very 
nearly the same mean hardness in Peru as in the United 
States. In a word, then, both the observed very rapid change 
in ionization with altitude and the observed lack of effect of 
latitude upon hardness, or mean penetrating power, especially 
at high altitudes, present rather serious difficulties to any 
sort of an incoming-particle theory even though these particles 
be secondaries which are required to account for but 7 percent 
of the observed sea-level intensity. 

From my point of view, however, these difficulties would 
be reduced and perhaps entirely avoided if it might be 
assumed that the earth's magnetic field does not have the 
generally assumed strength so that secondary particles of 
energies more nearly like those actually observed might be 
mixed, in temperate as well as in polar latitudes, with the 
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incoming photon beam, for the wide distribution of energies 
among these secondary-particle rays would help in under- 
standing the increased ionization with altitude, which 
must be attributed as well to the particle component as to 
the photon component of the rays. If this weaker magnetic 
field cannot be postulated, nor the low penetrability that 
Epstein suggests, then there is no further recourse other 
than to postulate some secondary influence of a non- 
magnetic sort to account for the observed dip in the 
equatorial belt. It is barely conceivable that the following 
hypothesis might have value. 

The earth’s negative electrical field must be taken as 
equivalent to an additional resistance to the inflow of 
negative particles, so that if these predominate, as they 
ought to do in view of Compton encounters, an increase 
in the electrical field would result in a pushing down of the 
curve and a decrease in the field to an elevation of it. Now, 
water vapor, plus the difference in mobility of newly formed 
negatives and positives, is presumably the great underlying 
cause of the electrical separations which result in strong 
electrical fields in the upper air, and it is in the equatorial 
regions that most of the water is lifted into the atmosphere 
and to relatively great altitudes. The absence of an im- 
portant effect at sea level might then be due to the fact 
that the rays that get down to sea level are so hard that 
the influence of the electrical field is negligible in com- 
parison, though not so at the higher levels where the rays 
are as shown relatively very soft. Clay’s strange results 
between Genoa and Suez might possibiy be due to the 
tropical temperatures existing at higher latitudes here than 
in the western hemisphere. It should be possible to test the 
correctness of this suggestion. A few months hence we shall 
have more and better information upon this point. 

The work herein reported has been made possible 
through a grant received for the purpose from the Car- 
negie Corporation of New York administered through the 
Carnegie Institution of Washington, to which bodies | 
wish herewith to make grateful acknowledgment. 
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discoveries in physics may be secured by addressing 
them to this department. Closing dates for this 
department are, for the first issue of the month, the 


twentieth of the preceding month; for the second 
issue, the fifth of the month. The Board of Editors 
does not hold itself responsible for the opinions 
expressed by the correspondents. 


Secondary Absorption and Fine Structure of X-Rays as Observed with a Double Crystal Spectrometer 


Several observers have noted, in the case of solids and 
polyatomic gases, a secondary absorption or fine structure 
extending for some hundreds of volts on the short wave- 
length side of x-ray absorption edges in addition to the 
well-known fine structure in the edge itself. The results 
of experiments on this fine structure are in good agreement 
with Kronig’s' theory of forbidden and permitted zones or 
energy levels in the crystal lattice. According to Kronig’s 
theory, electrons in an atom of the crystal can not absorb 
quanta of certain discrete energy values if the initial 
impulse given to the electron is parallel to a crystal axis. 
On this basis the maxima and minima of the secondary 
absorption should be more pronounced when the electric 
vector of the incident x-rays is parallel to a crystal axis 
than when it is in some other direction or oriented at 
random. It is therefore important to test whether or not 
the secondary absorption is dependent on the degree of 
polarization of the incident x-rays and on the orientation 
of the electric vector with respect to the crystal axes of 
the absorbing screen. 

The double crystal spectrometer with its high resolving 
power should be well adapted for this purpose if a combi- 
nation of crystals can be found such that the first crystal 
can act both as monochromator and polarizer for wave- 
lengths close to an absorption edge of one of the elements 
in the second crystal. Lindsay? has obtained secondary 
fine structure using one crystal for the double purpose of 
grating and absorber. Most of the data on secondary 
absorption have been obtained by the photographic 
method with a single crystal spectrometer, though Davis* 
and Hull‘ have obtained the fine structure in the main 
edge itself with a double crystal ionization spectrometer. 
It is therefore important to have the existence of the 
secondary absorption confirmed by the latter method. 
It so happens that the third order glancing angle with the 
(110) planes of calcite is approximately 45° for the wave- 
length of the K-absorption edge of bromine, thus giving 
complete polarization, while for second order reflection 
the polarization is some 75 percent. 


1R. de L. Kronig, Zeits. f. Physik 70, 317 (1931); 
75, 191 (1932). 

2G. A. Lindsay, Zeits. f. Physik 71, 735 (1931). 

3B. Davis and H. Purks, Phys. Rev. 32, 336 (1928). 

*H. L. Hull, Phys. Rev. 40, 676 (1932). 


One of us (Stephenson) is now engaged in experiments 
with a double crystal ionization spectrometer to test the 
effect of polarization on secondary absorption. The first 
reflection is from the (110) planes of calcite and the second 
from the cleavage planes of an artificial KBr crystal 
kindly sent us by Professor G. A. Lindsay of the University 
of Michigan. Because of the small intensity in the third 
order, the results so far have been obtained with the first 
crystal reflecting in the second order and the second 
crystal in the first, with a cube edge perpendicular to the 
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common plane of incidence. As can be seen from the 
accompanying curve, preliminary results show a _pro- 
nounced secondary fine structure on the short wave-length 
side of the bromine K-edge similar to that obtained 
photographically by others. This curve is a composite 
curve containing points from several different runs in the 
(2, 1) position. The scattering of the points is due chiefly 
to three causes: (1) variation of electrometer sensitivity 
from day to day for which correction is difficult, (2) 
temperature changes between runs, (3) the smallness of 
the current measured, being of the same order as that 
due to a single alpha-particle. The uncertainty of any one 
electrometer reading is quite high, being sometimes 
nearly ten percent of the mean. Several readings are 
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averaged for each point on a curve and the curves are 
repeated several times. All the points are not given in this 
composite curve, but those given are representative. 

The secondary absorption minimum C has been meas- 
ured four times in the (2,1) position and twice in the 
(2, —1) position; E and F have been measured twice in 
the (2,1) position. The peak D may be due to crystal 
imperfections. The fine structure in the main edge at A 
is quite definite. The data so far obtained are insufficient 
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to establish the effect of polarization of the x-rays on the 
fine-structure, but they do confirm the results obtained 
by others photographically. Complete details of the 
experiment will appear at a later date. 
C. D. Cooksey 
S. T. STEPHENSON 
Sloane Physics Laboratory, 
Yale University, 
March, 1933. 


Potential Energy of Diatomic Molecules 


Recently, in a paper by Ikehara and the present writer,' 
calculated results were presented giving the first-order 'S 
energy curves for molecules consisting of two identical 
atoms, each with an outer s-electron, for various values 
of the effective principal quantum number. The *2 energy 
curves were not given there although enough data were 
given so that they could be constructed. If one determines 
the *= curves one finds that for n=4 they are of the usual 
form corresponding to a repulsive force. For n=5, the 
curve has a minimum, however, and if one extrapolates 
to n=6 (which can be done fairly accurately by extra- 
polating the Coulomb and exchange energy curves ex- 
pressed as Morse? functions) one finds that in this case 
the minimum is below the asymptotic value at infinity, 
corresponding to the possibility of existence of a stable 
molecule. This effect is due to the rapid decrease in the 
importance of the exchange energy with increasing 2. 
In Fig. 1 these curves are presented together with the 'X 
curve for n=6. 

The *2 curves in the diagram are not particularly 
important in themselves since the corresponding molecules 
could hardly be obtained under ordinary conditions. 
However, they are of interest in showing that the first- 
order perturbation method may sometimes lead to energy 
curves of the “radioactive decay” type, i.e., having a 
maximum beyond the minimum. It is quite likely that 


1 Rosen and Ikehara, Phys. Rev. 43, 5 (1933). 
2 Morse, Phys. Rev. 34, 57 (1929). 
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for certain electronic levels, molecular energy curves 
similar to those discussed here may occur for smaller 
values of m, and hence with more pronounced maxima; 
the minima might be either above or below the axis. 
Such curves may be of importance in various types of 
predissociation, even if no rotation is present. 

N. Rosen, 

National Research Fellow 

University of Michigan, 
March 30, 1933. 


On the Permeability of Iron at Ultra-Radio Frequencies 


Recently there was published an interesting article by 
Hoag and Jones! who have worked out a method for 
exactly measuring the self-induction of wires in oscillating 
fields of wave-lengths between 22 and 64 cm. By this 
method they determined the self-induction of iron wire, 
calculated the permeability un», and compared it with my 
values of permeability ux, previously calculated? from the 
resistance of wires at a frequency of the same order. 
As was the case for my values of ux Hoag and Jones find 





1 J. Barton Hoag and Haydn Jones, Phys. Rev. 42, 571 
(1932). 

?W. Arkadiew, Phys. Zeits. 14, 561 (1913); Ann. d. 
Physik 58, 105 (1919); Zeits. f. Physik 28, 11 (1924). 


that the permeability uu, calculated by them greatly 
diminishes at shorter wave-lengths. On the other hand, 
their values of permeability are appreciably lower than 
mine. This can be explained in the following manner. 
The essential difference between us and yu, appears in the 
following equations: 


ur=(u2+p!2)! + p! tin = (u?+p"*)! —p’, 


where p’, an essentially positive coefficient, together with 
the magnetic permeability » serves to define the magnetic 
induction B. B is retarded with respect to the magnetic 
field H through the angle 6=arctg p’/u: 


B=pH_ cos (22t/T)+'Ho sin (2xt/T), 
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where T is the period. The value u, determines the coeffi- 
cient of refraction m of electromagnetic waves in a ferro- 
magnetic metal, while u, determines their coefficient of 
absorption? k: 


n=(oTu,)) and k=(o7yx)}. 


In practice ux is usually calculated from the resistance 
of ferromagnetic wires to an alternating current and up, 
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from their self-induction. The values ux, defined by me, 


were calculated from the resistance determined from the 


3W. Arkadiew, Phys. Zeits. 14, 928 (1913). 


coefficient of absorption of waves in parallel wires. In 
another paper I determined the magnetic permeability 
by a method which consisted in comparing the coefficient 
of reflection of the waves from a coarse Hertzian grating 
of magnetic wires, with the coefficient of reflection from 
a grating of non-magnetic wires of the same diameter 
and resistance. 

This method allows one to calculate the value of xz, 
and to compare it with wu, obtained by the method of 
absorption of the waves in wires.’ The results observed 
by me for \ between 1 and 73 cm for a wire 0.0245 mm in 
diameter, are shown by circles in Fig. 1. The values of 
un obtained by Hoag and Jones are shown by crosses. 
They relate to a wire 3.18 mm in diameter. The con- 
cordance of the results of Hoag and Jones with my values 
is better than expected. We see from the above, that the 
values obtained by them for u, were bound to be smaller 
than my values for ux. We see also, that the permeability 
of iron really decreases, approaches 1 at A=1 cm, and 
approaches the value corresponding to thermal and 
optical frequencies. 

W. ARKADIEW 

J. Cl. Maxwell Laboratory, 

(Moscow Magnetic Laboratory), 
University, 
January 21, 1933. 


4W. Arkadiew, Ann. d. Physik 45, 133 (1914). 

®5W. Arkadiew, Ann. d. Physik 81, 661 (1926); Table 
IV, Fig. 3; Festschrift d. Mosk. Magn. Laborat. p. 53 
(1931); Ann. d. Physik 11, 406 (1931), Fig. 3. 


An Attempt to Observe a Helium Isotope 


Langer' has recently demonstrated that Be® is radio- 
active and has suggested the possibility that He® might 
be one of the products of the disintegration and that, if so, 
He® should be found in old beryllium mineral. 

Professor Gruner placed at our disposal a large crystal 
of beryl from the pre-Cambrian dykes of the Black Hills 
in South Dakota (estimated age greater than 3108 
years). The gases were driven out by heating the crushed 
crystal and were passed through a charcoal trap immersed 


1R. M. Langer and R. W. Raitt, Phys. Rev. 43, 585 
(1933). 


in liquid air. A mass-spectrograph analysis of the gas not 
absorbed showed He‘. If He® was present its concentration 
was less than one part in 40,000 relative to Het. 

As a check on the apparatus, ordinary hydrogen was 
examined for the (H*H')* molecular ion. An analysis of 
the data, taken at different pressures to correct for (H;')*, 
showed H? present in a concentration of the order of 
1 part in 30,000 relative to H!. 

Joun T. TATE 
Puitip T. SMITH 
University of Minnesota, 
April 4, 1933. 


Collisions of Neutrons with Protons 


With the hope of learning more about the nature of the 
neutron I have been studying the collisions between 
neutrons and protons using a Wilson cloud chamber. 
This work is being held up to permit the completion of 
some measurements on nitrogen disintegration tracks. As 
the progress of the two problems is very slow because of 
the weakness of my neutron source I have thought it 
advisable to make this preliminary report on the neutron- 
proton collisions. 

Although the neutron itself leaves no track in a. cloud 
chamber the essential facts of the collision are known if 


one measures the direction in which the proton is projected 
relative to the initial direction of motion of the neutron. 
I have therefore arranged a source of neutrons in a small 
volume surrounded by a cylinder of paraffin of diameter 
12.5 mm, the whole being placed at the center of the cloud 
chamber. The tracks formed by the protons projected from 
the paraffin are photographed. Since the neutrons proceed 
in a radial direction the angle between the directions of 
motion of the neutron and of the proton may be measured. 

About 150 such tracks have been measured. Corrections 
are applied for the volume of the paraffin effective in the 





eee 














“Noe 














LETTERS TO 


scattering and for the solid angle subtended by the 
chamber. To eliminate the possibility of measuring tracks 
due to radioactive contamination (such tracks are about 
as abundant as the proton tracks, but much shorter) only 
those which extend to the walls of the chamber are con- 
sidered. On plotting the number of protons per unit solid 
angle projec.ed between @—5° and @+5° against @ one 
finds a distribution curve showing a decided maximum in 
the forward direction and falling toward zero at 90°. 

The behavior of two different models of the neutron in 
such collisions has been discussed theoretically by Massey.' 
He finds that both a dipole and a spherically symmetrical 


1H. S. W. Massey, Proc. Roy. Soc. A138, 460 (1932). 
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model should exhibit peaks around 90°, which is in dis- 
accord with these observations. The present data seem 
best to fit the cosine curve which describes the scattering 
of elastic spheres. Massey's models seem to exhaust the 
proton-electron combinations so that one is tempted to 
regard the neutron not as a composite particle but as an 
elementary particle. 
The details of this work will be described later in a 
paper to be sent to the Physical Review. 
F. N. D. Kurie 
Sloane Physics Laboratory, 
Yale University, 
April 4, 1933. 
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